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Die Auswélbungen der Kreisringmembran unter hydrostatischem Druck * 


Zweite Mitteilung: Inhomogener Spannungszustand** 


Von K. Karas 


1, Einleitung. Wahrend in der ersten Mitteilung! sowohl bei der Kreis- als auch der Kreisring- 
membran ein homogener Spannungszustand vorausgesetzt wurde, bei dem also fiir alle Schnitt- 
richtungen senkrecht zur Membran und fir alle Membranpunkte die gleiche SpannungsgréBe S 
vorhanden ist, soll hier ein wohl achsensymmetrischer, aber inhomogener Ausgangsspannungs- 
zustand vorausgesetzt werden, dessen HauptschnittgréBen mit S, fiir Schnitte senkrecht zum 
Halbmesser r und mit S, fiir Schnitte in r bezeichnet werden mégen. Wie in der fritheren Arbeit 
soll auch hier der AuBenrand der Kreisringmembran als fest angenommen werden, wahrend der 
Innenrand derselben als fest oder vollkommen frei beweglich vorausgesetzt wird, wodurch sogleich 
die auch in der ersten Arbeit berechneten, praktisch wichtigen Grenzfalle mit erfaBt werden kénnen. 
Sie konnten als Spezialfalle des allgemeineren Falles elastisch gelagerter Innenkreisplatte erwiesen 
werden, wobei elastische Lagerung wieder sowohl gegen Verschieben, als auch gegen Verdrehen 
derselben gefordert wurde, um méglichst allgemeine Verhaltnisse zu erfassen. Auch wurden fiir die 
beiden erwahnten Grenzfalle in zwei Zahlenbeispielen die Héhenlinien der Membranwélbflachen, 
sowie die Volumina ihrer Wélbungshiigel berechnet und die Héhenlinien und Meridianschnitte 
wieder in anschaulichen Schichtenplanen dargestellt. Alle Ergebnisse, die sich wiederum in ge- 
schlossener Form darbieten, konnten fiir den Fall eines homogenen Ausgangsspannungszustandes 
durch die der ersten Arbeit verifiziert werden. 


2. Der inhomogene, achsensymmetrische Spannungszustand. Bei ihm geniigen bekanntlich? 
die auf die konstante Membrandicke bezogenen SchnittgréBen S, und S, den Gleichungen 


Cc C 
S,=B—= S,=B+ 3 (La, b) 


roe, 


und die beiden Konstanten B und C werden durch die innerhalb gewisser Grenzen willkirlichen 
Werte von S, am AuBen- bzw. Innenrand der Membran bestimmt; diese sind 


fore TO 8. SVU) Ble) ruled aga (2a, b) 
Mit (2a, b) folgt® aus (la) 


Sar? — S;r? te} (ifs Te 
B = oo —, C == (S, ae S;) r2 mage ° (3a, b) 


Spaterhin wird sich die Einfihrung der abkiirzenden Konstanten 


Sa—Si) rar? 
te ta eae Tall (4) 
B Sgr2— Sir? 


* Meinem verehrten Kollegen, Herrn Prof. Dr. Alwin Walther, zum 60. Geburtstag gewidmet. 

** Der Verf. wurde bei der Zeichnung der Abbildungen und bei der Berechnung der Tabellen von Herrn 
Dipl.-Ing. K. E. Meier-Dérnberg in dankenswerter Weise unterstutzt. 

1K. Karas, Ing.-Arch. 25 (1957) S. 359. Im folgenden karz mit I zitiert. 

2 Man vgl. etwa A. Foppl, Vorlesungen uber Techn. Mechanik, Bd. III, 14. Aufl., S. 331—335. 

3 Die bekannte Formel og = pr/6 fiir die Normalspannung langs einer Erzeugenden eines Zylinders, der 
einem Innendruck p unterworfen ist, kann leicht aus (1b) und (3a, b) fur kleine Wandstarken 06 abgeleitet 
werden. Mit rj =1, ra =: + 0, Si = — p06 und Sq = 0 erhalt man 


, . pr iat ae : 
ri—r?=2ro, Le rior vay 5? C= pes =P a3 


damit ergibt (1b): Sp = pr, somit 09 = = Ps . Man vel. auch M. M. Filonenko-Boroditsch, Festig- 
keitslehre, Bd. 2, S. 76—82. 
11 
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als niitzlich erweisen. Neben den radialen SchnittgréBen S, und S$; an den Membranrandern sind 
noch die tangentialen SchnittgréBen S, und S,, ebendort von Bedeutung. Man erhalt + aus (1b) 


mit 3a, b) und fiir die Randwerte r, bzw. r; von r wie in (2a, b) 
oe Sa ("a +H) 258i fe 29672 — si (ra + ae 


p yao ? Q; Oo amareo 
a Ue Ls U Ya Up 


(5a, b) 


Wahrend bei der Scheibe und dem dickwandigen Hohlzylinder S, und S; in gewissen, durch die 
Festigkeit des Materials vorgeschriebenen Grenzen beliebig gewahlt werden kénnen, ist dies bei 
Membranen, die nur Zugspannungen iibertragen kénnen, nicht ohne weiteres méglich. Man wird 
also neben positiven Randwerten fiir S,, namlich fiir 


Su 0 Spe 0 (6a, b) 


v 


(S, kann in der Grenze auch verschwinden) auch weiterhin fiir alle Membranpunkte positive 
Werte der HauptschnittgréBen fordern miissen: 

r2 r? ar? 
Set? Sir? = Ges) Sy ri =5i7? (Se Shae 
s, = ————_ = ieee i) Ss, =—— — 


r2 
aa een er (esl), 
fine = 1 7, 

S,und S, sind nach (7a, b) monoton ab- oder zunehmende Funktionen. [Man vegl. auch Abb. 2.] 
Daher ist zunachst wegen (6a, b) auch (7a) verbiirgt. Auch (7b) ist erfiillt, wenn in (5a, b) So, > 9 


und S,. > 0 befunden wird. Mit der Abkirzung 


p=t>l (8) 
folgt daher aus (5a, b) 
2r? r2 | r? 
(S,; kann in der Grenze auch verschwinden). Mit (8) soll also 
2 wet 1 1 
ee Soy cy atuelimen ee ah 
8.2 Sie) Se = Si(z + za) (9a, b) 
sein. Man erkennt leicht, da® die Erfiillung von (9b) auch die von (9a) bedingt. Es ist namlich 
pel seis Vee ye (eae 0 


ue ma pe i eae 2 pu? (ue? ap 1) =) pe (we + 1) 
Ist also (9b) erfiillt, so weisen in allen Membranpunkten die SchnittgréBen aller Schnittrichtungen 
nur positive Normalspannungen auf. AuBer den Beschrinkungen (6a, b) und (9b) unterliegt die 
Wahl von S, und S; keinen weiteren Bedingungen?. Auf jeden Fall befindet sich ein Membran- 
sektor (Abb. 1) im Gleichgewichte gegen Verschieben in der r-Richtung; denn die in der y-Richtung 
ist wegen der vorausgesetzten Achsensymmetrie des Ausgangsspannungszustandes ohnehin erfiillt. 
Das Gleichgewicht in der r-Richtung erfordert (nach Weglassung des Faktors dq) 


S.ta—f Sal 5, hee : (10) 
WV ecerie || oe meee meres eemich eae (7h 
g =, ergibt sich aus (7b) 
: Sa 2 Sr? i a — Si att 
| Sear= (s Te S r?)+(S 92) ane ' 
Tate Ti 


Ts 
v 


und damit erkennt man die Erfiillung von (10). 

Ein praktisch manchmal auftretender Grenzfall eines inhomogenen statischen Spannungs- 
zustandes einer Kreisringmembran liegt dann vor, wenn ihr Innenrand frei ist. Mit S. = 0 und 
2 ees Ps 7 


7 : i DAs 
Auch die maximale Schubschnittgr6Be Tinax = 1/2 (Sp; — Si) am Innenrande r = r; kann manchmal 
interessieren; man erhalt mit d als Membrandicke 
(Sa — Sj) r? 
Dinax = d+ Tmax = ——* 
Chall 


2 Es sei ausdriicklich erwahnt, da®B nach (9) Sa = S; sein kann. 


rg 
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S, > 0 ist dann auch (9a, b) erfiillt, und man erhalt aus (7a, b) speziell: 


lie r? 
ees erat ae Sp =8 


t 


f 

2 ( + =| ; (11a, b) 
SchlieBlich erweist man, daB wegen (4) im ganzen Intervall (7) von r die Funktion 

a? Sar (r?—r8) + Sir? (2 — 79) 


r? ; (Sa r 22S: r?) ro > 9 (12) 


ist, sobald (9b) erfiillt wird. Denn der Zahler in (12) ist wegen (6a, b) stets positiv und fiir den 
Klammerfaktor des Nenners folgt wegen (9b) 


7 a i rts 
S,2—S,2 = S,( —— —r)=s," E> OE 
Statt (12) kann man auch 
ear (12a) 
schreiben. 

Hingegen kann im selben Intervall 
und unter den gleichen Voraussetzun- 
gen die Funktion 

r2 
et oat 
Sa t2 (r? —r2) + S; r?\(r2 —r?) 
(Sa 8) r2 r? a 
=0 (13) 
sein, je nachdem, ob wegen 


24 72 
s,— 8,2 8,55" =1) 


2 
27s 


L 
2 
Zr 


2 
. 
a a 
i 


<0 (13a) 


die Differenz (S,—S;) positiv oder 
negativ ist!; dann bestehen die Un- 
gleichungen 


Sal = De = S; 


oder (S,—S,) > 0, 
2 
Scull 2A S; = Sy = S;* a 
2 pb 
oder (S,—S,) <0. 


S., ist dabei durch die Festigkeit 
des Membranmaterials bedingt. Im 
obigen Grenzfall S; = 0 ist wegen (6) 
sicher S,— S; > 0 und daher nach 
(13) 1—r?/a® < 0; dies erkennt man 
auch daran, da dann (4) a? = r? er- 
gibt und also 1 —r?/a? = (r? — r°)/r? 
< Oist. Im andern Grenzfall gilt in(9b) 
das Gleichheitszeichen, wozu (5b) Sy, = 0 ergibt; dann aber ist nach (13a) S,— S; < 0 und somit 
nach (13) 1— r?/a2 > 0; dies bestatigt. man auch wieder durch (4), das dann a? = — rj ergibt, 
wonach 1 — r2/a? = (r? + r?)/r? > 0 ist. In beiden Grenzfallen ist natiirlich mit a? = + r? auch 
(12a) und somit (12) erfiillt. Im ersten Grenzfall verschwindet somit die Radialspannung S;, im 
zweiten die tangentiale Spannung Sy; am Innenrandkreise. 

Liegt schlieBlich als wichtiger Spezialfall ein homogener Spannungszustand vor, so ist 5,=S,=S 
und somit nach (la, b) B = S und C = 0, was auch aus (3a, b) folgt. Auch aus (5a, b) ergibt sich 
dann So, = Sy; = S, und zu demselben Ergebnis fithren auch die Gleichungen (7a, b). Nach (4) 


ist im Falle eines homogenen Spannungszustandes stets a = 0. 


Abb. 1. Ein Membranelement mit den an ihm angreifenden Schnittkraften. 


1 Man vel. FuBnote 2 von S, 158. 
Ulf 
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Fir alle zulassigen Membran-Spannungszustande ist somit der Inhomogenitatsfaktor a auf das 


foleende Intervall beschrankt: 
i‘ —rfid=4+F. (13 b) 


3. Die Grundgleichungen des Auswélbproblems. Diese sollen unter Betrachtung des Gleich- 
gewichtes eines Membranelementes zuniachst fiir irgend einen allgemeinen Ausgangsspannungs- 
zustand S,(r,~) und S,(r, p) abgeleitet und erst nachtraglich fiir den in Ziff. 2 angegebenen achsen- 
symmetrischen Spannungszustand spezialisiert werden. Mee 

Betrachtet man das Gleichgewicht des in Abb. 1 dargestellten Membranelementes, das in Rich- 
tung w einem Flachendruck p(r, ~) ausgesetzt ist, so erhalt man 


C OS; 5 ED oe obo 
— S, sin ea rdy + (S, S os ar) sin (Fe -F — dr) (r + dr) dp —S, sin ie %) dr 


1 dw 1 Ow 


OS, : i 
ati (s, +. a dy) sin( ie | roa dg) dr + prdpdr=0. 


Der Winkel po os der SchnittgréBe S, gegen die unverzerrte Membranebene vergréBert sich nam- 


on rag = Nanay dp. Ersetzt man die sin-Funktionen 
r oy r 0g? 


lich beim Fortschreiten um r dg um Pep ( 
durch die Argumente und multipliziert aus, so erhalt man nach Weglassung der sich tilgenden 
Glieder und nach Streichung des dann allen Termen gemeinsamen Faktors dr dp 


ow Ow OS; Ow Sy 02w CSoulanew 
Ser ee te ees Pag Oe Oe (14) 


r” Or? TOR r dg? 


Wie man leicht erkennt, gestattet (14) eine Umformung in die selbstadjungierte Schreibweise 


0 7) Ib Ow * 
ar (Se ar) r Op (S, ap | eas ‘ ee 


Die partiellen Differentialgleichungen (14) bzw. (14a) beschreiben die Auslenkung w(r,q@) fir 
irgendeinen, der Membran auferlegten Spannungszustand S,, S, und fiir irgendeinen Flachen- 
druck p(r, ¢). 

Liegt speziell der durch die Gleichungen (la, b) erklarte Lamésche Spannungszustand vor, so 
ergibt die Einfiihrung von (la, b) in (14) nach einigen Umformungen und Teilung durch r 


C\/ 1 ow 1 ?w CNG 
[2 . alle gas 25%) | [8 =) oe 5 ie Se 
Der Term mit 0S,/ép in (14) fallt namlich wegen (1b) fort. 


Gleichung (15) liegt nun allen weiteren Berechnungen zugrunde. Der hier durchweg voraus- 
gesetzte hydrostatische Druck p wird auch hier wieder gema der Gleichung 


P=p +p =ykh+yrcos@ (16) 
in einen konstanten Teil p’, der gleich ist dem Flissigkeitsdruck y h im Membranmittelpunkt 0 
(I. Abb. 1) und einen mit der Tiefenlage unter 0 um y = r cos @ linear veranderlichen Teil p'’ auf- 
gespalten, wie das in I bereits geschehen ist. Hierbei bedeutet y das spezifische Gewicht der 
driickenden Flissigkeit, h die Tiefe des Membranmittelpunktes 0 unter dem Spiegel, und r, ~ sind 
die Polarkoordinaten eines Membranpunktes, wobei der Winkel g vom vertikalen Membrandurch- 
messer (y-Achse) an gezahlt wird. Der Druck p’ bedingt achsensymmetrische, d. h. nur von r ab- 
hangige Auslenkungen w’, wahrend p’’ bzgl. der horizontalen x-Achse antimetrische Auslenkungen 
w'' bedingt, die auBer von r auch von ~ abhangig sind. Wegen der Linearitat der Differential- 
gleichung (15) ergibt sich die Gesamtauslenkung w durch Uberlagerung der Teilauslenkungen w’ 
und w’’ (die Striche bedeuten also hier keine Ableitungen): 


w(r,p) = w'(r) + w'"(r,@) . (17) 


* Ersetzt man in (15) B durch P, ferner C durch —c und p durch — ¢ 6%w/ét? mit ¢ als Flachendichte 
(Reinstein S. 9) und p durch p (Reinstein FuBnote 1 von Seite 102), so erhalt man 


Cy Ucw 1 w ON GHD Pw 
P Ei | eam A\\ | 4 
( salle or | r apt] [P 7) are” aie 
Diese Gleichung ist auf anderem Wege in der bemerkenswerten Géttinger Dissertation als Gl. (206) S. 124 


abgeleitet von FE. Reinstein, Untersuchungen iiber die Transversalschwingungen der gleichformig gespannten, 
elliptisch oder kreisformig begrenzten Vollmembran... Gottingen, 1911. 
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a) Ermittlung der achsensymmetrischen, Auslenkung w'. Wegen 


Ow’ Ow’ 
ap aga —° 
erhalt man aus (15) nach entsprechender Umformung analog (14a) 
d C\ dw 
ap es | yh (ior (18) 


Tn (18) 1aBt sich ersichtlich eine erste Integration leicht durchfiihren, und man erhilt mit D, als 
erster Integrationskonstanten 


C\ dw" 1 
Ceara eae uae oe (19) 
bzw. mit (4) 
a”\ dw’ 1 A 
Br | 7 5 yrs) = 0. (19a) 


In (19) kann die Trennung der veranderlichen r me w’ leicht vorgenommen werden und ergibt 
1 2 h r dr 


dw += rdr + Gaus: bh + D i) Br ion 


Mit D, als zweiter ieee erhalt man hieraus 
1 G D 
AT adh act ay We ny she 1 ee 
w cee t+ (grhge ty g)in (By C 
Man kann (20) unter Ausniitzung von (4) und mit etwas anderer Bedeutung der Konstanten D, 
auch folgende Gestalt geben: 


w= sags lage tre top) 


2=0. (20) 


+D,=0. (20a) 


In den folgenden Beispielen wird die Grundgleichung (20) mehrfach ausgenutzt und die Berechnung 
der Integrationskonstanten D, und D, durch die verschiedenen Bedingungen an den Randern der 
Kreisringmembran gezeigt werden. 


b) Ermittlung der antimetrischen Auslenkung w’’. Der Druckanteil p’’ = yrcos@ 
in (16) ergibt in (14a) gemeinsam mit (la, b) eingefitihrt zunichst, wenn man jetzt w'’ statt w 
schreibt, 
0 C\ Ow’’ 1 ON 
ie ; erase es 
. (Br : some | as [B + =| Wig tne cos = 0% (21) 


Da in (21) nur die zweite Ableitung von w’’ nach @ auftritt, so kann man diese Gleichung durch 
den Ansatz 

w'' = R(r) cos (22) 
befriedigen, in dem R(r) nur eine Funktion von r allein ist. Man erhalt dann nach Weglassung des 
Faktors cos y fiir R die totale Differentialgleichung 


# (0-2) 8) (Br +-£)8 yo. a 


Nach Ausfithrung der Differentiation erhalt man aus (23) 


C\@R C\ (dR R ee 
[2B r aed aa [2B -f- 7) (s *) on ec (23 a) 
Das allgemeine Integral von (23a) setzt sich aus dem der homogenen Gleichung 
(Br ae -(B | Ne RY 0 (24) 
pl) Ghee Mochi i 


und aus einem partikularen Integral zusammen, das aus dem ersteren durch Variation der Kon- 
stanten erhalten werden kann. Man erkennt, daB 
R=E,r (25) 


mit FE, als Konstanten die Gleichung (24) befriedigt, somit ein partikuladres Integral derselben 


darstellt. Mit dem weiteren Ansatz } 
dR dz d?R dz dz R te (26) 


[ia OG eee Ent 7a AG Er 


1 Man vgl. hierzu etwa R. Rothe, Héhere Mathematik Teil III, 4. Aufl., S. 189 f. 
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erhalt man aus (24), wenn man die Integrationskonstante E, sogleich weglaBt, 


d ad? C d 
[Br ale _ | Fi) | (2 | =)? _ = (0. (27) 
Mit 
dz 
HAUG) lia ar (28) 


erhalt man aus (27) nach einigen Umformungen und Trennung der Veranderlichen p und r die 
Differentialgleichung 


dp dr 2 Cae are! 
P aera 5 7 ser ips Se oa 
deren Integration mittels Partialbruchzerlegung sofort das Ergebnis 
prin) = E, (30) 


mit E, als neuer Integrationskonstante liefert. Wegen der nach (28) nochmals notwendigen Inte- 
gration von (30), die leicht mittels Partialbruchzerlegung durchgefiihrt werden kann, erhalt man 


E 2 
p= “ini 9 ee (31) 


wobei die abermals auftretende Integrationskonstante in (31) sofort mit E, bezeichnet werden 
konnte, da (31) im Zusammenhalte mit (26) zu 


R, = Brin /1—% + E,r (32) 


fiihrt, so daB also im zweiten partikularen Integral in (32) wieder (25) gewonnen worden ist. R,, 
in (32) stellt somit das allgemeine Integral der homogenen Gleichung (24) dar. Fiir die inhomogene 
Gleichung (23 a) kann man mit (4) schreiben! 


a a 
PRE lee ene ee ae 

le cee (33) 
dr? ry a dr r2 a? a> 


Fat man nun die GréBen E, und E, in (32) im Sinne der Lagrangeschen Methode der Variation 
der Konstanten als Funktionen von r auf, so erhalt man fiir E; und E; (Striche bedeuten Ableitung 
nach r) die inhomogenen Gleichungen 


a 


1 a 1 
In {1 a ee ee 
9 u2 n 7) re a Bibel 
r Tie 


mit der nach (12) fiir r > 0 nicht verschwindenden Wronskischen oder Systemdeterminante 


E,(r) 


i 
r fe 
A 1 i a r2 ees a2 
2 
und den Lésungen 
t 2 t 


Auch die Gleichungen (34a, b) lassen sich geschlossen integrieren; man erhalt nach langeren 
Rechnungen und unter Weglassung der neu auftretenden Konstanten, die im allgemeinen Integral 


der inhomogenen Differentialgleichung (33) in den Termen von (32) nach erfolgter Zusammen- 
fassung bereits enthalten sind, 


2 2 
E,(r) = 3 f In(1 a ee r? — q4 He 1 


| EG) = = Ose 


DTG aRAP Faktor der héchsten Ableit IR — ividi 
PERS AEE eitung zu erhalten, wurde durch (Br — C/r) dividiert. Man vel. 
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Setzt man die Funktionen (35a, b) fiir E, und E, in (32) ein, so erhalt man das gesuchte partikulare 
Integral R,, von (33) nach pateicebienier Kievene zu 


2 
eocee vhs 2 | aie 
- Sap 1 erin me 


| (36) 


Durch Zusammenfassung von (32) und (36) erhalt man das allgemeine Integral der inhomogenen 
Differentialgleichung (33) zu 

2 2 

R=R,+R,=E,r4 Ey 5 2 ln eal 5a +a? rIn/1 —5 


2 


r2 


(37) 


SchlieBlich erhalt man die Auslenkung w'' mit (37) aus (22) zu 


oa 


2 ne 
= [Fr ye pial —Z] sale | trial —Z)) 


cos @. (38) 


Fitr homogenen Spannungszustand wird B = S und a = 0, wobei der letzte Term in (38) ver- 
schwindet, wahrend der Faktor von EF, iibergeht in 


lim lr In (1 —s||= =limr aoe Soa ene qi : 


a—>0 a>0 ow. r 5. pre — at) Bip © 


(39) 
Somit erhalt man aus (38) 


E, y 
Ne 3, 8B pcos p 


2 


Hiermit ist aber mit etwas anderer Bezeichnung der Konstanten die Gleichung I (12) wieder 
gewonnen worden. 


4. Die Kreisringmembran mit fester Innenkreisplatte. Fiir den achsensymmetrischen Fall: 


gilt (20a). Zur Bestimmung der Integrationskonstanten D, und D, dienen die Randbedingungen 
(I, Abb. 3) 


WOT i.) oe (P= 4) = 0 (40a, b) 
Gleichung (20a) liefert mit (40a, b) die Lésungen 
2 
r2 In | soll r?In “ il | 
1 ‘i =P lyh {@ | 
D; ee 5 7, h a? aa Sti Sat 5 D, = eee a . a, “Fa =O (41a, b) 
mae "a 
Mit (41a, b) erhalt man aus (20a) fiir die achsensymmetrische Auslenkung w’ 
/r2 — a? Ve —a@ , r2— a2 
l whee is eres ae rin orm vai In ee 42 
w =—* ee = : - (42) 


4 B r2— @? 
i 
r?— a? 

L 


Dies erfiillt sowohl (18) als auch die Randbedingungen (40a, b). Die Wurzeln kénnten natiirlich 
auch weggelassen werden. 


Fiir homogenen Spannungszustand erhalt man aus (42) mit B = S und a = 0 


Yo r Ya 
—r In — + r?In— +77 ln 
if ri 


somit wieder I, (27). 


Fiir das Volumen des kreisringformigen Wélbungshiigels erhalt man 


Tan 


V=2f fw! rdrdp=22f w'rdr, ($8) 
mY Pe 


worin w’ aus (42) einzufiihren ist. Mittels der Formeln 


: In x 
[ried agen —F 46, [ otineae = oe + C (44a, b) 
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erhalt! man nach Einfiithrung der Substitutionen 


2 2 2 2 
| Rese | r a 
= x bzw. S — ie (45 a, b) 
rw ra 


die folgenden bestimmten Integrale als Teilergebnisse von (43) : 


x 
a 
—o-s 


Mit (46a, b) erhalt man aus (43) nach einigen Umformungen 


ite Lene te ta we te 2 j2=5 a a| AT 
Vi=wg Boe pea We ae Gy) 
im Vz — a? 
Vi yh r2—r? 1 (2 Sere eae ae rey? fe 
Writtel r2 7 4B re EB wae 2 | e 2 a” } In 72 — a M ( 4) 
In Ve = 
ig a 


Aus Griinden des Vergleichs mit anderen Fallen ist dabei w,,,;,,, wie in I auf die Au®enkreisflache 
bezogen. 
Fiir homogenen Spannungszustand folgt aus (47) mit B = S und a = 0 
we Ana A en a 
V ois: oer [tr + r?) In 


% 


somit wieder das Ergebnis I (40). 

Fiir den antimetrischen Auslenkungsanteil w’’ hat man (38) wieder den Randbedingungen 
(40a, b) mit w’’ statt w’ zu unterwerfen und erhalt nach langeren Rechnungen fiir die beiden 
Integrationskonstanten E, und EF, mit der abkiirzenden Bezeichnung 


a2 
—— 
ie r? (r2 — a?) r2— a? 
ake a aa aera uty a aes eee (48) 
as 


die folgenden Werte: 


y a 
b= SEN ie Im ( — | In 


e u | | 
—riin(L—F) +rtin(t FL 
L a 
3 (49a, b) 
ve y a? 2 a y a? ith a 
hk, ABN r2—r? + a? In a; || =aBN & r? + a In 2 Al 
[> wae! 


Mit (49a, b) erhalt man aus (38) nach langeren Rechnungen 


2 Tee ae 2 ee 
"= seR{—n(L— 3) [tn gaa t 2d] + n(1—F) [ot nS — 2 — ot] 
2 2_ 2 
a In(1 —f5) G ne | (r2 r)||r c08 (50) 


Dies erfiillt sowohl die Differentialgleichung (21), als auch die beiden Randbedingungen (40a, b) 


. . p 
wie man sich durch Ausrechnung tiberzeugen kann. 


1 Es ist 


SO arnt eh 
| moa Sra = 3 ym ee. 


2 @ 1 aoe el 
| mFS rar Fe atym Ae a) ce 
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Liegt homogener Spannungszustand vor, so ist B = S und a — 0. Hierbei verschwinden in (50) 
sofort alle Terme mit a?, wahrend die iibrigen mit N in (48) unbestimmte Formen ergeben. Man 


erhalt durch Grenziibergang ahnlich wie in (39) 


—2a —2a —2a 
2 2 2 2 2 2 
r r2 : r r?) ——____ | (r ia) = = 
: eas COs ( sat (a a ba Dies wage 
lim w’’ = lim a Se AGE a ee Nee ee 
ptt. % ee 68S (r? — a?) (— 2 a) + (72 —a?)- 2a 
(2 — a2) — a) 
2 2 aTe Eco potsan eo 
daca rei han th ete 
2 2 278 
ES ON a ee 22 ee ee Pe 
8S it : 1 a 8S ie ml L r2 
r2 r2. 


Damit aber ist Gleichung (32) der Arbeit (1) wiedergewonnen. 
SchhieBlich erhalt man durch Addition der Gleichungen (42) und (50) mit N in (48) 


20a a Sy i 


eay OE ces 2 
3 aa In(1 =) a? In 5 - | (r? rp| |. In(1 aH a In . a (r?2 — r?) 
2 r2 — q? 
le In(1 3 a In eae a (=r) 7 008 p . (51) 


Erteilt man w den konstanten Wert w,, so stellt (51) die Gleichung der Héhenschichtenlinien bei 
inhomogenem Ausgangsspannungszustand in Polarkoordinaten dar. Ein Zahlenbeispiel hierzu ist 


in Ziff. 7 berechnet. 


5. Die Kreisringmembran bei freier Innenkreisplatte. Es gelten fiir den achsensymmetrischen 
Fall wieder die Gleichungen (19a) und (20a). Die beiden Integrationskonstanten D, und D, ergeben 
sich hier aus dem Gleichgewicht der Innenkreisplatte gegen Verschieben und aus der Forderung 


des Verschwindens der Auslenkung w langs des AuBenrandes r = r, (man vgl. I Abb. 6a) 


Bile PT Te 9 hte Wt = T= 0 (52a, b) 
i ry Lek 
(52a) liefert mit (19a) D, = 0 und damit ergibt dann (52b) mit (20a) sofort D,: 
lyh lie 
Die ODS rp (iteln|s 1). (53a, b) 
Hiermit liefert (20a) schlieBlich wegen (13b) die Lésung 
ee yer 
wo = : aes G r2) + a? In = j (54) 


Fir homogenen Spannungszustand folgt aus (54) mit B = S und a = 0 sofort wieder I (44a). 
Die Verschiebung s der Innenkreisplatte folgt aus (54) mit r =r; zu 


ee 
aes Be tr r?) + @ In oa : (55) 


per Sie es 
die fiir homogenen Spannungszustand sofort wieder in I (46) iibergeht. 


Fiir das Volumen V des Wolbhiigels erhalt man mit (43) und der Formel 
eeihe: ee ees bs (56) 
x 


nach Hinzufiigung des Zylindervolumens 1? s unter der ausgelenkten Innenkreisplatte, wobei s 
aus (55) zu entnehmen ist, nach langeren Rechnungen' 


4 4 2__ g2 
Ve eee la : rt 4 @2 (r2 — 7?) 4 at In “2 F (57) 


1 Es kann auch die zweite Formel der FuBnote 1 von S. 164 hier ausgenutzt werden. 
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Fiir homogenen Spannungszustand folgt aus (57) mit B= S und a = 0 sofort wieder das 
Ergebnis I (47). Teilt man schlieSlich (57) durch rj 2, so erhalt man fiir die mittlere Auslenkung 
4 a2 


== 7h ha 
A beet re ree iy, (58) 


(oy SS 
W mittel 4B r2 2 a2 


Der antimetrische Auslenkungsanteil w'’ wird wieder durch (38) beschrieben. Bedeutet nun 
5 den Winkel, welchen die Ebene der Innenkreisplatte nach erfolgter Auslenkung mit ihrer ur- 
spriinglichen Lage einschlieBt, so bestimmen sich die beiden Integrationskonstanten E, und E,, 
sowie 6 aus den Forderungen des Gleichgewichts der Innenkreisplatte gegen Verdrehen, und der 
Anpassung von w” fiir r = r; an die Kontur der um den Winkel 6 verdrehten Platte und fiir r=r, 


an die Bedingung analog (52b), namlich w’’ = 0 (I Abb. 6b): 


Mom p" — Moms” — Moms y = 9, w''(r,)==0r, cosy. 0 '(r:) == 2 (59a, bay 


In (59) bedeutet M,,, p> das Moment der Driicke p’ = yr cos p bagl. der x-Achse und kann direkt 


aus I entnommen werden [Gl. (50a)]: 


1 
Mom p'’ ge pe (60a) 


S’ bedeutet die Komponente von S, am Rande der Innenkreisplatte, also von S; senkrecht zur 
u 


(x,y)-Ebene, also senkrecht zur Ebene der Innenplatte in ihrer urspriinglichen Lage. Nun folgt 
aber aus (38) 


3 In(1 =) an | y 3 eatin 1 — Bo Ta} eos. (38a) 


r? — a? 8B a? | a? — r? 
if 


Mit (38a) erhalt man schlieBlich als Moment der erwahnten Komponente S7,, die den Hebelsarm 
r, cos ~ besitzt 


n/2 a1/2 


a) i ho C Ow” 
Mom sy. = 4 | S; a ns \ a dy YT; cos PP —= 4. / [2 el or) r? cos Pp 
0 


j 
C E, a Bie" 
= rin [B 4 )| E, oa ln ( =) = er 


L 


4 | 2 | 2 a? r2 
| (3 ba? n 17 soe ar (60b) 
Die zur (x, y)-Ebene parallele noch verblicebene Komponente Sr, wird nun ihrerseits gemaB I Abb. 7 
wieder in zwei Komponenten Sr und Sry zerlegt, von denen nur die letztere ein Moment bzgl. 


der (unverschobenen) x-Achse besitzt. Da der Hebelarm von Sry = S; cos ¢p gleich ist r; cos ¢y 6, 
so erhalt man 


n[2 


Moms). = 4(B— 70 T; | r, cos? ~ dp =2(B—7\r80 = 7 B (r? — a*) 6. (60c) 
j L 


In (60c) hat in der zuerst vorgeschriebenen Gestalt der Klammerausdruck die Bedeutung S, 
Fithrt man nun (60a,b,c) in (59a) ein, so erhalt man nach Kiirzung durch 7? und Teilung durch B 


Ih ore a E. a2 2 a2 
4 B ( #) {7 Bae a(l | | oa 


i 


zy 2 i: 2 ar? f 
+ apart tein em ey + gaa} (2S (61a) 
Die Bedingung (59b) ergibt mit (38) nach Kiirzung durch r, cos y 
1 di 2 ) 2 | 
ES ln (1 7) p(t | a@ In alls (61b) 


SchlieBlich ergibt die Bedingung (59c) mit (38) nach Kirzung durch r, cos p 


E 2 
E, 4 yaln (I a iB(r Ein )=0. (61e) 


Die Auflésung der Gleichungen (61a), (61b), (61¢) ergibt 
> Y ii ee yee) 
Hs =a(" 4-2 In| 1 ); Bie 0G ae le: 1?) + a@In@—"] , (62a,b, ¢) 


Mit (62a,b) erhalt man schlieBlich aus (38) 


2 
ets 
a2 


Y r2—a 
w" = r2 2 ! 2 ale See ed 
8B [ e te ere te —a 


rcos@. (63) 
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Mit a = 0 und B = S erhalt man fiir homogenen Spannungszustand aus (62a,b,c) mit anderer 


Konstantenbezeichnung sofort wieder die Gleichungen I (53a,b, c), wahrend (63) wieder in I (13) 
iibergeht. 


SchlieBlich erhalt man aus (54) und (63) die Gesamtauslenkung der Kreisringmembran mit frei 
beweglicher Innenkreisplatte durch Superposition zu 
ae 


ww! + wl = 3B (2 — r®) + a4 In 4 — 


9 


a 


= (2h + rcosq). (64) 


Fir homogenen Spannungszustand (B =S, a=0) geht (64) sofort iiber in die fiir denselben 
geltende Gleichung in I (14). 


Ein FANE eal fiir diesen Fall vollkommen freier Innenkreisplatte ist in Zit. 7 berechnet. 


6. Die Kreisringmembran mit elastisch gelagerter Innenkreisplatte. Auch im Falle eines in- 
homogenen, achsensymmetrischen Ausgangsspannungszustandes soll, um méglichst allgemein zu 
peaben: elastische Lagerung mit der Konstanten c, gegen eine Plattenverschiebung in Richtung 
w’ und eine solche mit der Konstanten c, gegen Verdi tier um die horizontale Paar angenom- 
men werden. Erstere beeinflubt die achsensymmetrische PORE ORY: w’ unter dem Drucke p! =yh, 
letztere die antimetrische Auslenkung w’’ unter dem Drucke p’!=yr cos (man vgl. I, Abb. 9). 

Beim achsensymmetrischen Belastungsfall gelten fiir w' baw. dw'/dr wieder aie Gleichungen 
(20a) bzw. (19a). Bezeichnet man die Plattenverschiebung wieder mit s und bedenkt, daB das 
Gleichgewicht der Innenkreisplatte gegen Verschieben, sowie die Randwerte von w’ am Innen- 


und AuBenrand die Erfiillung der Bedingungen 
8,(Gp) 277 +yhrin—ce,s=0, wir =F) =\s; wi (r =i7,) == 0-<(65.a, b,c) 


erfordern, so erhalt man mit (19a) und (20a) fiir D,, D, und s die folgenden Gleichungen: 


aus (65a): —22D,—c,s=0, (66a) 
h 2 
aus (65b): ieee (a's at + 7] In Lie te (66b) 
lyh lyh D 2 
aus (65¢): es (3 Ma? typ) E 14 De (66¢) 
Fiihrt man nun die abkiirzende Bezeichnung 
Cs, r2—a? 
Rigas Oita +22B (67) 
ein, so ergibt die Auflésung von (66a,b,c) nach kurzer Rechnung 
J eer 
D, = is r soe | 
AN \ 2 r?— a? 
2 (68a, b) 
Doe i Lay [9S ina pea? a ea ED oa BOREL cae) 
Dice 4B ae as ORIN al ie r?— a? a? i 
a ry 
Ss ON & —ri+aIn = a ; (68c) 


Fiir c, =00 folgt aus (68c) wegen N, =oo sofort s=0, wie es dem Fall mit fester Innenkreisplatte 
entspricht. Fir c, =0 hingegen eG sich sofort (55), wie es dem Fall mit freier Innenkreisplatte 
zukommt. 


Mit (68a,b,c) erhalt man schlieBlich aus (20a) [man vgl. hierzu die spatere Gleichung (76)]: 


= Pm eae? Be 8) CO) aa B| In” = (69) 
re a 


= 1 
Ga ee iS « ine A B + [e, (r? — r3) —- a2 4. BY zal ; (69a) 
isha a* 
Die Ausdriicke (69) mit (69a) erfiillen sowohl die Differentialgleichung (18), als auch die Rand- 
bedingungen (65 a,b,c). Insbesondere folgt aus (69) mit N, in (67) fiir c, = co, wobei also bloB die 


Glieder mit c, beibehalten werden, 


, _ yh 1 l ra— a @ 72) __ (7? — 7?) In 
es Fee Pe er as 


=| hn 
2 rd 
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Nach kurzer Umformung erhalt man hieraus aber wieder Gl. (42), wie es dem durch ec, = oo ge- 
kennzeichneten Sonderfall mit fester Innenkreisplatte entspricht. Setzt man hingegen in (69) 
mit (67) c, = 0, so folgt sofort wieder (54), wie es dem durch c, = 0 charakterisierten Sonderfall 
mit freier Innenkreisplatte zukommt. La®t man andererseits a gegen Null gehen und ersetzt B 
durch S, so folgt zunachst aus (67) 

Ne = ¢,In* + is se 


U 


also N in I (57), wahrend aus (69) folgt 


‘oes re c, In = + 42 s) (r2 — r2) —¢, (r?— r?)-2In a 


oder 
i yh 


> ASIN 


, 


ea") 228 Ce Tos ina C In - iE ne 


Dies ist aber Gleichung I (60), wie es dem durch a = 0, B = S gekennzeichneten Fall mit homo- 
genem Spannungszustand zukommt. 

Um fiir den Fall mit elastisch gelagerter Innenkreisplatte auch das Volumen des Wélbhiigels 
zu erhalten, hat man (69) wieder der in (43) angegebenen Integration zu unterwerfen, wobei man 
(44a) ausniitzen kann, und zu dem Ergebnis das Zylindervolumen unter der Innenkreisplatte, 
namlich r? zs hinzuzufiigen, wobei s durch (68c) bestimmt ist. Nach langeren Rechnungen findet 
man 


"2 2 1 7 
Vasa gi ty 2 ane e aa 
2 Ese 


Fiir den Sonderfall mit fester Innenkreisplatte erhalt man aus (70) und (67) mit c,=oco unter 
Beibehaltung nur der Glieder mit c, 


V=— hn Vas 3) (" matin a a In eae (i a) : 


1, r?—da A 2, r? — q? 4 y 
Wig? aye eo ; 
ee 2 


Dieses Ergebnis stimmt aber voéllig mit (47) iberein, wie es dem nunmehr gekennzeichneten Falle 
mit fester Innenkreisplatte entspricht. Setzt man andererseits c,—0, so erhaélt man aus (67) 
zunachst N, = 2a B und aus (70) nach Kiirzung durch Ba 

ha r2— a? 
=fR|(¢—7) + 2 a* (r? — r?) + 2 a erase | 


a2 


V 


Dies stimmt aber véllig mit Gleichung (57) iiberein, die fiir frei bewegliche Innenkreisplatte giiltig 
ist. Liegt schlieBlich ein homogener Spannungszustand vor, so erhalt man mit B = S und a = 0 


aus (67) und (70) 
N,=¢,n242nS=N _ in I(57), 


Sy hal Tol A te : f Nas 
V = toms rj) In 5 (Ta rife, + tr rin S|. 


Dies stimmt nach Vorziehung des Faktors 1/2 aber véllig tiberein mit der Gleichung I (61), die fiir 
homogenen Spannungszustand und elastisch gelagerter Innenkreisplatte giiltig ist. Teilt man 
noch V in (70) durch 727, so erhalt man 


= yh Vr 2 9 2 r2? — @? 1 Gy? 
Wmitel ~ 7B N, | 4 = 7 (2st 20°) Iie ar ; a z c, 
Feng) 2 2 a* r2 — q? 
| “a 2 (TS ary yond") 2 In eae Bah. (71) 


Beim antimetrischen Fall mit der Auslenkung w” infolge der Belastung durch den Druck 


p’ =yrcos¢ hat man die Gleichung (38) mit den Integrationskonstanten E, und E, wieder den 
Bedingungen (59b,c) zu unterwerfen, wahrend statt (59a) nunmehr die allgemeinere Bedingung 


Mee —— Onn rea Mom Ss’ y Aye fo) —0 (72) 
v v 
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infolge der elastischen Lagerung gegen Verdrehen in Geltun 
g ist. Hierbei kénnen die Gleichunge 
(60a), (60b) und (60c ttelb ee 
Ce a ae elbar beniitzt werden. Fiihrt man diese in (72) ein, so erhalt man 
jaa Ais zal 


a? 


1 y r? a2 E a2 a a2 a) 
Pat (1 P ) E, 5) a in ( 7) r = | AB 3 ie L a? In eee 
ae a2 Cd a. 
(1 7° = in6=0. (72a) 


Die Gleichungen (61b), (61c) kénnen hingegen unmittelbar iibernommen werden. Im Verein mit 


(72a) ergeben sie fiir E,, E, und 6 die Eatin 
In (1 — iy 


E, = ane r( In z In (1 | In 


+r vee In (1— “al cot |r + a* In 


1 1a 


2 ax at, (73 a) 


eis a ; Age 
te 4B M, (" ie | a2 In 2 =) Gre (73 b) 
y A r2— a 
o=gpu,(2—i + oni )2ea B (73) 
" mit 

ae 
ME AO ay Bee leg So nach 4 
1 hea d a 7? rr a?) Ca (74) 

a) 


Fiihrt man (73a,b) in (38) ein, so erhalt man schlieBlich fiir die Auslenkung 


n __yreosp{|, a’ r? — a a r2— a? a r2— @ 
MS ep M, il jn (1 re)! Sane! In (! as In ra L In ( re) i Tan 


2 / 2 / 
+ (r2 — r?) In (1 3 (r2 len ul — =) — (r? — r?) In (1 — a] C4 


— 
+(n —r+a?In i ae (75) 
Bei der Berechnung von w"’ in (75) wurde von der pia 
if ibe 
In | 1 = In 1 =n a (76) 


Gebrauch gemacht, die wegen (12a), namlich wegen a? < r? nak ist. Gleichung (75) befriedigt 
alle Randbedingungen (59b,c) und (72) und natiirlich auch die Differentialgleichung (21) bzw. 
nach Weglassung des Faktors cos die Differentialgleichungen (23) oder (23a). 

LaBt man c, gegen Unendlich gehen, so bleibt von M, in (74) und ebenso in (75) nur der 
Term mit ¢,, and man erhalt sofort wieder die fiir feste Innenkreisplatte giiltige Gleichung (50). 
Setzt man andererseits in (75) cy = 0, so bleibt von M, in (74) und ebenso von (75) nur der 
Term 2 a?z B, so daB man nach Kiirzung durch diesen Faktor sofort wieder Gleichung (63) er- 
halt, die fiir den nunmehr erhaltenen Fall mit freier Innenkreisplatte giiltig ist. Setzt man schlieB- 
lich in (74) und (75) a = 0 und B = S, so erhalt man sofort wieder I (68). Teilt man namlich 
Zahler und Nenner in (75) durch a”, so erhalt man nach Weglassung der ersten drei mit lim a—>0 
wegen In 1=0 verschwindenden Terme des Faktors von c; und ebenso des mit a? behafteten 
Termes im Faktor von 22S, wenn man iiberall wieder B durch S ersetzt und (74) beachtet, 


F yr cos p 


S tat ae) mY : 
a8|iln a az) 4 T 2a | 
| 


2 2 2 2 2 2 2 »2 2 
Las a corals ES Se ee V5 as ; 2 2). 
Ag In ( =) ao ae In {1 ; = In 72 cy t+ (2 —7*)-2a 8}. 


Beachtet man nun die Grenziibergange [man vel. (39)] 


* 1 2 1 mar aE: l . I 1 
lim Ut er steer lim a In (1 7) = a ie % LE eae as 


2 
a—0O uf a—0 
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und schlieBlich 
; ti r? (r2 _ a”) Ta — r? 
lim - In 4-4 7 = 4S 
a TA 1a") aire 


° 
a—0 


so folgt nach Erweiterung mit r2ré sofort die Gleichung I (68), da 


2 2 
| (r* | te) | 


a 2 ae 
Ya ali 


2 2 2 2 2 2 2 2 
Pia) ia if Uf r Li LP Ti 72 
— = r2 


ist. SchlieBlich erhalt man durch Addition der Teilauslenkungen w’ in (69) mit N, in (67) und w’’ 
in (75) mit M, in (74) die Gesamtauslenkung w = w' + w' zu 


y_ fh toe , Ne (se Zoya\ et ate Be ee 
Ww = vain, (« Inge t 4B) Ce el 1, (re ri) da a By nae 


: 9 9 9 9 e 2 2 Pe pet 
r cos P 5 a” Te a’ Teo ; a ee a 
+S [fein ae) inte n(n + In(1 r ae 


a 
2 


¢ 2 
+ (r? — r?) In ( ") (r?2 — r?) In (1 -- 7) — (r? — 7?) In (1 — =) C4 


r2—a?\ 
an (" —r+ta?ln fe =a) -2a°x a} 3 (77) 


Erteilt man w in (77) den konstanten Wert w,, so stellt (77) fiir bestimmte Werte des Halbmesser- 
verhaltnisses 4 = r,/r; und der Riickfithrungskonstanten c, und cy, sowie des Inhomogenitats- 
faktors a des Spannungszustandes die Héhenlinien der Wélbflache der Kreisringmembran in Polar- 
koordinaten dar, deren punktweise Berechnung nach Wahl der Mittelpunktsentfernung r und dem 
daraus bestimmten Wert von cos @ leicht erfolgen kann. Man gewinnt so durch (77) eine fiinfpara- 
metrige Kurvenschar, von welcher natiirlich blof8 die Zweige innerhalb der Flache der Kreisring- 
membran physikalische Bedeutung besitzen. 

Ebensolche Kurvenscharen gewinnt man aus (77) auch fiir konstante radiale und konstante 
azimutale Neigungen des Wélbhiigels durch Bildung der Ableitungen (dw/dr), und (dw/r 0¢),, 
denen man nachtraglich konstante Werte erteilt, worauf der Zeiger c in denselben hinweist. Bei 
den Kurven w, = konstant sind w,, “4 =1,/r; v = h/r,, ¢, und cy die Parameter, wahrend bei den 


Kurven konstanter Neigung statt w, jeweils eine der oben angegebenen konstanten Neigungen als 
Parameter in Erscheinung tritt. 


7. Zahlenbeispiele. Wie im homogenen Falle empfiehlt es sich auch hier, durchweg dimensions- 
lose GréBen einzufiihren. Um die Auswélbungen, die zu verschiedenen inhomogenen Spannungs- 
zustanden gehéren, miteinander in Vergleich setzen zu kénnen, wurde bei allen diesen Spannungs- 
zustanden die gleiche Sparnung S, = S am Aufenrande der Membran vorausgesetzt. Im Falle 
eines homogenen Spannungszustandes wird, wie dies ja schon oft benutzt wurde, auch die GréBe B 
in (la,b) gleich S, wahrend a = C = 0 ist. Neben den bereits in I eingefiihrten dimensionslosen 


GroBen 


3 


* (78a—g) 


omit w mk 
somit = — 

% 8 S 
ist hier noch eine GréBe A einzufiihren, die neben a in (4) die Inhomogenitat des Spannungszustan- 
des charakterisiert und die am besten durch das Verhaltnis der Radialspannungen am Innen- 


und am AuBenrande beschrieben werden kann, wobei das beigefiigte Intervall fiir 2 aus (9b) oder 
(7b) mit r =r, folgt 


r ) a 
) =. —— —_ — == 0 Tt. SS — —_— = pp) = — 
Or ant GB iat are ane a} 


4 Si Si 2 
Wy = ‘Se a FS) ° Ane, = 0 = A Ss pe ip — Agee G (79 a, b) 
Im homogenen Fall ist A = 1. Die Gleichheitszeichen in (79b) beziehen sich auf die in Ziff. 1 
erwahnten Grenzfalle, bei denen einmal $;=0 ist (A =0), das andere mal So; = 0 ist 
[bei A = 2 w?/(u2+-1)] [man vel. (5b)]. 
ir die oftmals gebrauchte, die Inhomogenitat des Spannungszustandes ebenfalls charak- 
terisierende Konstante a in (4) mége die dimensionslose GréBe 


5: ar ae —A il ; 1 
f= = qe (+ a), = pz e Dili (80a,b) 


mit § = 0 im homogenen Falle eingefiihrt werden, wobei das Intervall (80b) durch (79b) bedingt 
wird und die beiden Gleichheitszeichen sich wieder auf die erwahnten Grenzfalle beziehen. SchlieB- 
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lich findet man aus (3a) 


Sl 1 w—l w+ 
pee Ty Seay ( oh S*s | we i Ble 2) 
Im homogenen Fall ist B = S und also x = 1. Zwischen fb? in (80a) und x in (81a) besteht der 
Zusammenhang " 
1—f2=x, (82) 


der spaterhin ausgenutzt wird. Ferner ergeben sich die folgenden, spater benétigten Beziehungen 


[man vel. (80a) ] 


1—2/ 5j 
r2— @? ees ue | 
SCE aera ieagee. Reoey lial Be (83) se a Grenzfall S.-0 
a yy x S A=0=A min 
pa @ nat 0 (Be =A) (1 —A) Fg? = BP 5 
Seer ape Ne pa re ia Re Ee = = 10 
ri—a A w—l Se ie) tee a 
Cat ae UA) OU) 
r—aq we—l =i me Ce, 
Fir den Spannungsverlauf des Ausgangsspannungs- 
zustandes findet man nach (7a,b) mittels (78) bis (81) 


in dimensionsloser Darstellung 


Sth (jc Sa ee 
S x o |}? S y ee 
hieraus folgt ‘ Homogener Fall 
= A=1 

Pee Oe 2 2) eA 

are > 2 ae a konst., 
also unabhangig von 0. 

Der Verlauf dieser Spannungen als Funktionen von 9 ie 
ist in Abb. 2 fiir die angegebenen j-Werte dargestellt. ea = S Grenzfall Sp=0 
Mit zunehmendem / nimmt S, + S, ab und erreicht ¥ A=18-A ma 
seinen Kleinstwert fiir 2 =/,,,, = 1,6 S (fir uw = 2) und 0=95 o=10 
den GréBenwert 2,67 Ss fur, = An 0) Allgemein gilt Abb. 2. Der Spannungsverlauf von S_ und S_, als 
9 fs we S; + So & 9 : we Funktionen von 9 fiir die angegebenen Werte von /. 
(eae YS s) =) p=)" 


a) Feste Innenkreisplatte. Die Einfiihrung der dimensionslosen GréBen (78) bis (85) in 
die fiir diesen Fall giiltige Gleichung (51) mit w, statt w ergibt 


a 
—_ —- VW, 0 
Ge y(0) 


cos y = ——__.—_,, 
Proxiit p2(Q) ; 
; pee ENG? =" 
| vile) = 2 (1 — 9") + 2 Tae ai (86a, b, c) 
und! A 
2 f8Iny + 1—— ay 
(oe) = 0\- ie In op 2 2 Ino + 1— @?}. 
1 Da nach (8) und (83) : 
Dieleee ON HE \ 
Neal & Spe: 7) hn es 2) oe Ind 


ist, so erhalt man aus (51) wegen (84) und (85) 
2__ Pp: Py eh) 2 l 
Whe gal Se Nad ted PE oa nC al hl a gl | In 2¢ 
vi=14 Ig D te 7 In x —In * In (1 — B? «?) — In o n r(@ ie 


} ip In AIn | 


1 2 pe 
+= ein — FA) —(1— J) n S| = Oa 


2 1 gr x 
+ (1— 0?) nA + ( =a] In (aa Q 
wenn man beachtet, daB wegen (80a) und (81a) x: (1 — f? u*) = 1: A ist. 


2 9 


Q 
In w In ES B 
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Die Gleichungen (86 a—c) erfiillen, wie man sich zur Kontrolle iiberzeugen kann, die Bedingung 
x — 0 fiir den AuSenrand g@ = 1 und den Innenrand 9 = 1/u fiir Fak nach (79b) zulassige /. 
Dazu ist notwendig und hinreichend fiir reelle p, daB fiir die obigen Werte von 0 y,(0) =Y2(0) = 9 ist. 
1, 01. 2Dab. hier y,(1) = 0 ist, erkennt man daraus, da nach (82) 
Oita De meen peg 
x 4 


ist. Entsprechend folgt y,(1) = 0, da 


ee 
Cpe 
wird. 
ena ame 
1 P é i \ ee a ate pe fie A 
22-0 air" Hier wird y, Gs == 2(1 ia] | 2 7, In 7 <= (Dy 
da ' 
( ae é*) - == Te 
we Ww 


OP Us 1 We pret vw. 
Wo fa == (\) 5 da op = 5 ist fiir Oe ey 
Im weiteren wurde 4 = 2, y = | die Spiegelgerade somit als obere Membrantangente gewahlt. 
Nach Einsetzen dieser Werte kann man die Hilfsfunktionen y,(0) und y,(0) in (86) wie folgt dar- 


stellen. 


Die Funktionen y,(@) und y,(0) in (86b,c) fir wu = 2, » = 1 und die in der ersten Spalte an- 
gegebenen Werte von A sind: 


Nie 
ae) = 


ef(1 — o) — 0,5 Ing] 


TE 2(1— 9) + 0,3423 In {1,3167 (? — 0,2405)} 
A = 0,05 
volo) = 9 [(l — o2) — 1,2044 In 9 + 0,3617 In {1,3167 (g2 — 0,2405)} ] 
(0) = 2 (1—p2) + 2,1640Ine =—=Z in I (43) (87) 
all 5 1 


N in I (43¢) Man vegl. die folgen- 


den Ausfiithrungen 


3 


NBO hetieie a0 ites 


2 (1 — g%) + 1,637 In {0,8 (0? + 0,25)} 
0 [( — e2) + 2,2200 In g — 0,8600 In {0,8 (o? + 0,25)} ] 


Driickt man in p,(@) (86c) x und f? mittels (81a) und (80a) durch J aus, so folgt 


In ge” (uv? — 1) 
1—A 1 op TET (Le) 5 lA 
polo) = @ pees Ings + 1 4 Ind a elon In 90 + 1— 93}. 
LaBt man hierin A gegen 1 gehen, wie es dem homogenen Fall entspricht, so verschwindet in der 
runden Klammer der Term mit In 4% und ebenso der mit In go, wahrend der Bruch einen unbestimm- 
ten Wert annimmt; seine Berechnung ie 


See 1—¢? 
lim @ ue wpe Sassy 
G4 In A e*(w?— 1) 


Es verbleibt 
; 1 1 — o? 9 5 1 
vx(ohi=1 = 0| (1 3) GED eed |= el (1 we #): 
Mit uw = 2 folgt 


5 1 
VAC) ee eee ee also —N in 1(43¢). 


Fir y,(@),-, ist dieser Grenziibergang noch’ viel leichter durchfihrbar. 
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In Tab. 1 sind i i i iir die i 
eee ae Funktionen Yr und W. in (87) fiir die in der ersten Zeile angegebenen Werte 
ee 2 ES ° ie erte von A jeweils in den ersten beiden Zeilen berechnet, wiihrend die nachsten 
u ; j i i 
4 r die Me er zweiten Spalte jeweils angegebenen a- Werte die nach (86a) berechneten Werte 
te gy entha ten, die dann zur Konstruktion der Abb. 3 bis 5 dienten. Darunter, durch einen 
‘ : A ao Seth 
a fens befinden sich drei Zeilen, welche die fiir cos y =—1,0, + 1 nach (86) ermittelten 
4 erte ent alten, die zur Darstellung der Meridiankurven in den Schnittebenen durch die &- 
zw. ¥-Achse in den Abb. 3 und 4 benétigt werden. 
- Der i ae 3 dargestellte Fall A =0, bei dem also am Innenrand S,= 0 ist, stellt einen Extrem- 
all auc i i i 
insofern dar, als er zwar den vorgeschriebenen Randbedingungen (40a,b) geniigt, wie 


oben gezeigt wurde, andererseits aber nach (86a) und (87) fiir A =0 zu jedem Winkel ¢p fiir 0 =0,5 
einen endlichen Wert ergibt. Man erhalt 


a =x [yy,(0,5) + y2(0,5) cos p] = 0,75 [1,5 + 0,5483 cosy] . 


Schnitt A-A 


0 
| 
| B 
| 
| 
05 
10| i = | ’ 
ae | : Ay 
A- | 
y | | Soni Bf" 


Abb. 3. Die Hohenlinien der Wolbflache der Kreisringmembran mit fester Innenkreisplatte fiir 2 — 0 und fiir die durch die 
angegebenen o-Werte bestimmten Héhenlagen der Schnittebenen, sowie die Meridianschnitte durch die €- und 7-Achse. 


Z. B. erhalt man so nach Tab. | fiir 7=0 und cosg=1 «=1,536. Dieser Fall mit sehr kleinem 
A 1aBt sich verwirklichen (Abb. 5), wenn sich z. B. die Membran an einen Zylinder vom Halb- 
messer r, anlegen kann. Sie schmiegt sich demselben und ebenso der Meridiankurve 1 =0 um so 
inniger an, je kleiner J ist. Das ganze Verhalten ahnelt dem in I, Abb. 4 dargestellten, bei dem 
allerdings nicht die Inhomogenitat des Spannungszustandes, sondern r; verdindert wurde. Wiirde 
in Abb. 5 auch noch r; variiert werden, so wiirde bei jeweils gleichem / die Anschmiegung an den 
Innenzylinder und die Meridiankurve um so inniger werden, je kleiner r; ist. 

In Abb. 5, in der auch der vertikale Meridianschnitt der Vollkreismembran fiir den homogenen 
Fall 4=1 (strichliert gezeichnet) aus I itbernommen worden ist, erkennt man den grofen Einflub 
der Inhomogenitat des Spannungszustandes sehr deutlich. Bei kleinem /, also kleiner Innenrand- 
spannung, wélbt sich die Membran wesentlich héher auf (in Abb. 4 und 5 fir A=0,05 wélbt sie 
sich bis iiber x =0,7 auf), als im Falle gréBerer 2 (in Abb. 4 und 5 fiir 2 =/,,,, = 1,6 nur bis nicht 
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Schnitt A-A 


O17 


05 


-08 


05 


Abb. 4. Die Héhenlinien der Wélbflachen der 
sae e Kreisringmembran mit fester Innenkreisplatte 
fiir 2 = 0,05 und 2 = LN Se = 1,6 (je zur Halfte 
dargestellt) und fiir die durch die angegebenen 
a-Werte bestimmten Héhenlagen der Schnitt- 
ebenen, sowie die Meridianschnitte durch die 

& und »-Achse. 


Schnitte A-A 
durch die n-Achse 


= \Valikresmembran 
as 


Abb. 5. Vergleich der Meridianschnitte der Kreieringmembran mit fester 
Innenkreisplatte durch die hag chse fiir die angegebenen Werte von /. 


12* 
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ganz 0 =0,4), sofern die AuB®enrandspannung in beiden Fallen dieselbe res Zugleich erkennt man, 
daB der Wolbhiigel fiir kleine Werte 2 am Innenkreise wesentlich steiler aufsteigt, als am AuBen- 
kreise. Auch dieses Ergebnis entspricht der Anschauung. ‘ 

SchlieBlich erhalt man fiir das Volumen V nach (47), wenn man die dimensionslosen GréBen 
B® (80a), x (81a) und (83) beachtet und wieder S, = S setzt 


il 1 
(pee eee 
V l vv re | ie me ‘ jie “2 see — alt a." pe (88) 
po gh Pe 2 Ble A. othe ek ee 
Tat 2 A 
oder 
A 
= il = 2] 
ae ees PO 
A Sp? (w—A) 2 (u2—A) ele ("") 
A 
Il 3 
nat A en ie a (eet fo a (88 a) 
4 S — pw? (u?—A)| 2 (uw? —A) ele (7) 
A 
Fir 2 = 0 ergibt sich insbesondere 
Ere a he Ng 88b 
k= esa a je ; ( ) 


Aus (88) ist ersichtlich, wie V von yw und J abhangt. Insbesondere erkennt man, da V mit ab- 
nehmendem / zunimmt, wie es auch die folgende Zusammenstellung erweist: 
A 1,6 l 0,05 0 


mi) 
yi 0,0416 0,0495 0,0955 0,166 


Fiir 7 = 1 folgt aus (88) fiir homogenen Spannungszustand 


Vore b2— 1 


[pe hae eee ee es 
= hSln eos [(u? + 1) Inu — (wu 1)] 
wie in I. Se 370. 


b) Vollkommen frei bewegliche Innenkreisplatte. Mit den dimensionslosen GréBen 
(78) bis (82) und (85) erhalt man mit w, statt w aus (64) 


vie —— Dy 
60s @ _ ae v wo) 


0 yo) ? WG) lO a re 


Fir homogenen Spannungszustand x =A =» = 1 erhalt man aus (89) y(0) = 1 — 0? und damit 
sofort wieder I (24). Mit den Werten «4 = 2,» = 1 und den bereits in (87) a 


erhalt man insbesondere aus (89b) 


(89 a,b) 


ngegebenen A-Werten 


A=Anin = 9: yo) = 1— oe? — 0,25 In [1,333 (o2 — 0,25] , 
Yen O0or v(o) = 1 — o? — 0,2405 In [1,3167 (0? — 0,2405)] , 
: 5 (90) 
Ava ls y(o) =1—¢?, 
A =Anax = 1,6: w(o) = 1— oe? + 0,25 In [0,8 (9? + 0,25)] . 


In Tab. 2, deren Aufbau ganz dem der Tab. 1 gleicht, sind fiir die oben angegebenen Werte 0 
und die in (90) hervorgehobenen Werte von J die Funktion y und in den nachsten Zeilen fiir die in 
der zweiten Spalte jeweils angegebenen ~-Werte die nach (89a) ermittelten Werte von cos p be- 
rechnet worden, die dann zur Konstruktion der Abb. 6 benétigt wurden!. Wieder durch einen 
Strich getrennt, enthalten jeweils die nachsten drei Zeilen die nach (89a) ermittelten «-Werte fiir 
cosy = — 1,0, + 1, die der Darstellung der Schnitte 4A— A, B—B und C—C in Abb. 6 und 7 
dienten. Aus Abb. 6 erkennt man ein analoges Verhalten wie in Abb. 4: Bei dem kleinen Werte 
A= 0,05, also kleiner Radialspannung am Innenrand, wélbt sich die Membran samt der frei 


' Dies wurde fiir 2 = 0 unterlassen, da sich der Hohenlinienplan dann nur sehr wenig von dem fiir 4 = 0,05 
unterschieden hatte, wie man auch aus den Meridianschnitten durch die 7-Achse in Abb. 7 erkennt. 
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beweglichen Innenkreisplatte viel héher auf, als im dort ebenfalls dargestellten Fall A= /,,,, = 1,6. 
Auch der Anstieg am Innenrand ist dann viel gréBer als bei gréBeren J-Werten, ein Ergebnis, das 
wieder der Anschauung entspricht. Noch deutlicher wird dies aus den in Abb. 7 in Vergleich ge- 
setzten Meridianschnitten durch die 7-Achse, die zu den verschiedenen j-Werten gehéren, wobei 
der strichliert gezeichnete Schnitt wieder aus I fiir homogenen Spannungszustand entnommen. 
wurde. Der singulare Grenzfall 2 =0 ist hier durch unendliche, praktisch sehr groBe Steilheit der 
Meridiankurve am Innenrand und unendliches (praktisch sehr groBes) Erheben der Innenkreis- 
platte gekennzeichnet, er ist also bei frei beweglicher Innenkreisplatte auszuschlieBen. Fir uw =2 
und y—1, wird also etwa 2 0,05 sein miissen, da sonst die Auslenkungen nicht mehr als klein, 
wie es die Theorie verlangt, angesehen werden kénnen. 


Fir das Volumen V des Wolbhiigels samt Innenkreisplatte erhalt man aus (57) mit (80a), 


(81a) und (83) 
ITS eae pA—1 1—A\ 2—1 1—A\2 Ww 
ire ei w2—A 2 ut | feed [tie (aoa) in ())- (91) 


Hierin ist die Abhangigkeit des Volumens V von den Parametern y und J ersichtlich. Fiir homo- 
genen Spannungszustand folgt aus (91) mit 7 =1 wieder 


Ty Te 


Piste 
} 8S E ie (91a) 
woraus man fiir 4=2 und »y=1 sofort wieder 
Neyae 92772 
i 1238S 


wie in I S$. 374 erhalt. Aus der folgenden Zusammenstellung erkennt man, daB V mit abnehmen- 
dem J zunimmt: 


A ail siehsb 1 0,05 0 


De 
V:5-2) 0,333 _ 0,368 0,532 oo 


Insbesondere folgt im singularen Grenzfall 4=,,,,=0 aus (91) in Ubereinstimmung mit Abb. 7 
Vi== OO; 


SchlieBlich findet man aus (55) die Auslenkung s und aus (62) den Winkel 6, fiir welche GroéBen 


nach Umschreibung in dimensionsloser Form folgt 


the 1 y re ee —Il i we 
ee ue . eA He (5 ) : (92) 
Re wlavite (ie) 1—A/ Tin 
y CaS) sae we | rs) In G5 ) : (93) 


Der Wert oo fiir s und <6, wenn A gegen Null geht, kennzeichnet abermals die Singularitat des 
Grenzfalles A=0 bei freier Innenkreisplatte. Man erhalt aus (92) und (93) 1: 


2 1,6 | 1 | 0,05 0 
yr 
sie = as 1,302 1,500 2740... |= 0 
r2 a 
6: re => 0,651 0,750 1,370 co 


Auch der Fall der elastisch gelagerten Innenkreisplatte laBt sich analog behandeln, wenn man 


neben den in (79) bis (85) eingefiihrten dimensionslosen GréBen und Funktionen noch die beiden 
dimensionslosen GréBen 


% Cs x Cd 


on 8” 1 taSr8 (94.a,b) 


* Die Bezeichnung o, und 6, in Abb. 6 und 7 soll auf die Abhangigkeit di 6 ie 
gewahlten MaBstaben fiir apBiniten, andi eHe cere uf die angigkeit dieser GroBen der Zeichnung von den 


. 
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en 9g Abb. 6. Die Héhenlinien der Wélbflichen der Kreis- 
-10 ——s a ringmembran mit frei beweglicher Innenkreisplatte 


- <= fir 2 = 0,05 und 2 = Dek = 1,6 (je zur Halfte dar- 


é gestellt) und fiir die durch die angegebenen «-Werte 
Schnitte A-A bestimmten Héhenlagen der Schnittebenen, sowie 
durch die 7 -Achse die Meridianschnitte durch die &- und 7)-Achse. 


Abb. 7. Vergleich der Meridianschnitte der Kreisringmembran mit frei 
beweglicher Innenkreisplatte durch die 7-Achse fiir die angegebenen 
Werte von A. 
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hinzunimmt, die die elastische Lagerung gegen Verschieben bzw. Verdrehen der Innenkreisplatte 
kennzeichnen. Aus (77) mit N, in (67) und M, in (74) 1aBt sich dann leicht eine Gleichung analog 
(86a) oder (89a) fiir cos p angeben, die dann Héhenlinienplane wie in Abb. 4 oder Abb. 6 zu kon- 
struieren gestattet. Da aber das zu erwartende Ergebnis im wesentlichen auch schon aus den in 
den vorhergehenden Zahlenbeispielen behandelten Grenzfallen dieses allgemeinen Falles hervorgeht, 
wurde von einer Durchrechnung desselben Abstand genommen, zumal in I, Zahlenbeispiel 4 die 
Auswélbung der Kreisringmembran mit elastisch gegen Verschieben und Verdrehen gelagerter 
Innenkreisplatte bei homogenem Spannungszustand behandelt worden ist und in einer spateren 
Arbeit iiber Laminarstr6mungen ein analoger Fall berechnet wird. Eine Singularitat wie im vor- 
hergehenden Falle wird wegen der elastischen Lagerung der Innenkreisplatte nicht in Erscheinung 


treten!. 


(Eingegangen am 10, August 1957.) 


Anschrift des Verfassers: Prof. Dr.-Ing. Karl Karas, Darmstadt, Technische Hochschule, 
Institut fiir Angewandte Mechanik und Technische Schwingungslehre 


_ + Eine hydromechanische Deutung der Auswolbung einer inhomogen gespannten Kreisringmembran wiirde 
nicht nur eine 6rtlich veranderliche dynamische Zahigkeit der strémenden Flussigkeit, sondern iiberdies noch 
deren Abhangigkeit von der Orientierung der Flachen voraussetzen, in denen sie wirksam ware. Hingegen ist 
eine solche Deutung bei homogen gespannten Membranen stets moglich, Fiir rechteckige und elliptische 
Konturen — auch in allgemeiner Lage gegen den Flissigkeitsspiegel — wurde dies vom Verf. in bereits abge- 
schlossenen Arbeiten gezeigt. Wahrend sich die Ergebnisse fiir elliptische Konturen in geschlossener oem 
angeben heBen, boten sie sich fiir rechteckige Konturen in Form auBGerordentlich rasch konver enter Reih 
dar, die eine zahlenmaBige Auswertung leicht zulieBen. ea 
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Einflu8 eines Kreisloches auf die Durchbiegung einer diinnen Halbebene 


Von O. Tamate 


1. Einleitung. Vor kurzem hat S. Woinowsky-Krieger! cine Reihe von Randwertaufgaben der 
Plattenbiegung mit Hilfe von Bipolarkoordinaten durchgefiihrt; es handelt sich dabei naturgema} 
um Bereiche, die von zwei Kreisbégen oder einem Kreisbogen und einer Geraden begrenzt sind. 
Dabei wird erwahnt, daB das zur gleichen Klasse gehérige Gebiet, z. B. der exzentrische Ring, in 
grundsatzlich ahnlicher Weise behandelt werden kann. Die durch einen Kreis und eine Gerade 
begrenzte Platte und die durch zwei kreisférmige Locher geschwachte Platte gehéren auch in 
dieselbe Kategorie. 

Wir betrachten nun die Aufgabe, die Spannungserhéhung an einem Kreisloch fiir die diinne 
biegebelastete Halbebene zu bestimmen, und behandeln dieses Problem auf Grund der Kirchhoff- 
schen Biegungstheorie elastischer Platten? und der komplexen Darstellungsmethode von Muskhelish- 
vili3, Zur Bestimmung der Koeffizienten der den Spannungszustand der Platte charakterisierenden 
Funktion wollen wir die Stérungsrechnung verwenden. 


2. Problemstellung. In einer isotropen elastischen Halbebene (x, y) befindet sich ein Kreisloch 
vom Halbmesser r, dessen Mitte den Abstand d vom geradlinigen Rand y = 0 der Halbebene hat 
(Abb. 1). Durch die komplexe Funktion 


y 


¢ 
6 
wird der Bereich von der (x, y)-Ebene auf einen konzentrischen 
Kreisring 1 < @<, der (0,@)-Ebene abgebildet. Der Wert 
o = 1 liefert die Gerade y = 0, also die x-Achse. Soll 0 = 0 


dem Rand des Kreisloches entsprechen, so mu gelten: 


e=x+iy—ao(C)=ia = mith sa €F 8 a(21) 


e 1 2 
Noth etl ier aig Sees (2.2) 
yy Al Ci 
Nun nehmen wir an, daB die Platte der Beschrankung der Abb. 1. Die gelochte Halbebene. 


Kirchhoffschen Biegungstheorie elastischer Platten folge. Dann 
erhalten wir die Grundgleichung fiir die Durchbiegung w(x, y) einer auf ihrer Oberflache kraft- 
freien, nur durch Randmomente und Randkrafte gebogenen Platte 


o Ch 


Bekanntlich kann man die biharmonische Funktion w(x, y) durch zwei analytische Funktionen @p,(z) 


und y,(z) der komplexen Variablen 2 = x + iy darstellen: 


w(x, y) = Re [z (2) + x1(2)] - (2.3) 
Hierbei ist Re [..... | der Realteil des Ausdruckes, der in der eckigen Klammer steht; ferner ist 
z=x—iy. Die Spannungen werden durch die in Richtung der Dickenerstreckung zusammen- 
gefaBten Krafte und Momente dargestellt. Fiir diese GréBen gelten die folgenden Beziehungen: 


M, + M, =—2D (1 +») [pil) + pi)]; si 
=e / : z 
M,—M, + 2iH,, =—2D(1—») [zg{(2) + yi@)] mit y,(2) = Es = (2.4) 
N,—iN, =—4 Do;{(2), 
wobei M,, M, die Biegemomente, H,, das Torsionsmoment und N, und N, die Querkrafte sind. 
Jede dieser CroBen wird auf die Langeneinheit bezogen und die x- und y-Suffixe bedeuten die 


Wirkung auf dem Querschnitt senkrecht zur x- bzw. ‘zur y-Achse. Unter D und y ist hierbei die 
Plattensteifigkeit bzw. die Querkontraktionszahl zu verstehen. 


1 §. Woinowsky-Krieger, J. Appl. Mech. 22 (1955) S. 129; Ing.-Arch. 24 (1956) S. 47. 

2 §. Timoshenko, Theory of Plates and Shells, Kap. 4. New York 1940. Re ‘ 

8 N. I. Muskhelishvili, Some Basic Problems of the Mathematical Theory of Elasticity, Ubersetzung vom 
Russischen ins Englische von J. R. M. Radok, Groningen 1953. 
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Beim Ubergang zu den orthogonalen krummlinigen Koordinaten 9 = konst., = konst. mit 
Hilfe der Funktion z = @(C) spect die Formeln (2. 4) fiir die Komponenten an Momente bzw. 
der Krafte in bezug auf die ee 0,0 


M, +M,=—2D(1 +») [0 + ©], 


M,—M, +2i Hy, =—2D(1—») - 2 [oo(C) (6) + o'() POT. (2.5) 
p sae Me C 1 / 
No Ng = 4D eerie ee 


wobei zur Abkiirzung 
Az) =Pilo)] =~), — wl) = plo)] = v(C), 
aaa wane a (2.6) 
a0) =2@ , Me) =e 
gesetzt ist. 
Bekanntlich lauten die Randbedingungen fiir freien Plattenrand 


if id Sa. ffx . 
— do) +2910) +¥i@ Hic + G. (2.7) 
Dabei ist C eine reelle, C, eine komplexe Konstante mit 


7 
Servarae, (2.8) 
Wegen der Willkiirlichkeit in der Wahl der Funktionen ¢,(z) und y,(z) kann man auf einem der 
Rander die Werte der Konstanten C und C, festlegen. 

Mit diesen Bezichungen und Bezeichnungen bestimmen wir den Spannungszustand in der Um- 
gebung des Kreisloches unter der Bedingung, daf einerseits auf den Lochrand 9 = @) und den 
geradlinigen Rand 9 = | keine auBeren Krafte wirken, andererseits die Spannungen in groBer Ent- 
fernung vom Loch den folgenden Werten entsprechen: 


M() = M,, Mi) = Hl) = Ni) — Ni») —0. (2.9) 


Dann lautet die Funktion der Durchbiegung der Mittelflache der durch kein Loch geschwachten 
Platte 


My 
2DU arya x? — yy?) , (2.10) 


An Hand der Formel (2.7) kann man die Randbedingungen in der Form schreiben 


Uy = 


= 1 — 290) + 2 GO) + He) = 0, 2.11) 
@ = 00: — = Pl009) + ey F' (000) + ¥ (00.09) = G. (2.12) 


Elierin ist: o:==e)? 
Die Funktionen ¢(C) und y(¢), die den Spannungszustand der geschwachten Platte charak- 
terisieren, lassen sich darstellen in der Form 


pC) = PC) + PME), 
pd) = pC) + pO) , 
so daB die Funktionen p(¢) und y((¢) dem Grundspannungszustand und gC) und p(s) dem 


erganzenden Spannungszustand entsprechen. Infolge (2.10) haben die Funktionen )(£) und 
pO(2) die Form 


3 . M, l+e 'e ap 1 re . Mlt+el4+i1 
(0) — —_ ey 0 S 
pC) La— ae pl) = —ia—® — ares : (2.14) 


(2.13) 


Nun nehmen wir die Zusatzfunktionen in folgender Form an: 


pO) ma fi 0 hi ay: re & tee (eG Ae bc 2 , 
ri (2.15) 
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wobei infolge der Symmetrie des Problems a,, b,, ¢, und d,, reelle Konstante sind. Es ist leicht zu 
beweisen, daB die zu ~) und wp) gehdrigen Spannungskomponenten nur dann im Unendlichen 
ioe) 


verschwinden, wenn die Reihe Dr (a, — 6,,) konvergiert. Also ist es klar, daB die Funktionen p() 


n=] 
und y™ alsdann auch den eindeutigen und im endlichen Bereich der Platte endlichbleibenden 
Verschiebungszustand angeben. 


3. Koeffizientenbestimmung. Die Koeffizienten a,, b,, c, und d,, durch welche pe?) und yw) 
_ bestimmt sind, ergeben sich aus den Randbedingungen fiir ¢@ = 1 und 9 = ¢) und aus der Kon- 
vergenz der Reihe. 

Fithren wir die Funktionen p), yp (2.14) und p®, p™ (2.15) in die Bedingung (2.11) ein und 
setzen die Koeffizienten verschiedener Potenzen von o gleich Null, so erhalten wir 


bn = an 
ee ares Che ae (3.1) 
In gleicher Weise ergibt die Bedingung (2.12) mit (3.1) 
1 
ee ae (0) (1 0) 
l (3.2) 
(1—9;) 4 ++ (1 — 0%) =x @5, 
1 2n 1 ~2(n— eae 
ne (1 ae Co” ) Ce (n—l)e (1 — Qo ( )) a4 te (1 7 05”) os ) (6, ae b,-1) 
= — (105%) 05%", (n=2) (3.3) 
1 1 eerie 
( aw 0,”) (a,, —, a,,_1) i ne ( r 0,°")b,, (n= l)e (1 a \ 05 ) oF by-4 = 0 ? 
mit 
Os ee (3.4) 


Wegen der Rekursionsform dieser Gleichungen finden wir leicht 


1 2 2n = Fe — 2(i— 
R=.) es (19,7) 0, ° 9b, = (1-2, )—(1 Go)» Xb: 05” ier, 


ne 


(n =2) (3.5) 


ile 5 n—1 1 ce % 
ee 0, ) (L957) 3) (1 = 0,7) 0; 208," - 
ne a We 
Demnach liefert die Konvergenzbedingung der Koeffizienten 
co a” ioe) iT oe 
x =—1—(1—9,7? 7 8, 5-9 = (—1) Dae (eo) be (3.6) 
i=!) i= 
Dieses Gleichungssystem (3.2), (3.5), (3.6) reicht zur eindeutigen Bestimmung der Koeffizienten 
aus. Eliminiert man zunachst die Koeffizienten a, und x aus den Gleichungen, so kann das folgende 
Gleichungssystem fiir Ermittlung der Koeffizienten b, gewonnen werden: 


1 Le rae i! a 
fe tees amor (De st] Fg, —— (tet 
n—| 


1 2n =% —2% 1 --27 
SL Oye) (1 = 05° ) yaw )@ Tes (1-— 05? ) b; 


= 


1 I —2n\ ,—-2n —2\2 ,-2n nt 
[ra — oP +55 — PO — 00 0 cilia, 5) 2s 5. (3.7) 


5 eP cee 
+ (1 ee ae +—e| 2 biog? 


i=n+1 


1 —2n = 
ae LESsigg |) Gg 07") G5" (n 22). J 


Fiir vorgegebene Werte von g und ¢ kann dieses Gleichungssystem leicht iterativ gelést werden. 


4. Anwendung der Stérungsrechnung. Bekanntlich verandert sich der maximale Wert des 
Biegemoments in der Platte mit Q) und e. Um die Abhangigkeit des maximalen Momentes von @ 
und ¢ zu finden, ist es indessen bequemer, zur Bestimmung der Koeffizienten die Stérungsrechnung 
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zu verwenden, in der 0) und ¢ die Stérungsparameter sind. Wenn Q als gro angenommen ist, 
so mégen wird die Ansatze 
00 
Cs y a®) Re +r), b, = yy b) pe 1) 
(4.1) 


1 = 


b,, = 5 b”) go ae (n = 2) 2 Kise Real Ona 


machen. Setzen wir dies in (3.2), (3.3), (3.6) ein und vergleichen die Koeffizienten der gleichen 
Potenzen von (9, so erhalten wir 


4 =—], x) = — 5), PAC) = bt) ae bi) = 6), x3) — bt?) LY bt?) b(°) b@) } 25), 
t t' t” : 
SAO) SS ome be) a plr-2) ee ptr) ey bias + 2) ae 22 bys 1) 4 pean (r 2 4) : (4.2) 
wobei 
1 
ee (r + 2) s zs Pa gerade 
= : , U=47 = 45 si wenn) b= ist; 
>: (Gael) > (peels oy (r —1) ungerade 

ferner 

u 1 1 U I EY Ae 1 0 1 0 
= a() =1, ve al== oe = B09) + b®) Sa a?) a a(!) | a al) b) | 60°) | bO) bt Lr 4 . 
1 1 1 os < a: a 
— al) = = a1) oe (av-}) a al) ee bt 1) ! b¢ 2) | bf 1) be 2) (r = 3) : 

1 1 re z 
= 6) = — al) 4 x, - 6) = — al) + af) + = bY) + xe) aa (4.4) 
A ye Mir Oey ORES Weer eM ES aT. 

Te nm n mM nm =~ 4 

NE es) aa ere ta ame ALI CEC) 

ne ™ ee ne ” ne ™ n n n+1 

1 ete ea 1 >2) (4.5) 
OR ae ea (n) @) 4 pil) 4 BA). — AO) (nz 
ne a Ne a, I (n a 1) € Ge bi | bt | er bY), ? 

1 1 
Seas (ni) ee (r—n) | ee 3 (r—1) ____ q(r—n—2)) ___ f(r—n + 1) (r—n) (r=n) __ f(r—n—1) 
ne a, ne a, 1 (n + 1) € (ayn eal ) bt + bY i Le oe 
en Poe 

ae OO) A) =a a OR OL LO) 

Je 2 Sa Ue 9 Dis SS sates ieee y) 
Th = 0) tae) gt) ee hee) eee) > 9 ge 
pela ees al ak rae Se a0. )earinants ieezee2) 
i fh) — Fal) i bY = a@) a) qa) + : b@) 
ie 1” a a ee ae m1 nT (n==1) 62” 
i bt”) = al’) ae) at-}) | a(’—2) ied dy bo), (r ZAG 2), (n = 4) ; 
es nm n (n—l1)eé : (n=>3) (4.6b) 
2=f(n—1) — (m=1) —— g(r—2) __ (n—2) _| (n—3) Sai h(n) ar (0) 
ne bt Ur oer a, iz ane a (n ea 1) & (OM, b”) ? 

1 1 1 
ak (r) 7 en ye) eee (r—1) ((7==2) i eee (r—n) uf (r) (r—n +1) > 

ne bY Oe ei a, or a, | ne bt | (n —_ 1) & (CONS bY ) (r = n) ss J 


Mit diesen Gleichungen kann man die Koeffizienten a, b©, x@ der Reihe (4.1) nacheinander 
bestimmen. 

Ferner erkennen wir aus diesen Gleichungen, daB die Koeffizienten a) und b© als Polynom 
von €, dessen Grad 2r + 1 bzw. 27 + 2 ist, dargestellt werden. 

Die vorhergehende Lésung ist nur formal, und um die Giiltigkeit der Lésung zu sichern, wire 
es nétig, die Konvergenz der Reihe (4.1) zu beweisen. Der direkte Beweis ist ziemlich schwer, 
aber das Ergebnis der numerischen Rechnungen zeigt, daB die Reihe (4.1) bis zum siebenten 
Gliede einschlieBlich eine praktisch befriedigende Konvergenz fiir 9) > 2 ergibt. 
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Die Werte der Koeffizienten, die in der Reihe (4.1) enthalten sind, haben wir bis zu r = 7 
berechnet, aber der Platz erlaubt leider nicht, die Werte hier in Tabellen beizufiigen. 


5. Entwicklung des maximalen Moments in Reihe nach ge, und €. Fiir die praktische Anwendung 
interessieren uns naturgema die Biegemomentskomponenten, und hiervon wieder die GréBtwerte, 
die aus Symmetriegriinden am Lochrand der engsten Stelle erwartet werden diirfen. 


Da am Lochrand M, = 0 ist, so ergeben sich die maximalen Biegemomente 


IY Revers 1 ss ns ee —(n 2 
My, mt +( pee 1) 2 Fos 1) (a, Qo oe b,, 05! ')) ( + Ns e 


L 


Nimmt man die Ansatze (4.1) hinzu, so kann das gréBte Biegemoment in eine Reihe von 0, und € 
entwickelt werden: 


Menvax/Mo 


Mmax 
My 


1 ais ' 
te Oe aa) ree e” 0 5° (5.1) 
2 r=2 p=1 


wobei die Grenze t folgende Werte hat 


= [r a J) wenn, 7 == pened | ist . 
. | lungerade| 


Die berechneten Werte von m,_,, sind in Tafel 1 an- 
gegeben. In Abb. 2 ist M,,,. 


in Abhangigkeit von 
é (oder v) und r/d aufgetragen. 


Im Grenzfall r/d—> 0 (0) > co), d. h. fiir die un- 
endlich breite Platte mit einem Kreisloch erhalt 
man mit 


=2—e=27 (5.2) 


das wohlbekannte Resultat!. 


Es ist leicht zu zeigen, das in der von beiden 


Begrenzungen gebildeten Verengung ein Biege- 
pale: eae rigee till gat Ful’ Sar ogy G20 025 030 
moment —ax M, iibertragen wird; dies gibt ein : 
durchschnittliches Moment Abb. 2. Maximales Biegemoment fiir reine Biegung 
um y-Achse. 
Min ae) 3 x” (5 3) 
My Gene ae 


Tabelle 1. Zahlenwerte der Koeffizienten m,, p 


24 2 na) oA eS eat ass | 8 | 9 sO a ea eae Seg 15 
| | | | | | | | | 
1 ly Ble 2) |) 0 10 ne OMe Oi ee Os 12 | 32 
2 Zale GN S TOP ==96,1 2-9 1 59 66 24119) = 129) 64) 184) 462 
3 3| 10| 3 | —44|—73 | 78| 331) 138)/—622)/—852) 104) 1128] 1591 1298 
4 | et tahet COr 1) eelS 1° 52s) eed 20) 16 F864) wit25s\ ) = 16i-—2370 320 
5 | 70: | 30 2 —100 40. 18 |—364| 430| 2506 | —1122 | —6502 5194 
6 | | eae 50 | 38|—106|—237|—276| 214] 1776) 314) —5304 
q | 27] 90] —38|—410|—101| 1012] ~ 451 |—3142|—4360| 10802 
8 | | | 27] 120] 131] 398] 256! —948| 1346| 8774 
9 | | 81| 270] 389| 834] 2218] 3440! 321) —5624 
10 | | 54) 270| 1101) 2354) —721| —8712 
11 | | 243/ 810) 897; 380) 3106) 13402 
12 | | | | 81) 540) 3211) 7500 
13 | | | | | | | 729] 2430) 1260) —4152 
14 | | | | | | 0; 810 
LS | te | | | | | $2187 |) 7290 


1 H, Neuber, Z. angew. Math. Mech. 20 (1940) S. 199. 
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und damit wird der Koeffizient der Spannungserhéhung 
Vilexace 


ead 


Die Abhangigkeit des Koeffizienten der Spannungserhéhung von r/d und y zeigt Abb. 3. 


15 02 04 06 08 ” 40 


Abb. 3. Koeffizient der Spannungserhéhung kp in Abhangigkeit von dem Kerbmab (r/d). 
Zugleich Lésung fiir biegebelasteten Streifen mit einem zentralen Kreisloch (vy = 0,3). 
Rechnung. ooo000 Lésung fiir gelochten Streifen. 


Nach der 6rtlichen Eigenschaft der Spannungsverteilung darf man erwarten, daB der erhaltene 
Koeffizient zugleich auch naiherungsweise die Kerbwirkung eines biegebelasteten Streifens mit 
einem zentralen Kreisloch wiedergibt. Diese Erwartung wird durch Vergleich mit den Ergebnissen 
der Rechnung bestatigt, die der Verfasser! fiir solche Streifen durchgefiihrt hat (Abb. 3). 

Meinem verehrten Lehrer, Herrn Prof. Emeritus Dr. S. Higuchi méchte ich an dieser Stelle 
fiir sein dauerndes Interesse an dieser Arbeit und fiir viele Ratschlage herzlich danken. 


(Eingegangen am 12. August 1957.) 


Anschrift des Verfassers: Professor Dr. Osamu Tamate, Sakurakoji, Sendai (Japan), Tohoku University. 


1 QO. Tainate, Trans. Japan Soc. Mech. Engineers 20 (1954) S. 687. 


XXVI. Band 1958 HH. Lippmann: Eine auf Kristallplastizitat beruhende mathematische Plastizitatstheorie 187 


Begriindung einer auf Kristallplastizitat beruhenden mathematischen 
Plastizitatstheorie 


Von H. Lippmann! 


1. Einleitung. Die mathematische Behandlung der Mechanik deformierbarer Medien stiitzt 
sich im allgemeinen auf den Kraft-, Spannungs- und Verformungsbegriff, angewendet auf Volumen- 
elemente. Das Newtonsche Prinzip von Wirkung gleich Gegenwirkung fordert eine gewisse Ein- 
schrankung des Systems der Spannungen und Krafte, welche sich durch die sog. Bedingung des 
Spannungsgleichgewichts [siehe Gleichung (3)] ausdriickt. Dariiber hinaus setzt man — je nach 
der Eigenart des betreffenden Mediums — eine Beziehung zwischen Spannungen und Verformungen 
voraus (,,rheologisches Grundgesetz“, ,,Stoffgesetz‘*) und versucht, simtliche entstandenen Glei- 
chungen unter Wahrung der Randbedingungen zu lésen. 

Die verschiedenen in der mathematischen Plastizitatstheorie gebrauchlichen Formen des Stoff- 
gesetzes sind phanomenologischer Art®?. Sie gehen in der allgemeinen Form auf Hencky* und 
v. Mises® zuriick. Und zwar iibertragt man die im ZerreifSversuch erwiesene Streckgrenzeneigen- 
schaft von Stahl [Verformung Null fiir geringe Zugspannung, Verformung zeitabhangig unter Bei- 
behaltung einer konstanten Schwellenzugspannung ky (.,Streckgrenze“‘, ,,Formanderungsfestig- 
keit™), spater Spannungserhéhung (,,Verfestigung“)| auf dreidimensionale Verformung eines Vo- 
lumenelementes. Den Weg zur exakten Begriindung dieser Ubertragung wies kiirzlich Ziegler’. 
Es sind Zusatzannahmen erforderlich. 

Demgegeniiber gehen wir vom physikalisch bewahrten (und bei physikalischen Untersuchungen 
allein zu Grunde gelegten) Kristallgittermodell aus und schlagen ein darauf beruhendes Stoffgesetz 
vor (Abschnitt 3). Wir versuchen, hierdurch einem Programm zu entsprechen, welches in jiingster 
Zeit wiederholt gestellt wurde’®. Wir kénnen verschiedene allgemeine SchluBfolgerungen 
ziehen, welche bisher in die Ansatze eingingen (Abschnitt 4). Bei den Anwendungen in den Ab- 
schnitten 5 und 6 beschrainken wir uns auf eindimensionale Spannungszustaénde und finden u. a. 
eine zwanglose Erklarung fiir die viel umstrittene Geschwindigkeitsabhangigkeit der Streckgrenze 
(Formanderungsfestigkeit), welche durch die Art der Formgebung — evtl. entscheidend — beein- 
fluBt wird. 


Zu Beginn prazisieren wir unsere Voraussetzungen. 


2. Bezeichnungen und Voraussetzungen. a) Die Punkte eines Metallkérpers K homogener 
Werkstoffzusammensetzung seien durch (i. a. krummlinige) orthogonale Koordinaten u = (u', u?, u?) 
dargestellt. Sie mégen dem Spannungszustand o = o;;, (i, k = 1, 2,3) unterworfen und im (von 
der Zeit t = 0 abhangigen) Gleichgewicht sein, in das sie aus dem spannungsfreien Ausgangszustand 
durch die Verschiebungen x‘(c, u,t) gelangten. Mit e = e'* bezeichnen wir den wie iiblich defi- 
nierten Verzerrungstensor. Auf die Volumenelemente von K wirke die Volumenkraft %, = %. 
Uber gleichartige ungleichstandige Indizes werde summiert. 

b) Wir machen folgende Voraussetzungen zur Herleitung des Stoffgesetzes: 

I. K setzt sich aus sehr vielen Kristallkérnern zusammen, deren Durchmesser klein gegen die 
relative Spannungsanderung sind. 

Als Faustregel gelte: Keine (oder wenige) Spannungssingularitaten. 

Exakte Formulierung: Stellt 4, das Symbol einer Anderung lings der Richtung r, d,s das 
Langendifferential in der betreffenden Richtung, D,,,, den maximalen Korndurchmesser und o,, 


1 Diese Arbeit entstand im Forschungsinstitut fiir bildsame Formung der Metalle, Zwickau, dem der 
Verfasser hierdurch aufrichtig dankt. Sie wurde auszugsweise auf der Gamm-Tagung in Hamburg 1957 vor- 
getragen. y 
R. Hill, The Mathematical Theory of Plasticity, Oxford 1950, 

W. Prager und P.G. Hodge, Theorie ideal-plastischer Korper, Wien 1954. y 

H. Hencky, Z. angew. Math. Mech. 4 (1924) S. 323 und Z. angew. Math. Mech. 5 (1925) S. 115. 

5 R.v. Mises, Z. angew. Math. Mech. 8 (1928) S. 161. 

6 H. Ziegler, Ing.-Arch. 25 (1957) S. 58; siehe auch D. C. Drucker, Quart. Appl. Math. 8 (1950) S. 411. 
7? Verformung und FlieBen des Festkérpers, Kolloquium Madrid 1955, R. Grammel, 5. 1, Berlin/Gottingen/ 


Heidelberg 1956. : 
© iy before und FlieBen der Festkérpers, A. Seeger, S. 90, Berlin/Gottingen/Heidelberg 1950. 
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eine bestimmte Komponente des Spannungstensors dar, so soll 


A Ont | BB d,s 
i ene eS 


gelten. 

‘ II. Herrscht in der Richtung r einer Gleitebene EF’ eines unverspannten Kristalls die Schub- 
spannung o, und die Normalspannung o,,, so gleiten die durch E getrennten Kristallteile in Rich- 
tung r abhingig von o,, o, und ¢ aufeinander ab. Wachsende o, und t kénnen den Betrag der Ab- 
gleitung nicht verkleinern. Die Spannung o, = 0 bewirkt keine Abgleitung. 

III. Innere Verspannungen im (viel-)kristallinen Gefiige erhéhen die zur Erreichung einer be- 
stimmten Abgleitung erforderliche Schubspannung o, gegeniiber dem unverspannten Gitter um 
einen ortsabhangigen, zeitunabhangigen Faktor m. 

IV. Die Abgleitungen der einzelnen Ebenen summieren sich proportional deren flachenmaBigen 
Anteilen zur Gesamtverformung, doch treten infolge von FlieSbehinderungen im Vielkristall gegen- 
iiber dem Einkristall nur die A(< 1)-fachen Betrage in Erscheinung. Der Faktor / ist wie m orts-, 
doch nicht zeitabhangig. 

V. Es existieren eine (von Ort zu Ort stetig veranderliche) Verteilungsdichte f(n, u) der Flachen- 
anteile aller Gleitebenen FE im Punkte u, dargestellt als Funktion der Normaleneinsvektoren n = ni 
von FE, und eine ortliche Dichte M(u) des Flacheninhaltes samtlicher Gleitebenenflachen von K. 
Als Normierung soll gelten 


(0) 
if fener) dowels, 


wobei O die Oberflache der Einskugel und do das zugehérige Oberflachenelement bedeuten, ferner 


(K) e 
f M(u) dV =M, 
wobei M die Gesamtflache aller Gleitebenen in K und dV das Volumenelement ist. 

c) Wir erlautern und erganzen Voraussetzung I.—V. wie folgt: 

Zu I. Eine zu (1) analoge Forderung ist stets dann nétig, wenn diskrete Erscheinungen durch 
stetige (oder differentierbare usw.) Funktionen oder statistisch angendhert werden sollen. Wah- 
rend sie meist stillschweigend benutzt wird, formulieren wie sie hier explizit. Auf Grund ihrer 
setzen wir den Spannungstensor értlich stetig differentiierbar voraus. Die Voraussetzung I ist 
Voraussetzung fiir V. 

Zu I. Es gibt eine int und o, monoton nicht fallende Funktion g(¢, o,, ¢,,), welche g(t, 0, 0,) = 9 
erfillt und die Abgleitung im unverspannten Einkristall beschreibt. Auf Grund der Kristall- 
plastizitatstheorie ist sie proportional zur Zahl der sog. Versetzungsquellen langs der Gleitrichtung, 
also wegen der homogenen Werkstoffzusammensetzung proportional zur Lange der Ebene in Gleit- 
richtung. Betragt diese d,s, die Abgleitung d,u (die Annaherung durch Differentiale rechtfertigt 
sich wegen I), so setzen wir 

du} — g(t, 0, On) |d,s| » (2) 
Verschiedene Ansatze fiir g untersuchen wir in den Abschnitten 5 und 6. 


Zu Milf. Im Sinne der Kristallplastizitat fassen wir ,,Abgleitung® als Bildung und nachfolgende 
Auflésung einer sog. ,,Versetzung“ auft. Der Gitterverzerrung und -verspannung durch gebildete, 
aber nicht aufgeléste Versetzungen messen wir keine makroskopische Wirkung bei, da durch sie 
infinitesimale Volumenanderungen auftreten, die sich im allgemeinen makroskopisch bemerkbar 
machen miiBten, aber noch nicht beobachtet wurden?. Daran andert nichts, daB fiir geeignete 
Korper die Verformung allein auf Grund des Versetzungsmodells ohne meBbare Volumenanderung 
beschrieben wurde (Balkenbiegung* 4.) Hingegen verursachen vorhandene Versetzungen, ins- 
besondere Verdrehung der Kérner gegeneinander usw. eine Verspannung, welche eine Erhéhung 
der kritischen Schubspannung k (die ttberhaupt erst merkbare Verformung des Kristalls bewirkt) 
um einen Faktor m verursacht. Sachs> und Taylor® errechnen fiir kubisch-flachenzentrierte Metalle 


2,200 =m = 3.00, 


+ P. Haasen und G. Leibfried, Fortschr. d. Phys. 2 (1954) S. 73. 

E. Siebel, Die Formgebung im bildsamen Zustande, Diisseldorf 1932. 

J. F. Nye, Acta met. 1 (1953) S. 153—162. 

W. T. Read, Acta met. 5 (1957) S. 83. 

5 G. Sachs, Z. VDI 72 (1928) S. 734. 

° Verformung und FlieBen des Festkérpers, G. I. Taylor, S. 3, Berlin/Gottingen/Heidelberg 1956. 
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der erstgenannte durch statistische Mittelung iiber die Kornorientierungen. Wir erweitern die 
Bedeutung von m auf den ganzen (nicht nur den kritischen) Spannungszustand, vernachlassigen 
jedoch eine Verformungsabhingigkeit, weil wir die Neubildung von Versetzungen vernachlassigen. 

Zu IV. Die Ortsabhangigkeit von J untersuchten (mit anderer Terminologie) z. B. Vitovec und 
Slibar!»?, Nach ihnen ist A an der Kérperoberflache gréBer als im Innern, was sich jedoch bei ver- 
ninftigen K6rpern mit relativ kleiner Oberflache technisch wenig bemerkbar macht. Verformungs- 
unabhangigkeit bewies Kochendérfer® fiir kubisch-flachenzentrierte Metalle. Bei der technisch 
wichtigen Warmformgebung diirfte sie infolge der Rekristallisation allgemein gelten. 

Die GréBe des Beitrages der betrachteten Ebenenrichtung E zur Gesamtabgleitung ergibt sich 
proportional zur Zahl der Versetzungsquellen in einem beliebig breiten Streifen der Ebenen E in 
Richtung r, also (wegen der Unabhangigkeit von der Streifenbreite) proportional zu dem Produkt 

Ebene Versetzungsdichte mal Streifenlinge 
~f(n,u) mal Streifenlange , 
welches beim Vielkristall in (2) fiir d,s einzusetzen ist. 

Zu V. Stetigkeit von f und M bringt lediglich mathematische Vorteile; jedenfalls kénnen auf 
Grund von I die GréBen f und M durch stetige Funktionen angenahert werden. Die Verformungs- 
unabhangigkeit rechtfertigt sich wenigstens bei Warmformgebung durch Rekristallisation. 

d) Neben dem Stoffgesetz besteht die oben erwahnte Bedingung des Spannungsgleichgewichtes, 
welche sich auf Grund der stetigen Differentiierbarkeit von o bekanntlich durch 

Div o=% (3) 
ausdriickt. 

3. Das Stoffgesetz. a) Eine gewisse Gleitebene EF im Punkte u besitze den Normalen-Einsvektor 


n —n'‘. Orthogonal zu diesem, also in E, denken wir die orthogonalen Einsvektoren v = v', 
w —=w’'. Der Einsvektor r = r‘ in EF besitzt die Darstellung 


r=vcosfh +wsinf, 


wobei f der Winkel zwischen r und v ist und 0 =f < 2z gilt. Wir betrachten die Schubspannung o,, 
welche der der Spitze von n zugewandte Teil des durch FE geteilten Kristallkornes auf den anderen 


Teil ausiibt: 


k 


G76, nr =o, nv" cos B +o,, n*-w* sin'B:, 


Es gibt eine Richtung r = r’, in welcher o, verschwindet. Gehért /’ zu r’, so findet man namlich 


Cik n* vk 
tg 7 — _~ ee 
gp Cit ni wk 
fiir o,,n' w* ~0 und fp’ = = fiir o,, n' w* = 0. Hieraus folgt r' = w im zweiten Fall und 


ca (oi, n* wk) v — (oj, n* VE) w 

oe 

ae V(Giz ni vb)? + (cin ni wh)? 

im ersten Fall. Es bezeichne r’’ einen der in FE gelegenen zu r’ orthogonalen Einsvektoren, fiir den 
oir né r”* > 0. Man kann dann fiir nicht verschwindende Nenner eindeutig 


Cin n* vk ye é 
fle er Os oF vy | ——__—_—_ —————— w, (4) 
Vig ni vj? + (cin, né wh)? Vik nt v*)? + (cin ni w*)? 
andernfalls etwa r'’ = v setzen. 
Wir wahlen voriibergehend v = r” sowie w = r’ und folgern 


0, = 0:5 ni v* cosf, 

also eine in f symmetrische Schubspannungsverteilung. . 
Bei Summation aller zugehérigen Abgleitungen heben sich die zu v orthogonalen Komponenten 
auf. Nur in Richtung r =r = », also in der Richtung maximaler Schubspannung, bleibt Ab- 


gleitung iibrig. Wahlen wir v und w wieder allgemein, so errechnet sich nach (4) 


Gor ni r"® — Vor. n' v*)2 + (6, n' w*)? . (5) 
r', r” und (5) sind unabhangig von den speziellen Vektoren v und w. 
A, Slibar und F. Vitovec, Schweiz. Arch. f. angew. Wiss. u. Technik 16 (1950) S. 76. 


F. Vitovec und A. Slibar, Schweiz. Arch. f. angew. Wiss. u. Technik 16 (1950) H. 3. ; 
A. Kochendérfer, Plastische Eigenschaften von Kristallen und metallischen Werkstoffen, Berlin 1941. 
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b) Wir wahlen als Achsenrichtung 1, 2, 3 eines lokalen Koordinatensystems die Richtungen n, 
r”, r’ und bezeichnen darinnen mit ¢ denjenigen Verzerrungstensor, welcher sich ausschlieBlich 
auf die Abgleitung in E bezieht. Stimmt E mit einer Gleitebene tiberein, so erkennt man aus der 
Betrachtung eines im betreffenden Korn gelegenen hinreichend kleinen Volumenelementes (Kanten- 
lange dv, Abb. 1), welches durch E halbiert wird, mittels (2) 


H4 m 
gt? = 0) fiir alle anderen Paare @; 0 (aa losda od ) 


Ein gréBeres, viele Kérner enthaltendes Volumenelement gleicher Lage verformt sich durch Ab- 
gleitung auf vielen parallel zu E gelegenen Ebenestiicken nahezu parallelepipedisch (ausgezogene 
schrage Linien in Abb. 1) nach demselben Tensor, jedoch proportional zum Anteil der Ebenen- 
stiicken parallel zu E im Verhaltnis zu den iibrigen. Nach 2c) (zu IV) miissen wir somit 

mn gtk — M feik 

betrachten. 

Bezeichnet 7‘ die Koordinaten irgend 
eines Vektors im jetzt zugrunde gelegten 
System, y' die Koordinaten im urspriing- 
lichen System, so gilt 


1 A Onn 
2 = ol = — te x = 8[F ; 10) 


dy 


dv 


y* eae Bane 


Born [3s es Bere 
Die Koordinaten e“ von é* im urspriing- 
lichen System lauten folglich 
a — Bi Bhat’, 

das heibt 

ek pas (ni rk alts ri n*) ele, 
Durchlauft n bei Summation iiber O alle Richtungen, so tritt jede Gleitebene genau zweimal auf; 
das Ergebnis ist also zu halbieren: 


Abb, 1. Abgleitungsmodell. 


a ae | f(n, u) g ¢ a ? °,] (né ye of ra n*) do ? 
oder wegen o, = 0;, n' n*, (4) und (5) 


(0) a 
g (« u V@jrni v!)? + (ojini w!)?, of. ni n' 
m 


on MOH Lain u) —— 


V(oji ni vl)? + (ojrni wl? 


X [(oj, n4 v') (nt v* + nk v*) + (oj, 0) w') (nr? w* + n*w')] ¢, (6) 
wo O den Teil von O darstellt, auf dem der Nenner nicht verschwindet. 

c) Zwar ist (6) unabhangig von v und w, doch treten v und w formal auf und miiBten bei An- 
wendungen ausgerechnet werden. Wir suchen eine Darstellung, in der das nicht mehr der Fall ist. 
Hierzu wahlen wir ein Hauptspannungssystem (c,, = 0 fiir a ~ b) als Koordinatensystem und 
beachten die Orthogonalitat der aus n, v, w gebildeten Matrix: 

(v?)? + (w*)? + (n2)? = 1, vv? + ww + n? n> — 0 fiir Chee li%, 
Hieraus folgt 
(oj, nF vo!) + (oj, nJ w')? 
= (04, n* 0} + Og 1 v2 + Ogg n3 v3)? + (oy, nw! + oo9 0? Ww? + oy n3 w3)? 
= (0, 0°)? [1 — (n')?] 
tH (Ogg n*)? [L — (n®)?] + (693 03)? [1 — (n3)?] — 2 oy; og9 (nr n?)? 
— 2 Og9 Ogg (n? n®)? — 2 d53 oy (n? n)? 
= — [ox (nt)? P — [oga(n*)? PP — [og9(n8)?]? — 2 o4,(n3)? og9(n?)? 
— 2 (1)? o59(n?)? — 2 6y9(n?)? o5,(n3)? 


= — [on(n!? + og9(n*)® + o5(n*)?]? + (Oy m4)? + (Gag n®)? + (G39 03)? 
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Driickt man dies in einer fiir allgemeine orthogonale Koordinaten legitimen Form aus, so sieht 
man sofort 
ae orn sem 12 feel 
(oj, nv’? + (oj, nr w!)? = > (6, n')? — (oj, nin!) (7) 
j 
Im Hauptspannungssystem selbst erhalt man eine symmetrische Form, wenn man nach dem 


zweiten Gleichheitszeichen der obigen Umrechnung ersetzt 1 — (n‘)? = (né+1)2 +- (n'+2)2 (Exp. 
mod 3): 


(On — Og9)” (n+ n)? + (Gag — Ogg) (n? n3)? + (033 — Oy) (n? n})? (8) 
SchlieBlich muB in (6) die eckige Klammer umgerechnet werden. Hierzu beachten wir 


(o;,, n' v*) (n* v? + n® v*) + (6,, n* w*) (n? w® + n? w*) 


= s & as : (oi ni 2 4 (Fix n' w)?] . (9) 


20ab | Oba 


Setzen wir fiir die eckige Klammer der rechten Seite von (9) die rechte Seite von (7) ein, so gewinnen 
wir die umgeformte linke Seite und schliefSlich 


(0) 


1 = re 
% g(t > (oj, n!)? — (oj, n' n!)?, oj) nt n!) 
n, w) ————————— a al 
( y= (oj1n )? — (oj, nt n!)? 
j 
x [nt (o,,n') + né (o,,n') — 2 n' né (oj, n) n')] do, (10) 
wo wie bei (6) Oc O durch 


in _ M(u) Au) 
hy ais 


Dy (oj, n')? — (6, n n')? > 0 
j 
charakterisiert wird. 


Gleichung (6) bzw. (10) ist die allgemeinste Form des Stoffgesetzes. 


4, Allgemeine Folgerungen. a) » e'' = 0 bedeutet (infinitesimale) Volumenkonstanz. Sie war 
auf Grund von 2 c) (zu IIT) zu erwarten und folgt aus »’é'! = 0: Orthogonale Riicktransformation 
liefert Y e** — 0, Integration die Behauptung. 

b) Hydrostatischer Druck o;; = p, o;, = 0 fir 1 Ak bewirkt keine Verformung. Dies folgt 
aus (10) wegen 

Dy (G1 n')? — (Fj nin’? = 0. 
j 

c) G4, = 9 fir alle a, b mit a~ b bewirkt e*° = 0: Die Achsenrichtungen des Spannungs- und 
Verzerrungsellipsoides fallen zusammen. Wir betrachten e*°(a ~ b) und setzen c gleich der unter 
a, b nicht enthaltenen Zahl aus (1, 2, 3). Alsdann zerlegen wir die Summation tiber O in Summation 
ber die beiden Teile 0,, O,, welche durch die (a, c)-Ebene getrennt werden, und gehen von (6) aus. 
Einem gewissen Dreibein (n, v, w), fiir welches die Spitze von n auf 0, liegt, ordnen wir eineindeutig 
ein zweites orthonormales Dreibein (n, v, w) mit der Spitze von n auf O, zu durch 


if Ne, n’ ——n?, NT! 5 
vy? =v", vy? = —v?, v=v, 
w* = —w’, ww, we =—w 
Es folgt 
o,,niv' =o,,niv', 0, w' =—o,,niw', a, nn! =o,,nIn', 
n°? vp? ae n? = (n* pert n? v") p n* w ee n° wt = n°? we ge n® w F 


(o;,niv')? + (oj; n/ w')?P = (o;, 0! v')? + (o;,n) w')?. 


Zugeordnete Integranden besitzen folglich entgegengesetztes Vorzeichen: Summationen iiber O, 


und OQ, heben sich auf. 
d) Ist g(t, o,,0,) = g(t,¢,) nicht von o, abhangig, so existiert eine nach o;, differentiierbare 
Funktion Fu, o) (,,plastisches Potential) mit 


ik — a F(u, 0) - (11) 


13* 
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Die Funktion F besitzt die unten angegebene Gestalt (12). Die genannte Bedingung ist unter 
gewissen Zusatzannahmen auch notwendig, z. B.: g(t, ¢,,0,) mége fiir wachsende o,, (zunehmende 
Reibung) nicht wachsen!. 

Wenn namlich g(t, o,,¢,) = g(t, ,), so erfiillt wegen (9) und (11) 


s 


(0 


) s s i 
F(u, 0) =e) [se de ic 5 doer | a(t me (12) 
0 0 


die Behauptung mit s ==) 2) (o;,, n*)? — (a, n' n*)? als obere Grenze. 
j 


Es werde nun (11) unter der gemachten Zusatzannahme vorausgesetzt, und es gelte etwa fiir 
6, <6, und bestimmtes o,: g(t, 0,,0,) < 8(t, 0,0’). Dann betrachten wir eine geschlossene, 
stetige, nicht zu einem Punkt entartete Kurve U im Raum der drei Hauptspannungen (als Unter- 
raum des neundimensionalen Spannungsraumes), welche einen Punkt P besitzt, in dem 


[2 Cunt? — (yn! nk? =0,,  o4,n'nk =a, 
i 


gilt. Die erste Gleichung ist nach (8) durch entsprechend groBe Hauptspannungsdifferenzen er- 
fiillbar, die zweite durch Addition einer gewissen Konstanten zu allen o;,. In P mége ferner U 
eine Strecke parallel der o,,-Achse durchlaufen derart, daB do,, > 0. Durchlauft o die Kurve U, so 
beschreibt 1; = 0; + (co, —o,,) 6; in einer Umgebung von P, wo do. = doz; = 90, Ti, = Oip, 
sonst — eine Kurve U’. Wegen (11) 

U 
Feit dG, == 0. 


Dann gilt wegen (8) und (10), da sich nur die Normalspannung andert, 

(YY) 

fe* dz, = 0 
im Widerspruch zu (11). Zwar ist U' an zwei Punkten unstetig und nicht zusammenhangend, 
kann aber ersichtlich ohne Anderung des Ergebnisses zu einer geschlossenen stetigen Kurve ver- 
vollstandigt werden. 

e) Ist g unabhangig von der Normalspannung ¢,,, so ist e‘* unabhangig von einem zusatzlichen 
hydrostatischen Druck (Erganzung zu b)?. Man erkennt dies nach (8) und (9) (Umrechnung der 
geschweiften Klammer von (6) bzw. (10) auf Hauptspannungen) im Hauptspannungssystem. 

f) Die Abschnitte a) bis c) enthielten Folgerungen, die wir auf Grund unseres Modells beweisen 
konnten. In der phénomenologischen Plastizitatstheorie dienten sie als Axiome. Lediglich die 
Einfiihrung des plastischen Potentials bedurfte einer Hypothese, welche naheliegt und etwa durch 
Abschnitt e) gerechtfertigt wird. Dies ist wiederum Axiom der gebrauchlichen Theorie. Gleichung 
(11) entspricht der Musesschen Form des Stoffgesetzes und unterscheidet sich von dieser 
erstens durch das Fehlen eines ortsabhangigen Proportionalitatsfaktors, dafiir Ortsabhangigkeit 
von F selbst (welche jedoch ebenfalls schon betrachtet wurde, vgl. Olszak*), und zweitens durch 
das Fehlen einer ,,FlieSbedingung* EF = konst. = ky (vgl. jedoch Abschnitt 6 a). 


5. Theorie der Drahtdehnung. a) Die physikalischen Anwendungen unserer Theorie unter- 
scheiden sich von den technischen des folgenden Abschnittes durch die Wahl von g. 

Der Verfasser* entwickelte auf Grund einer Hypothese von Engelhardt® iiber die Grundlagen der 
Drahtdehnung Formeln fiir die sog. ,,spontane Dehnung“ y(t) (Dehnung, die sich kurzzeitig nach 
Belastung auswirkt) und das sog. ,,Kriechen“ d(z,t) (7 die Zugspannung, t die Zeit) 


(2) =A f fla) de, (13) 
alt alt 
d(t,t) = A JH) dx — TF) Ga a) dae (14) 


1 Eine andere derartige Zusatzannahme: (11) moge fiir alle Verteilungsdichten f(n, u) gelten. 
. - Auch hier ist die Umkehrung unter gewisse Zusatzannahmen richtig, z. B. der in der vorangehenden 
uBnote. 
° W. Olszak, Elastisch-plastische Biegung von nicht-homogenen orthotropen Kreisplatten. Vortrag Gamm- 
Tagung in Hamburg 1957. 
4 H. Lippmann, Acta met. 4 (1956) S. 298. 
° W. Engelhardt, Z. Metallkde. 42 (1951) S. 117. 
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wo A und a gewisse positive Konstanten, f(x) eine statistische Verteilungsdichte, meist GauBisch 

F(x) = Be“#—o) (15) 
mit B als Normierungsfaktor, h als PrazisionsmaB, 6 = konst., und w(x) eine monoton fiir O<x<a 
wachsende Funktion mit w(0) = 0, w(a) = 1 darstellt. Die sich ergebenden Kurven werden dis- 


kutiert und gezeichnet (qualitativ in Abb. 2 und 3). Sie widersprechen nicht den Experimenten!. 
Der Ausdruck fiir die gesamte plastische Dehnung lautet 


D(z,t) = d(t,t) + y(t). 
Genaue Durchsicht der erwaihnten Theorie yom hier dargestellten Blickpunkt aus erweist, daB 
sie sich (wie die meisten derartigen) lediglich auf die Abgleitung einer einzigen Gleitebene unter 
der Schubspannung o, = 1 bezieht. 
Fiir uns ergeben sich nach allem bei passender Wahl des Faktors A die Ansatze 


? g=y bzw. gas (16) 


Abb. 2. Spontane Abgleitung. Abb. 3. Zeitabhingige Abgleitung fiir r-constant. 


b) Der eindimensionale Spannungszustand o,, =T, 0,, = 0 fiir a 1 oder b va 1 geniigt der 
Gleichgewichtsbedingung (3) bei Vernachlassigung von Volumenkraften. Er 1aBe sich durch einen 
langgestreckten zylindersymmetrischen Kérper von geringem Querschnitt mit an den Enden auf- 
gezwungener, tiber den Querschnitt konstanter Zug- oder Druckspannung realisieren, wenn der 
Zylindermantel spannungsfrei bleibt. Beispiele: Drahtdehnung, ZerreiBversuch, Stauchversuch 
schlanker zylindrischer Proben ohne Reibung an den PreBflachen. 

c) Wir betrachten ¢ = el! und setzen 

ni = cosy, n? = siny cosy, n = siny sing OSgs22, Ojo.) 


also 
do = siny dg dy, 

VX (6,, 0")? — (64, n' nv)? = 7 |cos yp sin y| , 

n'(o,,n') + ni(o,,n') — 2(n')? (6, ni n*) = 2 1 cos? y sin? y . 
Unser Draht sei isotrop, d. h. f = (4), und homogen, d. h. m, M,A = konst. Es folgt nach (10) 
fel + |cos p sin | , r cost y) eos yj sin? y dy . (17) 

m 

i 
Wir setzen nunmehr Unabhangigkeit von o, voraus und formen weiter um: 


_MA 


e 4 


Dan 


ee” [s(t a5 8im2y)sin2y V1 — cos 2 p dy 


‘14 

= Ms g ¢ z. sin x) sin x yl —cos x dx 
8/2 m 

m/2 


m/2 ba 7 
_ Mi {| a (tea ide py cin x) sin.x {71 cos x 4 V14 cos x }dx. 
~ «By2 8/2 m 
ae 0 m/2 Y 0 


1 G, Miiller, nach persénlichen Mitteilungen. 
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Da tak! . ifn ieee cid tis a 
V1 — cos x 4 V1 4 cos x =) (V1 — cos x -- 1 + cos x% \P = }2 yl sin x 5 
so ist 
n/2 
2 [ (+. jpsinx)sin x V1 + sin x dx (18) 
m ‘ 


6 
unabhangig von ul’. tat 
d) Setzen wir g=d, so erkennt man aus der Positivitét des Integranden und weil g fiir jedes 
feste t und x einer Kurve der Abb. 3 folgt, daB diese qualitativ auch fiir unsere jetzige Dehnung 
richtig ist. Eine ahnliche Entscheidung fiir g =y fallt schwerer. Wir umgehen sie durch numerische 
Auswertung. Dazu setzen wir die Proportionalitatsfaktoren 


ABMiA\% _ogipm=—1 [vgl. (13), (14), (15)] 


20 zi 16h 
r=05 und erhalten bei Benutzung des tabellierten Gaup- 
o schen Fehlerintegrals E(x) sowie mit r = hb 
’ n/2 1 
: rath = elk ing ))sine YI + sin xd. 
45|— 02 T sin x 
, (19) 
of 
Der Wert r =0 entspricht dem bei Lippmann allein 
pee behandelten Fall b=0. Wir verzichten auf exakte 
-0f Ausrechnung fiir verschiedene Werte r, nahern viel- 
7 mehr fiir kleine r linear mittels der Taylor entwick- 
lung an: 
ee e = F(t) + rG(z), 
wobei 
/2 
Eq) pene) ae sin « V1 + sin x dx 
“\¢ sin x i 
6 
m/2 (ae ) Wye, re ss 
, E(t) = [ e \tsinx/ | sin x 1 + sinxdx, 
AT MCR RE OM SY TT BEE) Wi é 


Abb. 4, Theoretische spontane Dehnung eines Metalldrahtes 
in Abhangigkeit von der Spannung t und dem Parameter r. iP 


pee 

Vu 
Eine daraus numerisch ermittelte Kurvenschar findet sich in Abb. 4. Der allgemeinste Fall (be- 
liebige Konstanten) entsteht durch kombinierte lineare Verzerrung beider Achsen. In der Tat 
entsprechen die Kurven qualitativ denjenigen von Abb. 2. 


6. Zerrei®- und Stauchversuch. a) Wird beim Drahtdehnungsexperiment im allgemeinen die 
Spannung konstant gehalten, so beim ZerreiB- und Stauchversuch die Geschwindigkeit. Wir 
befassen uns mit dem eindimensionalen Spannungszustand (5b) und benutzen im wesentlichen 
die Terminologie des vorangegangenen Abschnittes. 

ZerreiB- und Stauchpriifungen finden vorwiegend an Stahl statt, wo y (Abb. 2) zu einem eckigen 
zeitlich verzégerten Verlauf ausartet, wahrend das Kriechen zu vernachlassigen ist: 


© fur ohh, 
*~ | Funktion von ¢ und o, furto = ke 


Die kritische Schubspannung k mége nur vom Werkstoff und dessen Beschaffenheit (Temperatur) 
abhangen (Abschnitt 2 c) zu IIT). In Abb. 2 und 4 wurde t nicht beriicksichtigt. Sie beziehen sich 
auf too oder stationare Vorginge. , 

Wir denken die zweite Zeile unserer neuen Definitionsgleichung von y in eine Potenzreihe nach 
o, entwickelt und brechen diese nach dem ersten nicht-verschwindenden Glied (Exponent N =0) 


' Vel. FuBnote 4 auf Seite 192. 
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ab. Dies ergibt den Ansatz 
0 fiir Grae ke 


g(t, o,) = O, (20) 


nce ( ab fir oo. =k. 


Uns interessiert die Streckgrenze (Formanderungsfestigkeit) ky, d.h. diejenige Zugspannung, 
welche zur Aufrechterhaltung einer gewissen Dehnungsgeschwindigkeit (Stauchgeschwindigkeit) 
ndtig ist. ky wird sich unabhangig von der Verformung, aber abhingig von der Geschwindigkeit 
ergeben. 

Formel (20) ahnelt der phanomenologischen Flie8bedingung (vgl. unsere Bemerkung in 4. f)) 
und gibt ein Analogon zum ideal-plastischen Kérper'. 

b) In (17) ist jetzt nur iiber ein Teilintervall von (0, 2) zu integrieren, dadurch definiert, daB 
dort 


9 
= |cos p sin y| = a jsin2 y| =k 
oder 
a are sin aut SS + -— _ arc sin fg : 
Ae ++ 5 are cin? ™E<y <n — are sin euUE : 
Verformung tritt also erst fiir tT 22mk auf. 
Mit (20), (17) und der Abkiirzung 
¢ = are sin Z (21) 
folgt 
a/2—c aw—e 
MA t \N fi 2 MANNE: 
é = hii) (ne) | | + i |siny |cos y sin y|N** dy , 
c n/2+¢ 
MA ¢ vere + N+2 Nt1 
oi — Ae) (a) | sin\** p cos\** p dy . (22) 
c) Wir berechnen das unbestimmte Integral 
f= f sin’? » cosN** yp dy 
1. Ist N gerade, so gilt 
le IN [ZN ee 
cos’ y = (cos® yp)? = (1 — sin? p)N? = Y (— 1pe( / Jsintey . 
e=0 @ 
also 
N/2 
Nj2 N/2 ( 
=) (— Tae sinNt+?e+? wy cosy dy = JY; (— 1)¢ Woe - aE) Big Pasea te 
e= Q o= 
2. Ist N ungerade, so gilt wie oben 
(N+1)/2 
sinN+1 a ae (— 190 (PY coste 
1 QS @ 
also 
N + 1)/2 
were (— 1)? (' ace 
T=— NED GHD — cosNt+2et2 y , 
Q 


o=0 
Einsetzen der Grenzen c und 7 ae liefert in beiden Fallen Faktoren von der Gestalt 
een 


cos® ¢ — sin® ¢ = (cos c — sinc) 3’ cos’ ¢ sin®—!“ ¢ (R = 3, ungerade) 
4=0 


. 2 . 2 i A sath 
= (cos c — sin c) |cos csin c+ y jcos ¢ sin 
1=0 


—1— R—1—A 
1 ‘al 


A oy 
c + sin’ ¢ cos 


1 Vgl. FuBnote 3 auf Seite 187. 
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Wegen 


| ee ——— és = / ——- 
cose = 1 V1 + fl —sint 2 ¢ : sinc = J 1 yl gin? De 


folgt nach Einsetzen von (21) 


; mk 
cos ¢ — sin ¢ = (cos ¢ — sin c)? =)i— - 
und 
x R—1 
Ro) Bonen cos ¢ — sinc Dn ee 2mk\n Sa ALE bas, 2m k\2\» 
cos® ¢ — sink ¢ = (RDP le | 2 ee fe 1 Z ; 
v= v 


Man beachte hierbei, daB keine Wurzeln auftreten, weil sich diese in der Differenzbildung weg- 
heben. Ferner handelt es sich wegen (3) = 0 fir B> A nur um endliche Summen. Wir erhal- 


Imk- 
me und einem Poly- 


ten aU gs als Endresultat, da®B sich ¢ als Produkt von(, ci ya — 


nom in ee i dayerellt: 


Rik ae 
MAN) By ——— S(N,t) mit t22mk (23) 
¢ Tt 

und | 


fiir gerade N, 


1)? (NV + 1)/2 Neel N+1 
yi v( aoe. mhz +e (“FH +e) So pemme S/— — + ee 
<— N+20+2 i / ( Tt a nee 

ess lt y= 


<p—esyy} et 
fiir ungerade N. 


d) Die umstandlichen Gleichungen (23) lassen sich fiir groBere t asymptotisch annahern. 
Ks ist 


ees Shei wo /X 
Cae Til) hese 2s lie aw ookeele 
= 


2» 
also 
o /N cbs 
Sos oe Pe 
vy=U y Q2y 
entsprechend : | 
So INSET =i 3 | 
bet bined 
Dron eae nang 
oN 2» 


Demnach strebt ¢ fiir jedes feste t und 1 —> co asymptotisch gegen eine Parabel N-ter Ordnung: 
N/2 
N/2 (— y*( 
Q 


WY 8O'6 <3" 
Fee ee a Pumeae fir nS Setade 
a (es 1° (a Ae oa ungerade 


(24) 


nh ee aCe 
| o=0 N+2Q+2 
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Spezialfalle sind 


MA 
Ni O's é~— h(t) = C(t), 
fol Ae T 
seal ace un es eng Me 


Genaue Auswertung nach (23) gibt 


iesim e= Cj p1—22* (1 +=). 


t 


ae ei4 T Conk mk | 3 /mk\* 
N=1, ¢ = Ct) 25 V1 Soller, eee 


e) Die Existenz der Streckgrenze (Formanderungsfestigkeit) bei Stahl besagt: 
FlieBgeschwindigkeit ¢ = konst. = v/L (v = Ziehgeschwindigkeit, L = Augenblickslange ~ An- 
fangslange der Probe — entsprechend bei Stauchung) fiir 7 = kp, also t= 0. Nach (23): h(t) ==) 
h(t) =ct. Dies in (23), (24) usw. eingesetzt, ergibt nach Differentiation und Substitution 
die jeweilige Geschwindigkeitsabhangigkeit von kr. 


a 
| 
= 


ie 
amk 


Abb. 5. Theoretische Geschwindigkeitsabhangigkeit der Streckgrenze fiir N = 0 (Kurve 1) und N = 1 (Kurve 2). 


Mit C= M/c/6 ist in Abb.5 —k;/2mk gegen v/CL aufgetragen. Fiir N=0 existiert eine maxi- 
male Ziehgeschwindigkeit. Die Naherung durch die jeweilige Parabel N-ter Ordnung [nach (24)] 
ist gut. Das Experiment muf iiber die GréBe von N entscheiden. 

Besonders wichtig scheint uns folgende abschlieBende Bemerkung: 

Die Geschwindigkeitsabhangigkeit von ky erhalt man bei konstantem k — bei geschwindigkeits- 
abhangigem k erst recht — unter wesentlicher Benutzung der Kérperform und Verformungsart. 
Sie kann bei Anderung dieser Verhaltnisse durchaus anders als hier ausfallen. Infolgedessen sind 
beispielsweise in dynamischen Freistauchversuchen ermittelte Werte ky nicht fiir das Gesenk- 
schmieden, Strangpressen, Blechbiegen, Walzen usw. verwendbar. Hieraus erklart sich zwanglos 
die schlechte Ubereinstimmung bei den Experimenten verschiedener Autoren. Wir formulieren: 
Die GréBe k, ist nicht nur von Werkstoff und Geschwindigkeit, sondern auch von der Verformungs- 
art und der erstrebten Kérperform abhiangig. 


(Eingegangen am 29. August 1957.) 
Anschrift des Verfassers: Dr. Horst Lippmann, Hannover-Dohren, Am Lindenhofe 37. 
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The Elastic Problem of a Homogeneous Circular Ring Acted Upon 
by Equally Spaced Concentrated Twists of Equal Magnitude * 


By I. Malkin 


1. Introduction. The problem of a circular ring under twist occupies in a certain sense of the 
word a solitary position among the one-dimensional problems of Applied Elasticity. Indeed, it 
would not be simple to find, among the problems of the group just referred to, another one of the 
nature of a paradox like our present problem. Referring even to its most elementary special case 
as represented in [1]!, which is based essentially upon Grammel’s work on the subject (see below), 
we see that a circular ring acted upon by uniformly distributed twist, is actually in bending, 
while the torsional stresses are zero; while the stresses produced in such a ring are thus virtually 
bending stresses, the bending deformation is zero. We will see below that the unusual nature 
of the problem becomes even more emphatic in the non-elementary case, when the twisting action 
is not uniformly distributed along the entire periphery of the circular ring, like in the special case 
treated in [1], but concentrated to a finite number of equal and equally spaced twist moments. 
We arrive at a situation which calls for a thorough and careful study of a seemingly simple problem. 
An attentive investigation is very necessary in order to avoid considerable errors in such practical 
applications, in which the ring under twist is an important element. 

A satisfactory solution of the title problem will be 
derived in the present paper on the basis of Grammel’s 
Ir, theory of the circular ring under twist by concentrated 
moments [2], [3]. In the interest of greater clarity, which 
is always desired and welcomed by practical engineers and 
designers, Grammel’s work will have to be reproduced here 
quite elaborately, before the necessary ultimate conclusions 
can be drawn from it. It should be mentioned in this con- 
nection that the theory developed in [2], [3] was the sti- 
mulating factor for further valuable studies by F. K. G. 


od 


i 


Fig. 1. Ring Acted Upon by Concentrated Equally Fig. 2. Beam in Bending 
Spaced and Equal Twist Moments. by End Couples. 


Odqvist [4], in which the general case is treated of a homogeneous circular ring under loads conti- 
nuously, but in general non-uniformly, distributed along the circumference. It must be stated, ho- 
wever, that the studies carried out in [4], are only the beginning of a full discussion, and it is to 
be hoped that these studies will be continued and completed. Attention must be called furthermore 
to an analysis of the twist problem forming a part of the subject covered in [5]; the results deri- 
ved in the present paper will make necessary some modifications in the analysis just referred to. 

We start with a derivation of the loading action diagram of the ring under a system of con- 
centrated twisting moments equally spaced and of equal magnitude. 


2. Derivation of the Loading Action Diagram. The equations, which determine the loading 
conditions of the elastic system under consideration, can be derived in various ways. The shortest 
and best derivation is given in [2], [3]. We shall follow it essentially in the present Section 2. 

Assume the ring to be acted upon by n (n =2) twisting moments equally spaced and of equal 
magnitude T (fig. 1); the spacing angle of the twists is 2F =2z/n. 

It is already at this point important to recognize the general type of deformation, which the 
ring will undergo under the action of the n twists T; G=1 to n). Unjustified analogy with the 


* The present paper is dedicated by its author to the memory of his teacher Professor Dr. Richard von Mises 
(April 19, 1883 — July 14, 1953). 
* Numbers in brackets refer to the bibliography at the end of the paper. 
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case of a straight beam in bending by end couples (see fig. 2) might mislead some reader to the idea 
that the displacement components of the ring under its loads, shown in fig. 1, in its original plane 
must be of a higher order of smallness than those in the direction normal to that plane. This 
would mean that in the case of sufficiently small deformations, to which, by the way, we confine 
ourselves in the following, the bending of the ring in its original plane is negligibly small as com- 
pared with the bending out of that plane. Such a conclusion would be, in general, erroneous and 
the analogy with the case of the beam (fig. 2) is faulty: in consequence of the flexibility of the 
ring in the radial direction as well the possibility of appearance, under certain circumstances, of 
radial displacements cannot be discarded. This will become evident in later parts of the present 
paper. 

After this preparatory remark we return to our immediate objective of deriving the loading 
action diagram of the problem from our fig. 1. 


Separating from the ring (fig. 1) the part N P, we have to apply to each of its ends N and P, 


in addition to the twisting moment 5 F, not only a bending moment B, but also, in accordance 


with the preceding remark, a bending moment B, (fig. 3). 
Using the right-hand screw rule in the presentation of mo- 
ments by means of vectors, we easily verify the orientation 


Fig. 3. Equilibrium Sketch for Ring Part NP . Fig. 4. Equilibrium Sketch for Ring Part NS. 


se ae 1 
of the vectors B as shown in fig. 3. Considering namely the resultant of T and B at each of the 


points IV and P of the ring part N P we first see that these two resultants must be vectors of equal 
absolute amount and opposite direction; this leads to an inward vector B at the point N and an 
outward vector B at the point P. Now we secondly conclude that both resultants just mentioned 
must be parallel to the center line OC of the ring part N P (see fig. 3) by reasons of symmetry. 
This leads to the conclusion 


Boe=—— 1 cot F. (1) 


bol 


Ae 1 Mes 
The bending moments B and B, and the twisting moments i T, all shown in fig. 3, represent the 


system of loads, which determines the deformation and the stress distribution in the entire ring. 
This system of loads can be subdivided into two parts, the one consisting of the moments By and 


the other of the moments - T and B. The first part will be taken care of in later Sections of this 


+ ae 1 ; : ; , 
paper; the second part, consisting of the moments of and B, is the subject of the immediately 
following considerations. 


Separating from the part N P (fig. 3) a smaller part NS, corresponding to a central angle @ 
(fig. 4), we state for the latter again the equilibrium conditions. They require to apply at the 
point S a twisting moment T, and bending moments Bg und By. Since the resultant of Bg and 


1 é F A : 
To must equal that of the moments B and 5 T, appearing at the point P, we easily find from the 
vectordiagram at the point S the relations 
1 cos (O — F) ‘ 
B a, —_ —__— » 
Tease 5; near Te) saan, 2) 


1 en (OF). 


2 sin F 


(3) 


These two formulas give us the distribution of the bending and the twisting moments, produced 


1 é 
by, and representing the action of, the moments B and 71 T, applied at each of the ends N and P 
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of the ring part N P. This distribution of the moments (2) and (3) along the circumference of the 
center hee (see below, Section 5) of the ring is visualized graphically in fig. 5; the constant bend- 
ing moment By of thus far unknown magnitude appears in fig. 6. The two figs. 5 and 6 represent 
the loading action diagram for the dlacue system under poceileratae From this diagram we 
will have to derive all abe information necessary for the calculation of the stress distribution in 
the ring. 

Pecine for the time being, as just stated, the bending moments B, from consideration, we 


1 
can say that, as far as the remaining loading by the moments B and 7 T is concerned, both the 
external and internal mechanical actions exerted on the ring are now fully described by fig. 5. In 


UV= apparent 
discontinuity F. 
of To (see Sect 4) |~ 


U @ twisting moments 


re - 


(0) bending moments Ba 


FT cosec le 


Fig. 5. Illustrating Formulas (2) and (3). 


by 
Fig. 6. Diagram of Bending Moment B,. 


analogous conditions conventional systems of the one-dimensional class would be easily accessible 
to a “restricted” (By =0) investigation of the stress distribution, but not so in our present problem, 
see Section 1. In a certain sense of the word, the real problem now only begins. 


3. Examination of the Nature of the Elastic Problem of the Ring. Suppose the number n of 
the twisting moments is large. Then F =z/n is a small angle, and sin F ~» F, so that, according 
to (2) and (3), 


(a) twisting moments 


5 Tae: 

mill Bo & ee es =—= CONSts (4) 

bending moments (0) a 1 

t 0ST, 2 5T, (5) 

Tn 

ert 4 

06 (6) 

ee licane bee wee Somtul ties The diagram, figs. 5a and 5b, assumes in this case 


the form shown in figs. 7a and 7b. These results 
shall be compared with those of p. 138 of [1]. 

The uniformly distributed twist, denoted there by M,, is in our notation Tn/27 R, where 
R=radius of the ring. Thus equ. (a), p. 138 of [1], is, with M replaced by Bo, 


Tn 
Buy =M,R= 7_ (7) 


Fig. 7. Illustrating Formulas (4) and (5). 


in agreement with our equ. (4). So we see that the presentation (2), (3) remains valid even for the 
limiting case of n—>oo treated in [1], as far as the calculation of the bending moments 
Be from the applied twists is concerned. The torques — if any at all are produced by 
the concentrated twists — must obviously approach the limit zero with n—>oo, as the comparison 
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with [1] shows. The comparison of the latter elementary case with the general form of the problem 
will be continued below. 
With Increasing n the general case, illustrated by fig. 1, evidently approaches continuously 
. Se case just referred to, provided, of course, that the product nT does not exceed certain 
et eee : 5 ee 
# a ee this Ss ee and interesting to have at least some idea about the rapidity 
ormation of the general to the special case in question. With reference to fig. 5 such 
an idea can be derived from a simple and useful formula for the ratio o of the minimum value of 
the bending moment Bg to its maximum value. According to (2) we have 


ine cos (0 — F)6_0) - 1 
0 = oe (Gear, =—"cos f — cos = (8) 


re ratio is represented graphically in the chart fig. 8, showing @ as a function of n. The rapidity 
: the eute 0 with increasing n is remarkable. From the diagram fig. 8 we see that starting, say, 
rom n=8, the bending moment Be of the ring is practically almost constant and equal to its 


7 2 3 Y § 6 7 8 g 70 11 12 13 14 15 
Fig. 8. 9 = f(n). 


maximum value. This suggests the idea of the probability that starting approximately from 
n=8 the condition of the ring is already very near to that of a ring under uniformly distributed 
twist, characterized by the absence of torques. We shall see below that this absence of torques 
extends over a much wider range of the values of the number n. 

A further perhaps still more emphatic illu- 
stration of the somewhat surprising negligi- 
bility of the torques acting in the ring is 
obtained from the following preliminary con- 
sideration. It is obvious that the function 


D = 00) (9) 


where @ denotes the twist angle of the ring, 
must be of the general type shown in fig.9, in 
which A,, Ag, ... are the points of applieation of the twists T,, T,,.... From this diagram we con- 
clude immediately that the torques are zero exactly at the points of application of the 


oy aoe eet ere 


0(0) 


Qy ax 


A, 2 


Fig. 9. Illustrating Type of Function (9), 


twists (i. e. at the points, where the twists acting in the ring reach their maximum amount -5- 7) ; 


because the torques are proportional to the derivative dD/dO (see following Section 4), and this 
derivative is zero at the points of maximum twist. 

On the basis of the preceding introductory considerations of this Section 3 we may state the 
characteristics of the intended analysis as follows: 

The loading action diagram of our elastic problem is given by figs. 5 and 6 for Tg, Bo and By; 
with exception of Bj, which is an unknown constant, both the bending moment Bg and the 
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twisting moment J'g are known functions of the azimuth angle 9; the torque, however, which 
produces and determines the torsional stresses in the ring, is an unknown function of 0. The 
determination of this function is the important step implied in the concluding sentence of Section 2. 
To achieve this step we must know the function (9), the twist angle P(O), shown qualitatively 
in fig. 9. In an attempt of determining this function we would have to start from an expression 
of the form 
. oe: eee “th 
OG) = » 4; oe 2 (0S@Os8s FP) (10) 
j=0 
where the a; (j=0;1;...;.N) are indeterminate constant coefficients. An approximate presen- 
tation of the function ®(@) can be obtained by means of the formula 


D(O) = a, + a, cosn O (x = F) 020s Ff), (11) 


Since formula (11) fulfils already the boundary conditions d®/dQ=0 at O=0 and at O=F, 
there is good reason to believe that it must lead to results of practically entirely satisfactory ac- 
curacy. 

Thus our problem reduces essentially and mainly to the determination of the constants do 
and a, in (11). To determine them, as well as the constant bending moment By, we will have to 
proceed in a manner based upon Grammel’s theory of the ring under twist (see Section 1). Before 
reproducing here Grammel’s theory (Section 5) we shall insert in the following Section 4 a few 
remarks for the purpose of removing possible sources of misunderstandings. 


4. Remarks. It should be noted that, speaking of the substitution of a uniformly distributed 
twist T, for a system of a large number n of concentrated equally spaced and equal twists, the 
definition of T) must be given by the formula 


IPI 


5 12 
—~ O9nR (12) 


Ty 
(see second paragraph of Section 3), and not by means of the diagrams of the twisting moments 
in the figs. 5 and 7. This follows from the fact that the twisting moments shown in the latter two 
diagrams are produced moments, while both the concentrated and the uniformly distributed 
twists T’ and 7), respectively, are supposed to be externally applied moments, of course. 

Another source of possible confusion arises in the sign determination for the twisting moments 
(see figs. 4 and 5a). It is easily seen that in the vector diagram of the bending and twisting mo- 
ments, constructed as shown in fig. 4, the twisting moment Tg undergoes a change of sign at the 
point O=F, where Tg=0 (see fig. 5a); it leads to discontinuities, visualized in fig. 5 by the 
broken lines, of the function T,(@) at the points of application of the twists T. In other words, 
our full line diagram fig. 5a could be considered as being in some contradiction with the actual 
distribution of the (internal) twisting moments, a misunderstanding which is, however, easily and 
quickly removed as follows. The absence of sign changes in fig. 5a expresses only that along the 
circumference of the center line of the ring the moments T tend to twist the ring everywhere 
inward, This statement in connection with the indication of the magnitude of the twisting mo- 
ments T’g, given in fig. 5a, defines the twists completely and removes any possibility of con- 
fusion. 

The exact determination of the twist deformations is more complicated in the problem of the 
ring under twist by concentrated moments than in the case of a straight beam in torsion, see [3], 
pp. 222—224, and [4], form. (32) of p. 101. The first of the two sources just referred to indicates, 
on the basis of a study by K. Marguerre [6], that the necessary modification of the twist defor- 
mation formula is of the character of a correction; the magnitude of the error incurred by neglect 
of that correction is shown in fig. 16 as a function of the number n of twists. Since we are concerned 
here primarily with a particular fundamental aspect of the problem of the ring under twist by 
concentrated moments, namely with the question of the general way of computing the torques 
appearing in such a ring, we shall omit here any discussion of the modification just mentioned, 
using in the following the approximation by means of the corresponding straight beam formula. 

Finally, we once more call attention to the important fact that the moments, shown in fig. 5a, 
represent the twists, acting on the ring, and not the torques, which we must know in order to 
determine the torsional stresses in the ring. ‘To determine the torques, we must find the twist 
angle @ as a function of the azimuth angle 0. This is the characteristic difficulty of the problem. 
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5. Grammel’s Theory of the Ring Under Twist. Grammel’s study of the subject offers not only 
an instructive piece of analytical work, it presents a solid foundation for an effective quantitative 
discussion of the problem. More specifically, it provides us with all basic relations necessary for 
determination of the twist angle ®(@) along the circumference of the center line of the ring. This 
determination will prove to be, in a certain sense of the word, the decisive step in the process of 
applying Grammel’s investigation to the solution of the ring problem. A detailed acquaintance 
with this investigation is therefore essential for a thorough understanding of the solution of the 
problem in all of its practically important aspects. 

Referring to our fig. 10, virtually reproduced from p. 434 of [2], we note the following parti- 
cular features of the treatment in the two sources referred to in the concluding paragraph of Sec- 
tion 3: 

1. The position of the “neutral point” (called ‘“‘center of ro- 
tation” in [1]) is the subject of a careful study. A detailed 
derivation shows that the neutral point nearly coincides with 
the centroid of the cross sectional area of the ring, if the cross 
sectional dimensions are not too large as compared with the 
diameter of the ring. The result just mentioned is developed 
explicitly for the case of the uniformly distributed twist, and 
the conclusion is then used for the case of n equally spaced 
concentrated twists of equal ab- 
solute magnitude as well. In 
going over that derivation care- 
fully, we actually see that this 
question of the uniformity of 
twist distribution does not enter 
into it. 

2. The cross sectional princi- 
pal axes §, 7 of inertia through 
the centroid 0 of the cross sec- 
tional area enclose angles of 
amount ft with the outward ra- @ 
dial direction x and the upward 
axial direction y, respectively 
(see fig. 10). In many cases the 
angle t may be zero, but Gram- 
mel’s analysis applies more gene- 
rally to the case t~ 0. Vz 

Under the action of the loa- = _| — 
ding conditions, indicated in the Fig. 10. Orientation of Angles and Fig. 11. Illustrating Definition of 
diagram figs. 5 and 6, the cen- Axes According to [2]. Curvature Vector ¥. 
troid line of the ring will be- 
come a curve of double curvature. The curvature properties and characteristics of this curve 
must be given some special consideration. 

Using the terminology of differential geometry, we know that at any point of a space curve 
there is an “osculating plane” to such a curve. A “line element”, or infinitesimal arc of the curve 
at that arbitrarily chosen point, and the radius r of curvature of the line element, both lie in the 
osculating plane. 

Now let us consider the following elementary accessory problem. 


The curvature of a line element at a point C (C,, C,) of a space curve (see fig. 11), measured 
in the osculating plane (4,0, C,) of the curve at that point, is 1/r, where r =radius of the curvature 
circle A,0,C, of the curve at that line element. It is required to determine the curvature I/rp 
of the projection of the line element just mentioned on a plane (/’) enclosing an angle I’ with the 
osculating plane. 

We have represented the curvature circle AOC in fig. 11 by its projections A, 0, C, (ellipse) 
and A,0,C, (straight line) in plan and elevation, respectively, with the plane through C and the 
curvature center 0, normal to the osculating plane [and thus to the plane (/’) as well, see fig. 11], 
taken as the vertical projection plane, while the (/’) plane becomes the horizontal projection 
plane. 
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Since the curvature circle of point C(C,,C,) assumes the shape of an ellipse in the plan view, 
the curvature radius rp of the line element at the point C, of the horizontal projection of our line 
element is, with 0, C,=c, obviously 


sempre ee teks [3a 
LU Aa op c cos I” ( ) 
so that 
deeeh Le arm (13b) 
Taye Me 


The last relation shows that the curvature 1/r of a curve in general can be treated as a vector 0 

(Grammel in [2], [3]) defined in the following manner: 
The absolute magnitude of ) is the reciprocal of the radius r of curvature, while its direction 
is taken to be perpendicular to the osculating plane and combining, in line elements of outward 
convexity in the osculating plane, with the direction of the increase of the azimuth angle 0 
to a right hand helical system. 

This definition of the curvature vector ) is visualized in the elevation of fig. LI. 


osculating 
plane 


0; 


Fig. 12. Cross Sectional Plane Before and After Deformation. Fig. 13, Equilibrium Sketch. 


In order to avoid possible confusion and misunderstandings we call attention to the fact that 
the direction of measuring the azimuth angle @ is different in [2] and [3], which leads to differences 
in sign in derived formulas; this must be kept in mind in comparing the two sources just referred 
to with each other. 

Now we apply Grammel’s vectorial curvature definition to the space curve, into which the 
center line of the ring (fig. 10) is transformed under the action of the twists. More specifically, 
we are interested in the curvature properties of each of the two plane curves representing the pro- 
jections of the space curve under consideration on the planes §O¢ and 7O€¢ (fig. 10). 

According to our preceding detailed presentation the vector f of our space curve determines, 
for each line element, the curvature of each of its two projections just referred to: with the con- 
vexity (or concavity) of the line element remaining unchanged in the projection (see definition 
of »), the curvature of the projection on the 70¢ plane is the component k, of the vector f in the 
& direction, the curvature of the projection on the §O¢ plane is the component k, of the vector £ 
in the 7 direction. 

Let us apply these statements to our ring of fig. 10, reproduced in fig. 12. The curvature of 
the originally plane circle is represented by the constant vector f, of the absolute amount 1/R, 
while its direction is that of the y axis. The components of this original vector f) in the € and 7 
directions, respectively, are, according to fig. 12 


Bay 
ko: = Rp sine, (14) 


kon = 2 cost. (15) 


Considering the same originally plane circular center line of the ring after deformation, we find 
its curvature to be represented by a vector f, which is in general, of course, a function of the azimuth 


XXVI. Band 1958 I. Malkin: The Elastic Problem of a Homogeneous Circular Ring Acted Upon 205 


angle © both with respect to absolute amount and to position in space. From the definition of 
the curvature vector (see above) it is easily seen that the principal axes OF and Ov and the vector f 
are all situated in one plane, namely in the cross sectional plane of the deformed ring at the line 
element considered (see fig. 12). The components of £ in the principal directions O& and Oy shall 
be denoted by k, = k;(@) and k, =k,(@), respectively, so that the curvature changes of the 
plane projections of the line element of the space curve on the planes 7 0¢ and £0€ are, with re- 
ference to (14) and (15), respectively, 


Ak, = k{O) —— sint, (16) 


Ak, = k,(@) —= cost. (17) 


It is necessary to keep in mind that the curvature changes (16), (17) must be defined relatively 
to the principal planes §O¢ and 70€¢ disregarding the movement of these planes in the space, 
because the bending problem must be treated with respect to the principal axes of the cross section. 

Introducing the modulus of elasticity E and, with reference to fig. 10, the moments of inertia 


Df aA, (18) 
I= {dA (19) 
where dA = element of the cross sectional area of the ring, we define the flexural rigidity constants 
ELCs 5 (20) 
EI,=(C, (21) 


for the bending action in the 70¢ and §O¢ planes, respectively. 

The next step is to establish the relations which exist between the curvature changes (16), 
(17) and the moments acting on the elastic system under consideration. 

Refer to fig. 13, representing again our ring with a part O,Q cut out of it. This part 0,Q is 


i : 
acted upon by the moments of, af ey F and By at the point O, and by the moments Tg, Bo 
and By at the point Q; the resultant of oak and af cot Fis a vector equal and opposite to the 


' vector, representing the resultant of Tg and Bg at the point Q; both vectors are parallel to OP, 
the angle 0,OP being equal to F. Draw now a plane through OQ normal to the osculating plane 
of Q; looking at that plane in the direction of the arrow « we find the part QO, to be acted upon, 
at the point Q, by the moments Bg and B, as shown in the lower part of fig. 12. Note that our By 
is identical with Grammel’s — Nt, . From fig. 12, lower part, we find the total bending moment at 


the point Q to be 
M, = By sin (t + ®) — Bg cos (t + ®) in the 7 O¢ plane, (22) 
M,, = B, cos (t + ©) + Bo sin (t + ®) in the OC plane. (23) 


These bending moments are responsible for the curvature changes (16) and (17), respectively, 
Therefore, by virtue of (20) and (21) and with the usual sign conventions for the bending moments. 


M, =— C, Hk, , (24) 
M,, = — C, Mk, . (25) 
Substitution of (22), (23), (16), (17), (20), (21), into (24), (25) gives 
— By sin (t + ®) + Bog cos (t + ©) = EI, [k; — (1/R) sin t], (26) 
— By cos (t + ©) — Bo sin (t + ©) = E I, [k, — (1/R) cos t] . (27) 


Before proceeding further with the derivation of the necessary formulas, we reproduce here Grammel’s 
simple and instructive check of the relations (26) and (27). Grammel calls attention to the fact 
that the latter two formulas take care of the two kinds of bending which take place in the ring, 
namely of the bending proper (bending caused by the actual shape of the deformed center line 
assumed to be produced by bending only without twist) plus bending produced by twist (like in 
the elementary case treated in [1]); these two kinds of bending can be termed ,,flexural bending“ 
and “twist bending’’, respectively. ; 

If the actual shape of the center line of the ring after deformation were free of twist, then the 
actual total angle t + © would have to replace the original angle t, as far as the computation of 

14 
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the curvature changes (16) and (17) is concerned. The corresponding bending stresses (“flexural 

bending”’) are, as easily seen, ’ 
0, = — £7 {hp RS em (tJ) |, (28) 
o; = EE [k, — R-“cos (C= ®)]- (29) 


To explain the signs in (28), (29) we refer to the lower part of fig. 12: if the curvature change (16) 
is positive, then filaments, characterized by positive values of 7, undergo contraction; if the cur- 
vature change (17) is positive, then filaments, characterized by positive values of €, undergo ex- 
tension. This takes care of the “flexural bending” stresses in the ring. 

To obtain the expression for the “twist bending” stresses, we return to fig. 10, a part of which 
is reproduced in fig. 14. Referring to the notations of the latter, we see that a rotation of a cross 
section of the ring through the twist angle ®(@) produces a local extension of the radius of the 
filament through the point P by the amount 

r [cos (U + ©) — cos U] = — x (1 — cos Bb) —y sin @. (A) 
In the case of a ring, whose cross sectional dimensions are small as com- 
pared with the radius R of the ring, the “twist bending” stress is on the 
basis of (A) given by 
p E 


oy) = — BR [x (1 — cos d) +ysin O]. (30) 


a From fig. 10 we have the relations 
y x =é&cost—y7sint, (31) 


Fig. 14, Illustrating Derivation 


of (30). y =€ sin t + 7 cost (32) 


which are to be substituted into (30). The sum o of the three expressions (28), (29), (30) is a function 
of € and 7. Consider now the integral + fo 1 dA, taken over the entire area A of the cross section 
of the ring. Since {7 dA, taken over the entire area A, is zero, only the items (28) and (30) are 
to be taken into account, so that 


i fon dA =—E [k, — RD ain (¢ + )| | n* ee = | n {(E-cost —n sina) (1 — cos D) & 


+ (Esint + 7 cost) sin @}dd =— ET, E en (33) 


Entirely analogously we find 
— [asda == — E [k, — R™ cos (¢-+- 0)] [era -- R | EE cot sin t) (1 — cos ®) + 


+ (€ sin t + 7 cos t) sin @} dA = — E ups ‘, = = cos t| . (34) 


In other words, the integrals representing the left hand members of (33) and (34), which are the 
integrals representing the bending moments in the 70¢ and in the €O¢ planes, respectively, are 
correspondingly equal to the right hand members of (26) and (27) with sign reversed, which is 
correct in view of the relations (24), (25). This completes Grammel’s check. 

The curvature of the projection of the deformed center line of the ring on the coaxial circular 
cylindrical surface through the undeformed center line is 


k, = ke cos (t + ®) —k, sin (t + ®) (35) 
and similarly the curvature of the projection of that deformed center line on the plane perpen- 
dicular to the axis of the cylinder just indicated is 

ky = kz sin (t + ®) + k, cos (t + ®). (36) 
Multiplying (26) by C, cos (¢ + ©) and (27) by C, sin (t + ®) and subtracting then the second 
result from the first we get equation (37). Multiplying (26) by C, sin (¢ + ®) and (27) by C,cos(t + ®) 
and adding the results, we obtain equation (38): 


; 1 
C, C, (k, + R sin OD) = = (Gi — C,) By sin 2 (¢ + ®) + Bo [C, sin? (¢ + D) + C, cos? (t + O)], 
(37) 
1 1 ; 
C, C, (kz, — R cos 6) = — 5 (C, — C,) Bo sin 2 (t + &) — By [C, sin? (t + ®) + C, cos? («+ B)]. 


(38) 
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Considering the projection of the deformed center line of the ring on the cylindrical surface, 
we denote by u the deflection of the projection positive in the + y direction. Progressing along 


the part QO, in fig. 13 and remembering that by definition k, is positive outward, we have the 
relation 


du ; 
a rR ky : (39) 


Denoting furthermore by v the radial contraction of the projection on the original plane of the 
center line of the ring, we have the differential equation (see [7]) 


d’v ee 
coat e = R(ky—R4). (40) 


Thus far we have followed closely Grammel’s work. The set of equations reproduced in this 
Section from [2] is the foundation for the concluding steps, which will be carried out below. 


6. Remark on the Integration of the Differential System of the Problem. The further treatment 
of our ring problem must be based upon the choice of an appropriate hypothetical expression for 
the twist angle @ as a function of the azimuth @. This has been emphatically anticipated already 
in Section 3. We shall use in the following the hypothetical formula (11) with a, and a, as indeter- 
minate constants. This formula (11) is to be combined with the equations (37), (38), (39), (40) 
representing a system of four equations for the four unknown functions 


key» Ka, u,v (41) 
of the azimuth 0. The integration of this system of equations must lead to a solution appearing 
in the form of a system of functions of the argument © and of the indeterminate constants Bo, 
dy, a, and the constants of integration. The constants will have to be determined from the boundary 
conditions in connection with the condition of minimum strain energy [8]. The entire work thus 


sketched will be carried out in the following for the case of small twist angles ®, for which @ can 
be substituted for sin ® and | for cos @. 


7. The Solution. Returning to the equations (37), (38), we treat the case 
GezCG, Osmall t>06 so that t--Oat. (42) 


This case shall be discussed here in a manner very similar to that used in [2], pp. 436—437. Thus 
we introduce the constants 


= Green euzie KR. Groom’) Cycoser—L, Cysint + C, costt = NV. (43) 


With these notations the equations (37), (38) can be approximated by 


GG(h + R)=K Bo +E Bo. — ©, 6,(hy— ge) = K Bo +N By (44) 

Rodithis iaby yaiue of (2) and (11) 
CC, i a * (Guo 08 p ) ee (45) 
a7 (i =z) eT (46) 

Substitution of (39) and (40) gives 
pe Ce K peer te ==R (a) + a, cos nO), (47) 
ee ei oo Wy Bea: cle at | (48) 


where a prime denotes differentiation with respect to @. The general integrals of (47) and (48) are 


1/R?K B, ps ne iy ee se OL ar, 6) eal 4O 
u— CG 2,0 xf Cas Ra) @ Come Jie n? ( ) 
R? N B, KTR 


(O— F) + cos (0 — F)—cos F cos QO), 
(50) 


v =c, cos O + c,sin@ CG 1C,C,snF [O sin 
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where c, to cy are constants of integration. For reasons of symmetry we must have 


du/dO = 0 dv/dO = 0 when 0 = 0 and when O=F. (51) 
Therefore, according to (49) and (50) 
Cpe ROE Ti 2a, (52) 
'(C,-CaIR) op Bees Pet 20h (53) 
c, =—(K T R28 G, G, sin? F) (2 F + sin 2 F) (54) 
Cpe) (55) 


while the constant c, is determined by the choice of the reference line for the deflections wu. 
The strain energy of the ring part between O = 0 and O = F is given by the expression 
F 


R [((Mi | Mz aa 
Q {| G = Cy F G, a6 eo) 
0 


where M, and M, are defined by (22) and (23), respectively, with t + ® replaced by t, because 
of (42), while M,, representing the torque in the ring, is determined by (11) as follows. Denoting 
by G the modulus of rigidity of the ring material and by I, the polar moment of inertia of the 
cross section of the ring, and putting 


G I, = C, (torsional rigidity) (57) 


we have the formula 


M, = C, (d®/R dO) = — (n G,/R) a, sinn@. (58) 


The strain energy (56) must undergo the process of minimization with respect to the constants dy 
and a, occurring in (11). The indeterminate constants appearing in (56) are a,, occurring in M,, 


_ equ. (58), and Bo, occurring in M, and M,, equ. (22) and (23). By virtue of (53) we have, however, 
d/da, = (C, C,/K R) (A/0B)) . (59) 
Thus the process of minimization produces the two equations 


F F 
f Mt (0M,/0a,)\dO = (n C,/R)? a, f sin? n O dO = 0 (60) 
0 0 


leading to a, = 0 (61) 
and 


F 
[1671 M, (@M,/9B)) + C1 M, (@M,/2B,)] dO 


0 
= {[C,? (By sin t — Bg cos t) sint + Cy! (By cost + Bosint) cost]|dO =0. (62) 
0 
By virtue of (2) the last equation becomes 
F B, (CD sin?t + C,' cos*t) = (C,*— C4) i T sin t cos t (63) 


from which finally 


(Cm —=C, |) VT sint cost (64) 


; 2 F(Cy" sin? t + Cy" cos? t) 


Substituting (64) into (53) we find ay, so that all of the three indeterminate constants become 
known; this determines the displacement components (49) and (50), the curvature components k, 
and k,, defined by (37) and (38), respectively, and the moments M;, M,, M,, defined by (22), 
(23), (58), respectively. 


We conclude the above considerations by a brief discussion of the results just obtained. 


8. Discussion of the Results. Perhaps the most important result is expressed by equ. (61): 
under the assumptions, which represent the basis for our considerations, namely the assumptions 
of small initial curvature and small deformations, especially small twist angles, and with the 
reservation expressed in the semifinal paragraph of Section 4, the twist angle @ is, according 
to (11) and (61), a constant, namely a, independently of the number n of twists; the torque in the 
ring is, therefore, zero according to equ. (58). 
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Turning to equ. (64) we substitute F = z/n into it; if we then assume T'n to remain constant, 
then By remains constant as well, provided that C,, C, and t > @ do not change either. Further- 
more we note, with reference to the third paragraph of Section 2, that the condition for the dis- 
appearance of the bending moment B, is C, = C, = G; this, however, automatically involves an 
angle t of arbitrary magnitude (since any pair of mutually orthogonal axes can be used then as 
the pair of the principal axes; see also below), so that in particular we may assume t = 0 in the 
case just stated. This case of C, = C, = C will be considered here somewhat more in detail. 

If C, = C, = C, then, according to (43) and (53) to (55), 

RT RT 


Koa 0 i Eee eG fea oral PS eG ea — 1S ey ree (65) 


Substitution into (49) and (50) gives then, with B, = 0 and Cy = (2 C/R? T)c,, 


DEG i. N op I ale I 
ep cto glee pr at ral oth Jet a (66) 
v=0 . 
Equ. (45) and (46) give now, on the basis of the first, the second and the fourth of the equ. (65), 
Ff Et 1 
h=F¢|iar°SO—F)—1; =z, CF 


It is easily seen that, as soon as F' becomes small of the first order, k, becomes small of the second 


order, as it should be. 

Since ® = const = ay, where ay is given by the fourth of the equations (65), independently 
of the number n of twists, the twist angle ® must be the same as in the elementary case treated 
in [1], a conclusion which affords a valuable check of our calculations. Substituting # = z/n into 
the formula for a) and using our formula (12), we find 


2 
© =a, RT PN (68) 


2Cx%  C 
in complete agreement with form. (126), p. 139 of [1]. The remaining question is this: while our 
formula (68) is derived on the basis of the assumption that C, = C, = C, such an assumption 
is not used in [1], a fact which might suggest that the identity of the two cases, just compared with 
each other, is somewhat questionable. The situation will be clarified presently. 
Let us return to the analytical formulation of our problem by means of the system (11) and (37) 
to (40). We consider two special cases, namely 
(Il) C, =C,=C, t arbitrary, so that we may assume t= 0, 
(LD) Coe2 iG, f= 6 
while ® is assumed to be small in each of the two cases, like in all of our preceding considerations. 
While the special case (I) is covered by the solution given in Section 7, the special case (IT) 


has not been treated yet. | 
In the special case (II) the system (11), (37)—(40) obviously transforms into 


@D 
@® =a,+4,cosnO, Cy G + R) = Be. (69) 
1 
C, (Hs —3 ss By ? (70) 
we Re (71) 
vo’ ty = R? (Hs —z) ‘ (72) 


In the special case (I) we obtain the same system of equations with the only difference that now 
eC ge C, 


From (70) and (72) we get, with c; = constants of integration, 


R2 
poe ope v = ¢, sin @ + cy cos O — ~ Be (73) 


2 


v' = c, cos 9 —c,sinO. (74) 
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The symmetry conditions v' = 0 when 9 = 0 and when @ = F lead here to 


Cs — Cg a 0 (75) 
and 
R2 
( — yim By . (76) 


Note that the result applies to the special case (II), in which C, A C,, t = 0; this special case 
includes, however, the possibility C, = C, = C, t = 0, as well, which is the special case (I), and 
we have incorporated this latter special case into the general case treated in Section 7, see above. 


Therefore (75), (76) lead to the conclusions [cf. equ. (65) to (67)] 
Bo e=0, =O), | all Wee (77) 


The equations (69) and (71) are to be treated in an entirely analogous manner. We give here only 
the results, leaving a detailed derivation to the interested reader: 


R°T cs(O—F) RT -.., Ry eRGE P a 
iow ED eae Sey 7CF? Pega 00 TO cs 2C, Owe? (ee 
RT 
ee OM oe & i 
1 


Comparing now the fourth of the equ. (65) and equ. (79) with form. (126), p. 139 of [1], referred 
to in the two paragraphs following our equ. (67) and (68), we actually see that the identity indicated 
there is genuine; the elementary case treated in [1] applies to the case t = 0, and its’form. (126) 
is actually identical with our form. (79). This completes the check indicated above. 


Me Error in °/o 


1 7 7 
od gt Zu 


Fig. 15. u = f(@) for n = 2. Fig. 16. Illustrating Error in Twist Deformation Formula. 


Since our equ. (58), (60), (61) remain valid in the case considered here, the torque in the ring 
is again zero in the special case (II) as well, and furthermore equ. (78) becomes essentially identical 
with the first of our equ. (66), with the only difference caused by the constants C, and C. This 
takes care completely of the special case (II) as well. 


As an illustration we give in fig. 15 a diagram of the function (66), for the case of n = 2, for 
which that function assumes the form 


2C : @? 
er uso 0 — sinO ote (80) 


9. Conclusions. With the reservation and assumptions stated at the beginning of the preceding 
Section 8 the results derived in the present paper can be briefly formulated as follows: 


A homogeneous circular ring acted upon by equally spaced concentrated twists of equal magni- 
tude is always free of torques. The bending moment By, represented by a vector normal to the 
original plane of the ring, vanishes, if one of the two cross sectional principal axes is situated in the 
plane just mentioned; the case that the two principal moments of inertia of the cross sectional 
area of the ring are equal to each other is a special case of the one just taken care of. If the angle t 
is essentially different from zero [cf. definition of case (I) in Section 8], which means that C, 
and C, are actually different from each other, then the bending moment B, is different from zero 
and given by our formula (64). This is basically the solution of our problem of a circular ring 
under concentrated twists, because the effects of the bending moments Bg, see formula (2), can 
be studied by means of elementary methods. 
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Torsionsschwingungen von Scheiben mit Exponentialprofil 


Von S. Tameroglu 


1. Einleitung. Die Gleichungen der Torsions- und Dehnungsschwingungen ohne Knotendurch- 
messer von Scheiben mit beliebigem Profil lauten bekanntlich wie folgt* 


ao 4 (s 1 ¥f alt y ae Dye: 00 4 efit 
C 


re dr y r} Or eG \ or ot r (1) 
Go , (yt Ucona ae : untae) (One 5 or es 0; 
Or?! (a: ¥ aN r |ér dr mry 7 gemE \ av ot : 


dabei sind r und @ die Polarkoordinaten der Scheibe (r;<r<r,) vom Profil y(r); ferner bedeutet 
Or, t) den Zentriwinkel des Schwingungsanteils der azimutalen Verschiebung, o(r, t) den Schwin- 
gungsanteil der radialen Verschiebung, y das spezifische Gewicht der Scheibe, g die Erdbeschleu- 
nigung, G den Schubmodul, m die Querdehnungszahl, EF den Elastizitatsmodul und @ die Dreh- 
geschwindigkeit der Scheibe. 

Diese Gleichungen lassen insbesondere stehende Eigenschwingungen von der Frequenz « gemaB 


dem Ansatz 
P= (rein 27506, 0 = hr)cos 2ornt (2) 


zu, deren Amplituden Y und R Funktionen von r allein sind und nach (1) den Differentialgleichungen 


dy 5 i 1 ; es L 42 | ® yioeeg Bi } 
My Aare r 


dr dr y I Fa) he 


dR dy 1 ; 1\dR : dy {| Ue oF De ayy 
Gr | ea if | (1 i dr mry ay eens) gat oka 


(3) 


mit den Abkiirzungen 


J, == AGE Ae p= 2n0]/ 


Tae 
gmk 

gentigen miissen. Bei rotierenden Scheiben besteht fiir torsionale Stérungen in erster Linie mit 

den Torsionseigenschwingungen Resonanzgefahr, in geringerem Mafe aber auch mit den Dehnungs- 

schwingungen; jedoch darf man die fiir die Resonanzgefahr maBgebenden Eigenfrequenzen x bei 

beiden Schwingungen ausreichend genau an der nichtrotierenden Scheibe, d.h. fiir m = 0 er- 

mitteln®?. Damit kommt statt der Gleichungen (3) einfacher 


So 


dy dy 1 Py Kae 

dr? (a y #4 | dr Ee 0, o 
d?R dy 1 1\dR Se! 1 s 
dr? =r (ae y | dr r (« dr m ry =) 1 } 


Diese Gleichungen sind fiir @ = 0 nicht mehr gekoppelt und kénnen jede fiir sich gelést werden. 


In der vorliegenden Arbeit soll nun zunachst die Gleichung (4) fiir Exponentialprofile von der 
Gestalt 


(ya aC Cm (r/c)P (6) 


behandelt werden, wobei a, c und p reelle Zahlen bedeuten. 


2. Differentialgleichung. Aus (6) berechnen wir 
aye ae r\pl 
dr y Ie ( | r (7) 


are C. B. Biezeno u. R. Grammel, Technische Dynamik, Bd. 2, 2. Aufl. S. 33, Berlin/Géttingen/Heidelberg 
Dos 


2 Vel. C, B. Biezeno u. R. Grammel, a. a. O. S. 35. 
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und erhalten damit statt (4) die Differentialgleichung 


dy — ?P CG | dy se 
: E CES | ar Ties 
dr? 1 L r J dr Aa = 0 (8) 
Fihren wir durch 
( rei : p= 100) — 72 2 2.27% 
apa i es mt = 0) Ged PSP OS 47 oe (9) 
dimensionslose GréBen ein, so kommt statt (7) 
Sones 
7! aL sone 7 = u2 ” — 0 , (10) 
S 


wenn wir weiterhin die Ableitungen nach € durch Striche bezeichnen. 

Sofern wir fiir den Exponenten p in (6) nur natiirliche Zahlen p = 1, 2, 3, 4,... zulassen, kann 
man partikulare Lésungen der Differentialgleichung (10) auf folgende Weise finden. Wir nehmen 
an, daf} sich diese partikularen Lésungen Npi() nach Potenzen von & entwickeln lassen: 


7 pr(é) =I1+ A596 Ae a436° +4- oe -}- nec’ sees, (11) 


Dabei haben wir wegen der Linearitat von Gleichung (10) a,9 = 1 gesetzt; ferner ist a,; = 0, 
da nach (10) fiir eine regulare Lésung 7/'(0) verschwindet. Die Koeffizienten Ap2, Ap3,-.~+ bestimmen 
sich rekursiv aus (10). Um anzudeuten, da sie noch Funktionen von p und u? sein werden, schreiben 
wir statt (11) pragnanter 


Np1(S) = 1 + fpalu®, p) & + fys(u®, p) & + fra(u®, p) 4 ++. (12) 


Um weiter aus den partikularen Lésungen 7p) (€) die allgemeinen Liésungen7,(&) zu finden, 
gehen wir von der Darstellung 


Np = (5) Np (E) (13) 
aus. In Gleichung (10) eingesetzt, kommt 
/ vr 3—pé Il , 
wo Mp1 + 2 Op 1 8 Ypi ee —(v' Mp1 + Mp1) + We Nn = 9 (14) 
oder, da 7,,, der Gleichung (10) fiir sich schon geniigt, 
, Se Bh, 
vw Heg t (20 == ra ps) es (15) 
Dies ist fiir 
qo (16) 


eine Differentialgleichung erster Ordnung 


oi 3=pe 
a! tps +9 (2793 aaa tp) = 0 


mit dem allgemeinen Integral 


HeeéP 
/ 
q=v= (17) 
: Pup 
und der Integrationskonstanten H. 
Eine weitere Quadratur ergibt mit der zweiten Integrationskonstanten Q 
& 
, ey 
vee H | rae pes (18) 
Ss Np1 
und die allgemeine Lésung von (10) ist damit in der Form 
Hn, | reas (19) 
= = — ce 
"p QO Np. -| Npt Bae S 


a 


gefunden. 
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Um das Integral in (19) auszuwerten, entwickeln wir 


egP 


why bt ber ae Pree Gees ei (20) 

Damit kommt aus (19) mit der weiteren Konstanten U 

: 1 b &2 &3 
Np =9Np1 + Hpa( 9 &2 = Opa We bog Se Coase Ops te Bee v) (21) 
oder mit 
D=Q+HU (22) 
1 b 1 

tip = Dtpr + Hips ( rp Se F Opagln g's .07, 6 oy Seed (23) 


Wir haben somit das allgemeine Integral von (10) mit zwei willkiirlichen Integrationskonstanten 
D und H als Kombination des regularen partikularen Integrals 7, mit einem zweiten partikularen 
Integral dargestellt. Die hier auftretenden Entwicklungskoeffizienten b,,, 6,9, 6,3... sind durch 
(20) eindeutig bestimmt. Man ermittelt sie aber bequemer, indem man mit (23) in die Differential- 
gleichung (10) eingeht. Man sieht dann, da auch diese Koeffizienten noch von p und u? abhangen: 
daher schreiben wir statt (23) 


Mp =N(D+ Hops) (p= 1,253.45...) (24) 


mit der Abkiirzung 


1 Fyi(u, f) 
boalt®, p) = — gy — PHP) 4. Fr afu, p) Im E+ Foo(tt, p) E+ + Fpaluts p)et+--+-. 25) 
In dieser Form stellt sich die allgemeine Lésung der Differentialgleichung (10) fiir bestimmte Werte p 


und uw? dar. 


3. Randbedingungen. Die allgemeinen Randbedingungen fiir Torsionsschwingungen lauten fiir 
federnde EKinspannung am Innenrand r = r; und am Auf®enrand r = r, 


(Bom B, Wr) . (26) 
(nr = Be Er » age 


Hier sind die B; und B, jeweils gegebene oder leicht berechenbare Festwerte. Auf die dimensions- 
losen GréBen (9) umgerechnet, hat man statt (26) und (27) 


vp(6:) =cB,; 1 p(S;) > . (28) 
Np(Eq) = ¢ Bn, (E,) - 


Beziehen wir noch den einen Parameter ¢ des Scheibenprofils (6) und den Innenhalbmesser r, 
gemaB 


en ee. (29) 
auf den AuBenhalbmesser r,, so wird 
Vy r Xe 
to my ; 
fo ta ©. CY) 
a Ta 7 Hy ° 
und die Randbedingungen (28) lauten mit (24) 
, x a x “4 ae ut. H. os 
a am 
we 1 erie al 1 co 
io +H a(2n )~e Rn) +H 
oder geordnet 
A [(p1 Spo)’ — ¢ B; (“pr Sp2)]. au + D [(yp1) — ¢ B, (Npa)I, ies 0, 
S 2/1 ; S = Kq/ Hy 39 
H Uno Cpa)’ — © Balpr Spal), _, + Plt) —-¢ Beryl), _,, =9- aa 
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Damit dieses in H und D lineare, homogene Gleichungssystem eine nichttriviale Lésung hat, 


mu seine Determinante verschwinden. Wir haben somit die Frequenz- oder Kigenwertgleichung 
fiir u?, p 


[7 pt a =a €Bi()p1" pe) |. = we), [pt gt B; "pil: = xy/ry | 
19 tet ; eine s()s, 33 
[Mp1 Sp)’ —¢ Ba(nps Cpadle — ie Apt =e B, Haale ‘cel i) 


Wie diese Gleichung numerisch gelést wird, soll an den folgenden Beispielen im einzelnen gezeigt 
werden. 


4. Beispiel I. Wir wenden das angegebene Lésungsverfahren auf Scheiben an, die am Innen- 


rand starr eingespannt (B; =co) und am Aufenrand frei sind (B, = 0), so da® die Randbedin- 
gungen einfacher 


Nip(Si) = Np (Ea) = 0 (34) 
lauten. Fiir die Parameter des Scheibenprofils (6) wahlen wir mit x, = 2,5 
. p=, C= OT. (35) 
also 
y=ae P's, (36) 
Die Differentialgleichung (10) vereinfacht sich damit zu 
i SF a + uy =O) (37) 


und kénnte mit konfluenten hypergeometrischen Funktionen allgemein gelést werden. Dadurch 
wirde aber die Kigenwertsermittlung schlieBlich ebenfalls auf Reihenentwicklungen hinfihren. 
Wir bestimmen den niedersten Eigenwert © 


ie tt e Te ee 
De 35S DIU 0 eee (38) 
in Abhangigkeit vom Halbmesserverhaltnis x, = 1,/r, . 


Zuniachst erhalten wir durch Einsetzen der Entwicklung (11) in (37) fiir die Koeffizienten aj, , 
die Rekursionsformel 


[n (n — 1) + 327] a1,n—(n— 1) ),,-1 + wa,» = 0, (39) 
woraus fiir die Partikularlésung (11) im Falle p = 1 
u? u2 u?(2—5 u?) w(2—=9u*) u2 (8 — 50 u? -+ 35 ut) 6 
(Nes Sane eae ee Ao bay 35 240 - 28-48 . 


u? (40 — 314 u? + 463 u4) a ror (40) 
240 - 28-40-63 


kommt. Fiir die zweite partikulare Lésung kénnen wir das zweite Glied der Gleichung (23) ver- 
wenden. Geht man mit diesem Glied in die Differentialgleichung (37) ein, so ergibt sich 


(21 1) (a t a4 “8 + baa + buat + base? + big? + bia é! peg 


im ( a “at PEt bag + 2OgE + 3 bygd# + Abs +o) =0 (41) 


oder mit (40) 
3 Tu? us u2(6—55 u?) .. _u?(34—293 u*) 4 u?(72 — 674 u? + 875 ut) 
eo 8 g 40° 120-8 a3 ~ 40-24 - 35 240 - 28-80 
u2 (520 — 5426 u? + 12067 u4) 0 
240 - 28-72-70 


5 
é 


|e = Pas | bis | bi § | bys Bt bag + o-| 


u? u(2—5u*),, w(2—9 u*) £5 u2 (8 — 50 u? + 35 ut) 


u2 7 
3 seats 
+ [1 ge F — GF 40-24 ° 240 - 35 240 - 28 - 48 5 


fa? 2 3 uA 3 26 b 2 = 
i 7 ao 2 = | | $4 Bs ra t+ Bag + 2 by5F + 3 dy G? + 45,6 +--+) = 0. 
(42) 
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Dies fiihrt auf die Identitat 


1 u*\ 1 5 uu uw 1 
(b4, — 1) Bai (2 er ae ae [3 bys bis 8 by + i) E 


6 u2 2 4 Tu? u2 u? (45 u2 —2) 
r ‘ by — bys ee a ms a zi) a7 5 bis— bu — 9 bss 359 812 T 9p. 8 


u? u? (25 ura 2) u? (221 u? — 18) 


2(37 u2—2 | ?—18) , 
5 Na 7 a\eae 6 Bie Sabir ats agg ei S08 a eae 
2 
u?(8— 106 u? + 175 u*)I .. ow u? u> (6 — 55 u*) 
940 -28-20 Sonate ease 8 bis 24 bua 120-8 bis 
u? (26 — 257 u?) u? (136 — 1522 uw? + 2275 u*) u? (280 — 3542 u? + 9289 u4) ele. Som 
) 240535 4 © 240 - 28 - 240 1 240 - 28-72-70 
(43) 
aus der sich die Koeffizienten b,,, nacheinander zu 
2+ uw? " Dee line 
by = Fy, p) = 1, b,. = F,,(u*, p) = Ke 12 bs = F,,(u?, p) = Tey ee 
5 40 + 156 u? + 35 u4 & s a 280-+- 2116 u2 + 1573 ut 
big =a Pyy(u’, P) == : 960 9 bis sieaic F,5(u ’ P) 940 - 140 a) (44) 


F 1120 + 14984 u2 + 24734 ut + 3325 us 

big = Fya(u’, p) = 240 - 280 - 12 x : 

23520 + 528584 u? + 1622894 u*4 + 708701 u® 
7 240 - 280-21 - 84 


biz = F,,(u’, P) — j 


bestimmen. 


In (24) mit (25) eingesetzt erhalten wir als allgemeine Lésung der Differentialgleichung (37) 


l Ug Ue wr ow) we 8) £5 u?(8 — 50 uP + 35 u') £6 
ee fies 60° Mio BAS 240 - 35 240 - 28 - 48 
u? (40 — 314 u? + 463 u*) : Le eee ig u2 Pe ae cig NOEs Ze 
O10 20 r40mGs Shes (| paras 262" € 5 4 Ings 60 . 
_ 40+ 156u2 + 35 ut... 280 + 2116u? + 1573 ut, | 1120+ 14984u?+ 24734 ut+ 3325 u8 ,, 
r 940-8 ° 240 -140-3 o 240 - 280-48 $ 
_ 23520 + 528584 u? + 16228964 ut + 708701 uf ., (45) 
=, 240 - 280-21- 84-5 ese ‘ 
Die Frequenzgleichung (33) lautet wegen B,;= oo, B,=0, x, =2,5 jetzt 
£42(0.4 x5) 1 
10( 2 is 0 : (46) 


| (iu Cie) e = (hf 14(0,4) | 


Um sie aufzulésen, wahlen wir folgendes Naherungsverfahren. Wir beriicksichtigen von der 
Entwicklung (45) nur die ersten Glieder gemaB 


u2 oie 1 1 2 + wu? 10+ 17 u? 
so da wir also 
ur. ue 
Nn = 1 Se Gta ee 
(48) 

a 1 1 2+u 10+ 17 v2 

jaa 9 £2. = In € 3 

a 2 & E 4, 60 


haben. Damit wird (46) eine Gleichung von endlichem Grad in u2, die numerisch gelést werden 
kann. Den so erhaltenen Naherungswert fiir den niedersten Eigenwert u2 werden wir hernach 
verbessern, indem wir von der Entwicklung (45) weitere Glieder hinzunehmen. Der Gang der 
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Rechnung im ei i ir das 7 Wea —005 ; * 
g einzelnen sei nun fiir das Zahlenbeispiel %_=0,05, also 0,4 x%,=0,02 vorgefiihrt. 


Aus (48) kommt 


€12(0,02) = — 1301,952679 — 0,972339 u2 , 
C12(0,4) = — 6,016479 — 0,115740 u2 , | 
14(0,4) = 1 — 0,021067 u2, (49) 
C12(0,4) = 23,291667 + 0,908333 u2, | 
nia(0.4) = — 0,108000 u2 . J 
Mit diesen Werten lautet die Frequenzgleichung (46), wie folgt: 
— 1301,952679 — 0,972339 u2 | 
23,291667 + 1,067427 u? — 0,006636 w! — 0.108000 u2} Sy 
oder ausgerechnet 
0,111649 ut + 139,543465 u? — 23,291667 = 0 (51) 
mit der positiven Wurzel 
u? = 0,166911 . (52) 


Wir priifen jetzt die Giite dieses Verfahrens, indem wir in der Entwicklung (45) noch ein weiteres 
Glied beriicksichtigen und somit also 


moni Be— Be Heat 
Ce =—; 2 =  Eaia eys hie Ee Oa a pe 
2 § 4 60 240-8 : 
ansetzen. Damit kommt 
€,,(0,02) = — 1301,952671 — 0,972306 u* + 0,000007 ut , 
C12(0,4) = — 6,013146 — 0,102740 u? + 0,002917 ut, 


Hy (0,4) = 1 — 0,021120 u? + 0,000133 ut , 
C5,(0.4) = 23,308334 + 0.973333 u? + 0,014583 u1, 
73(0,4) = — 0,108533 u? + 0,001333 ut . 
Die Frequenzgleichung wird 
—1301,952671—0,972306 u?+-0,000007 u* 1 
23,308334-+-1,133686 u2-+-0,00026 Ll u’—0,0000663 u® +0,000006 us —0,108533 u? +-0,001333 ut 


oder explizit 
h(u2) = 23,308334 — 140,171146 u? + 1,630237 ut + 0,000663 us +- 0,000007 uw’ = 0. (54) 


Der Vergleich mit (51) zeigt, daf§ u? (52) schon eine Naherung fiir ihre kleinste positive Wurzel 
u2 darstellt. Daher setzen wir 


Dic git : : ? 5 
u2 == u? -+ h2 mit h? < v3 (55) 


und haben nach dem Newtonschen N daherungsverfahren 


2 h(u*) 7a 0,0408 ; 
h? = A 5(u2)] = 139,6269 0,0003 . (56) 
d(u*) j 


Aus der Kleinheit der Korrektur h? (56) schlieBen wir auf eine befriedigende Konvergenz des 
Naherungsverfahrens. Wir begniigen uns mit der so erhaltenen Naiherung 


uz = 0,1666 (57) 
zum niedersten Eigenwert u? im Falle x,=0,05. Die Durchfihrung dieses Verfahrens fiir verin- 


derlichen Wert des Parameters x, ergibt die Werte von Tabelle 1 fiir den niedersten Eigenwert 


(Tabelle 1): 
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Tabelle 1. Niederster Eigenwert fiir p = 1, ¢ = 2,5 ra; Scheibe am Innenrand starr 
eingespannt, Aufpenrand fret 


u? ue 


of 
* * eae, We 
(2. Naherung) (3. Naherung) OW 


Der Tabelle entnimmt man, da die Konvergenz des Verfahrens mit zunehmendem Halbmesser- 
verhaltnis r,/r, schlechter wird. Daher sind die Eigenwerte fir r;/r, = 0,20 und 0,25 noch in einer 
dritten Naherung berechnet, die man dadurch erhalt, daB man in den Entwicklungen fiir 7,, und 
Cy. je ein weiteres Glied beriicksichtigt. Das Ergebnis ist, auf die bezogene Higenfrequenz 


2H OT, /v/gG umgerechnet, in Abb. 1 dargestellt. - 


5. Beispiel II. Das zweite Beispiel unterscheidet sich von dem soeben behandelten Beispiel I 
lediglich dadurch, da® wir hier in der Profilkurve (6) p = 2 wahlen, also von der Differential- 
gleichung 


al | 


ausgehen, wozu fiir die Koeffizienten a2,, der Partikularlésung die Rekursionsformel 
[n (n— 1) + 3n] a2,y —2 (n — 2) ag, n-2 + Wag, np 2 = 0 (59) 


gehért. Daraus folgt mit aj, = 1, ag, = 0 


U2 


Yer =1 8 e 


wi(4—u"*) 3 
B24" 2 


(4 — u?*) (8 —u?) u? 
8-24-48 


(4 — u?) (8 — u?) (12 —u?) u? 
8-24-48 - 80 i 


&6 


go... (60) 


Die in der zweiten Partikularlésung auftretenden Koeffizienten b,,, werden wie oben berechnet 
und ergeben sich zu 


by = Fy (v?, p)=9, 


boy = Fy9(u?, p) = 4 » b= F,,(u*, p)=0, 


P Tut + 56 u2 + 96 
bog = Fy,(u*, p) 48 ai : tg F,,(u*, ele 

61) 
19 uS + 228 u* + 800 u2+ 768 ( 
bog = Fyg(u*, p) = i= 6 - 48 mt : 


boy = F,,(u’, p)=9, 


__ 298 u® + 4768 u® + 25568 ut + 51968 u? + 30720 
48 - 24+ 80-8 ; 


bog = Fy,(u?, P) 


Die weitere numerische Auswertung liefert schlieBlich die in Tabelle 2 und Abb. 1 angegebenen 
Naherungswerte fiir den niedersten Eigenwert bzw. den bezogenen Wert der entsprechenden 
Kigenfrequenz. Die Konvergenz ist fiir technische Bediirfnisse noch ausreichend. 


Tabelle 2. Niederster Eigenwert fiir p = 2,c = 2,5 ra; Scheibe am Innenrand starr 
eingespannt, Aufenrand frei 


u2 


(1. Naherung) 


2 
uu“ 
* 


(2. Naherung) 


2 
** 


(3. Naherung) 


u 


i 
2m 
Be V; G 


0,05 0,1453 0,1398 0,1397 0,1495 
0,10 0,5857 0,5679 0,5675 0,3013 
0,15 1,3343 1,3123 Ea 0,4582 
0,20 2,4085 2,4226 a 0,6226 
0,25 3,8208 3,9718 3,9972 0,7997 


XXVI Band 1958S. Tameroglu: Torsionsschwingungen von Scheiben mit Exponentialprofil 219 


6. Zusammenfassung. Um die Eigenfrequenzen der Torsionsschwingungen von nichtrotieren- 
den Scheiben zu ermitteln, kann man, ausgehend von der Differentialgleichung (4) fiir die Schwin- 
gungsamplitude, bei vorgegebenem Scheibenprofil mit Hilfe von Reihenentwicklungen die Fre- 
quenzgleichung fiir die ersten Eigenfrequenzen in Form einer algebraischen Gleichung erhalten. 
Bei Scheiben mit einem Exponentialprofil (6), deren Innenrand starr eingespannt und deren AuBen- 
rand frei ist, konvergiert dieses Verfahren hinsichtlich der niedersten Kigenfrequenz befriedigend 


LUOCTY a 
09 96 = 


0 G05 G10 GT5 G20 G25 


Abb. 1. Niederste Torsionseigenfrequenz von Scheiben mit starr eingespanntem 
Innenrand, freiem AuBenrand und dem Profil 


y=a ee Ta)? fur p = 1 und p= 2. 

und liefert die in den Tabellen 1 und 2 und in Abb. 1 dargestellten Ergebnisse. Da die Dif- 
ferentialgleichung (5) fiir die Amplitude der Dehnungseigenschwingungen im wesentlichen von 
derselben Art ist wie (4), laBt sich das Verfahren grundsatzlich auch auf die Bestimmung der 
Kigenfrequenzen fiir die Dehnungsschwingungen anwenden. Es zeigt dann aber eine mangelhafte 
Konvergenz — eine Beobachtung, die man auch in anderen Fallen schon gemacht hat! —, so daB 
es sich hier fiir die numerische Rechnung nicht empfiehlt. 


(Eingegangen am 11. September 1957.) 


Anschrift des Verfassers: Dr.-Ing. S. Tameroglu, Istanbul (Turkei) 
Nigantas, Giizelbahce, Yeni Kardesler Sok. No 1/7. 


1 Vgl. L. Collaiz, Eigenwertaufgaben mit technischen Anwendungen, S. 396f. Leipzig 1949. 
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Erganzungen zur Theorie des Voith-Schneider-Propellers 
Von W.-H. Isay 


1. Einleitung. In drei vorangehenden Arbeiten! wurde eine Theorie zur Behandlung der 
Strémung durch Voith-Schneider-Propeller entwickelt; diese umfaBt sowohl freifahrende Propeller 
als auch aloes im Nachstrom eines Schiffsrumpfes. In der vorliegenden Arbeit sollen zur weiteren 
Abrundung der Theorie noch einige Fragen untersucht werden. Diese lassen sich im wesentlichen 
mit den bine entwickelten Methoden iohaudele: so dafB wir uns auf eine knapp gefaBte Darstellung 
beschrinken kénnen. 


2. Zwei Voith-Schneider-Propeller nebeneinander. Viele Schiffe haben zwei nebeneinander an- 
geordnete Voith-Schneider-Propeller (Abb. 1), und es diirfte interessant sein, festzustellen, inwie- 
weit sich diese gegenseitig beeinflussen, und ob es wenigstens naherungsweise méglich ist, beide 
Propeller als Einzelpropeller zu behandeln. 

Um den Einflu8 des Nachbarpropellers auf den betrachteten Propeller zu ermitteln, ersetzen 
wir ersteren durch eine 4quivalente Zirkulationsdichte (vgl. Arbeit II, Ziff. 2). Somit ist die vom 
Xo Propeller 2 und seinen freien Wirbeln in der Kreisbahn des Propellers 1 


fe 2 induzierte Geschwindigkeit gegeben durch 
\ pe \ - 2, ae do 
j A ae A) Bere pA) ek poem 1 (6 
\ DS 2 i uy Uv, On 3 I, ‘ele 
lg : 0 
—_— \o 220 3 TNS), ef? are 
aie: | | an | | 
f gs omit y=4/R und 0)>2. Analog erhalt man die vom Propeller 1 in der 
\ | i \ Kreisbahn des Propellers 2 induzierte Geschwindigkeit u®, v®, indem 
me a man in (1) 6) durch (—d,)) und J’, durch J) ersetzt. Formel (1) gilt 
Bae fiir einen 4fligeligen Propeller; bei 6 Fliigeln ist lediglich der Faktor 2/3 


wegzulassen. Die Randbedingung an den Fliigeln der beiden Propeller 
ist natiirlich jeweils die gleiche wie in Arbeit I, Ziff. 4; lediglich sind zusatzlich die vom Nachbar- 
propeller stammenden Geschwindigkeitskomponenten u®, v® bzw. u®, v® zu beriicksichtigen. 
Man erhalt dann ein System von zwei simultanen Integralgleichungen fiir die Zirkulations- 
verteilungen y,(y,t) und y,(y,t). Die aus der Randbedingung am (4-fliigeligen) Propeller 1 ent- 
stehende Cbtaaa ddl lautet [vgl. (II, 20)] 


wr, cos (wt + A) — uy cos (p + wt) + uy 3 cos (wt + A) sin (py + t) 
n 2a 
ah bilge [ { r ae sin ® + sin (® — #) — dy cos © + A(t) [xy cos B + dy sin B + cos (®— #) — 1] 
Anw 3 2). Seah Ce i ew SNES: 5) 
Mhig 0 506 a c08 (D0) — 7 ton Pte y 0s = Og sin Es Oy ata 
x di dy 


22 
1 2 [( ,, {#\sin (® — #&) — 0) cos ® + A(t) [59 sin ® + cos (PW — 9) — 1] 
ee ere 3 Spi 2 te eee eal 
uw 3 a) 
0 


1 
We 3 0 — cos (® — #) — oy sin ® + 6, sin 3 


4 a dy 


aff (2 sin D + sin(® — 9) + h(L) [x cos ® + cos (b— 9) — 
1+ a4! cos (® — #0) —x cos B+ xcos 


4 
gna ee ; q—l1a 
t mR OOS (wt + A) : re +- ee a 
q 


4 
l 1 I —1 
=;)> [afeet?s : oy te 3 (# y — a) dy 
q=1 


I 


‘ : EL Ing.-Arch. 23 (1955) S. 379 (Arbeit I); 24 (1956) S. 148 (Arbeit IL); 25 (1957) S. 303 
Arbeit 
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Dabei ist 


P=ypt+ot, hit)= ne cos(wt+A), I(t) = swsi@n t) dy. 


Die aus der Randbedingung am Propeller 2 entstehende Integralgleichung (2’) ist entsprechend, 
nur ist y, durch y, und 6) durch (—6o) zu ersetzen, so daB wir darauf verzichten, sie extra aufzu- 
schreiben. Die bei einem 6-fliigeligen Propeller geltenden Integralgleichungen lassen sich unter 
Beachtung der friiheren Gleichung (I, 6) sofort ablesen. 

Fiir die Auflésung des simultanen Integralgleichungssystems (2), (2’) wird das bisher benutzte 
und in Arbeit I ausfiihrlich beschriebene Iterationsverfahren genau entsprechend angewendet. 
Die Kerne der Integralgleichungen (2), (2') sind ja gegeniiber (II, 20) baw. (I, 6) unverandert 
geblieben; die neu hinzugekommenen Glieder stehen alle auf der linken Seite. 

Bei der Durchfiihrung des Iterationsverfahrens verwendet man fiir die neu hinzugekommenen 
Oa die nachstehend aufgefiihrten Integralformeln, deren Beweis wir in Ziff. 5 nachtragen: 


z{ ic pine % sin @ =e sin Sates cos ® + h(t) [x cos Abie Oo sin ® + cos (QO) = Z| di) "| 


| 3 Vaal i 62 — cos (®B — #) — xy cos © + yx cos & — by sin B + 6, sin J 


: ag A E lS ee —i®d\» 
etter), a ee Se = ay eal =| (ale, eae | 
n (eae epee n ou a 109 j 
27 
Ng { pin sin (® — 9) —d4 cos ® + h(t) [dp sin ® + cos (P—9)—1] yo | 
ge 1 + 503 — cos (@ — 8) —6, sin ® + dy sin # L (4) 
— (i-- A) e—*? — i + h(t) "\n+1 ,—i® < ae eat! a > G 
ate 16)" onel (=1) A 2 v 1 Oo (mel, d9>2), 
22 5 ; 
1 i sin (P — oe cos PD + h(t) [sin ® + cos (®—8)—1] 4 G-phye*? G—hye® 
xD . ! id r * 
e. 1 + 4.53 — cos (@ — 8) — by sin @ + dy sin 8 ¢ + 1d, oe? — id, 


(5) 
Die Berechnung der resultierenden Fliigelkrafte erfolgt genau nach dem friitheren Verfahren 
(Arbeit I, Ziff. 8), wobei natiirlich die vom Nachbarpropeller induzierten Geschwindigkeitsanteile 
u®, vy) baw. u®), v), zusatzlich zu beriicksichtigen sind. Zur Berechnung dieser zusatzlichen 
Geschwindigkeitsanteile benétigt man die nachfolgend angegebenen Integralformeln, fiir deren 
Beweis auf Ziff. 5 verwiesen wird: 


co 2a 


l [ i e'"” (sin & — sin © + d)) dd dy i nor ae 
c Sr Tee 
ake oe 42 + 4.03 — cos (P — 0) —7.c08 ® + 4.08 9 —d, sin ® + 6, eo, ne a 
as ein? 
1 (cos ® — cos ea ONY en 255 1 
va | | ~ 63 D—# D Vc ta ante "Gy tl, 
gece =H a: cos ( )—yx cos B + x cos & — dy sin ie (7) 
20 
ein? . 
Z | : (sin # — sin B + 6g) dd eee Bees 7 oe ae (8) 
3 1+ = 9 — cos (P — #) — 6, sin B® A Se (e + #65) 
22 
ind 
3 | 1 se PT eee ae 5 yet l (n> 1), (9) 
oe as 62 — cos (BD — #) — dy sin D + dy sin 0 (e + ¥0o) 
22 , 
1 | (sin ® — sin DB + 0,) dd gue 2 Oo asin © (10) 
2: T 1+ 62—26)sin ®’ 
7 14 ; 62 — cos (® — #) — 6, sin B® + oy sin # r 6 rt 
1 ie (cos B — cos #) dd i 2cosD (11) 
eal 1 5 Don ae ge 2d, sin D 
wa 5 62 — cos (D — #) — 6y sin ® + dy sin B 


15 
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Auch zwei nebeneinanderliegende Voith-Schneider-Propeller kénnen ohne weiteres im Nach- 
strom eines Schiffsrumpfes liegend behandelt werden (Abb. 2). Entsprechend der in Arbeit IIT 
Ziff. 2 entwickelten Theorie haben wir dazu lediglich die homogene Anstrémung uy zu ersetzen 


durch 


up + aka | : R) 


zta+b 1005 


Fiir die Randbedingung am Propellerfliigel benétigt man die Werte von u* und vj auf dem Umfang 
der beiden Propellerkreise. Es ist 


io 2 21@ 
Hobe : ie 2Re 3 Re (12) 
poi tan 


* ee ee 
L bap = 
Uy UUs j 


1 2 
a+b—— dR (a +5 —+ 4, R) 


auf den Umfang des Propellerkreises 1; und entsprechend auf dem Kreis 2, wenn 6) durch (—d,) 
ersetzt wird. Die Zerlegung von (12) in Real- und Imaginarteil ist ohne weiteres méglich, und die 


~—a 
- b a Ae © 
aie an 
/ \ 
4 ues 
\ | 1 ») 
lie? 


Abb, 2. 


auf der linken Seite des Integralgleichungssystems (2), (2) zusatzlich auftretenden Anteile sind 
dann leicht angebbar [vgl. Formel (III, 4)]. 


Die tiber den Propellerumfang gemittelte Nachstromziffer betragt nach (12) fiir jeden der 
beiden Propeller in x-Richtung 


(a 4b) = 58 R 
PA Reabe i 5 (13) 
te + b+ 83 Ri 
Der iiber den Propellerumfang gemittelte Nachstrom in y-Richtung betragt 
+ uy B 5p R (a + 6) i bei Propeller 5 (14) 


C + b+ 4 8 rel 


3. Zahlenbeispiele. Wir teilen hier kurz die Ergebnisse von zwei Zahlenbeispielen mit, die mit 
Hilfe der Theorie aus Ziff. 2 berechnet wurden. Dabei handelt es sich um einen Afligeligen Voith- 
Schneider-Doppelpropeller mit folgenden Daten: Fliigeltiefe 2x = 20°, Steigung r}/R = 0,5, Fort- 
ee up/o R = 0,25, Fligelstellung 2 = 0 (fiir beide Propeller). Ferner sei der Abstand 

es oe 

Es sei dann Beispiel 6 der eben beschriebene Doppelpropeller freifahrend, und es sei 
Beispiel 6a der gleiche Propeller im Nachstrom des durch b—=3R, a=2R gegebenen Schiffs- 


rumpfes. 
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Alle Geschwindigkeiten und die Zirkulation sind wieder in Einheiten von Uy angegeben; alle 
Krafte sind in Einheiten von 0 u2R angegeben!. Von den einzelnen Geschwindigkeitsanteilen, aus 
denen sich die resultierende Geschwindigkeit w (relativ zum Propellerfliigel) zusammensetzt, 
geben wir hier lediglich die gegeniber friiher neuartigen ul), v), ul), v(?) sowie die Nachstrom- 
komponenten u*, v* wieder. 


Tabelle 1. (Beispiel 6) ’ 


Propeller 1 Propeller 2 
wt | u@) | vc) | wt vee amt | u@) v2) 
( + 0,000 | + 0,043 | 180° | + 0,014 | + 0,042 0° |—0,039 |—0,084 + 0,042 |—0,092 
S08 + 0,002 | + 0,046 | 210° | + 0,012 | + 0,039 30° |—0,024 |—0,076 + 0,055 |—0,118 
60° +0,006 | + 0,047 | 240° | + 0,007 | + 0,037 60° |—0,011 |—0,073 + 0,045 |—0,163 
90° +0,006 | + 0,046 | 270° | + 0,006 | + 0,036 90° |+ 0,001 |—0,072 —0,017 |—0,187 
120° + 0,008 | +0,048 | 300° | + 0,004 | +0,037 | 120° |+ 0,014 |—0,074 — 0,062 |—0,143 
150° | + 0,013 |.-+ 0,047 | 330° | + 0,001 | + 0,039 | 150° |+ 0,027 |—0,080 —0,056 |—0,102 
Tabelle 2. (Beispiel 6a) 
Propeller 1 Propeller 2 
wt | ugh vo uk oe ot uf?) v@) uf vi 
0° + 0,031 | + 0,057 | —0,208 | + 0,110 0° — 0,030 | —0,108 | —0,208 | —0,110 
30° + 0,037 | + 0,062 | —0,241 | + 0,086 30° — 0,014 | —0,097 | —0,185 | —0,144 
60° + 0,046 | + 0,061 | —0,284 | + 0,067 60° + 0,002 | —0,092 | —0,175 | —0,184 
90° + 0,048 | + 0,055 | —0,350 | + 0,068 OOF + 0,017 | —0,090 | —0,174 | —0,228 
120° + 0,050 | + 0,053 | —0,422 | + 0,122 120° + 0,032 | —0,091 | —0,195 | —0,287 
150° + 0,051 | + 0,046 | —0,440 | + 0,229 150° + 0,049 | —0,096 | —0,266 | —0,337 
180° + 0,049 | + 0,041 | —0,372 | + 0,322 180° + 0,068 | —0,109 | —0,372 | —0,322 
210° 0,043 | + 0,040 | —0,266 | + 0,337 210° + 0,086 | —0,140 | —0,440 | —0,229 
240° + 0,038 | + 0,040 | —0,195 + 0,287 240° + 0,076 | —0,197 | —0,422 | —0,122 
270° + 0,034 | + 0,042 | —0,174 | + 0,228 270° + 0,003 | —0,228 | —0,350 | —0,068 
300° + 0,031 | + 0,045 | —0,175 | + 0,184 300° — 0,051 | —0,178 | —0,284 | —0,067 
330° + 0,030 | + 0,051 | —0,185 | + 0,144 330° —0,047 | —0,131 | —0,241 | — 0,086 


In den Tab. 3 und 4 ist die resultierende Geschwindigkeit w nach GréBe und Richtung (gegen 
die positive x-Achse) enthalten. Ferner die Profilzirkulation /’, und die Fliigelkrafte K zerlegt in 
ihre x- und y-Komponente. Dabei haben wir wie in Arbeit II S. 158 angenommen, daf} der Druck- 
mittelpunkt mit dem Profilmittelpunkt zusammenfalle, und da es sich um abgerundete symme- 
trische Profile? handelt. 


Tabelle. 3 (Beispiel 6) 


Mittelpunkt Propeller 1 Propeller 2 


des 
Fliigelprofils 


DY; Ky Kyy a Bs I", | Kx Koy 

0° 4,29 |— 67,0°, —0,32 —1,24 | —0,53 | 4,35 ;— 681°; —0,41 | —1,64 — 0,66 
30° o,l1 | — 42,9°| —0,49 —1,71 —1,84 | 5,12 |— 44,0°| —0,61 — 2,19 — 2,206 
60° 5,44. | — 23,8°| —0,62 —1,37 — 3,10) | 5,43 | — 25,0°) —0,75 | —1,73 — 3,10 
905 4,95 |— 3,9°| —0,23 — 0,08 —1,13 | 4,93 /— 5,2°} —0,35 — 0,16 —1,73 
120° 4,92 | +. 19,7° | 0,46 — 0,76 +2,14 14,85 |-+ 184°) + 0,32 — 0,48 + 1,45 
150° 4,69 | + 44,8° + 0,93 —3,06 | + 3,09 | 4,58)-+ 43,5°) + 0,77 — 2,43 + 2,56 
180° 4,25 | + 72,3°| + 0,88 — 3,94 +1,13 |4,10|-+ 71,2°| + 0,74 —2,88 | + 0,98 
210° 3,75 | + 103,3°) -+- 0,54 | —1,98 —0,47 | 3,52 | + 103,1°| + 0,48 — 1,66 — 0,38 
240° 3,30 | + 138,7°| + 0,11 — 0,24 —0,27 | 3,04 |+141,1° | + 0,22 —0,42 | —0,53 
270° 3,11-| + 177,8°| —0,23 | +,0,03 + 0,71 | 3,01 | —178,0°}) -+ 0,02 0,00 — 0,07 
300° 3,03: | —139,1°| —0,32 — 0,64 + 0,74 | 3,09 | —137,5°| —0,20 —0,41 | + 0,45 
330° 3,47 |— 98,9°| —0,24 | —0,82 | +0,13 | 3,56 99,4° | —0,28 — 0:97 + 0,16 


1 Die Durchfiihrung der in der vorliegenden Arbeit enthaltenen numerischen Rechnungen lag in den 
Hiinden von Frau Marzke und Fraulein Endemann (beide Rechner bei der Deutschen Akademie der Wissen- 
schaften zu Berlin). : é ; f 

2 Dazu mag hier noch ausdriicklich bemerkt werden, daf ja die Gesamtzirkulation I’ eines diinnen un- 
gewolbten Profils in guter Naherung mit derjenigen eines dicken Profils gleicher Skelettlinie (Stromungs- 
bedingung) iibereinstimmt. Das laBt sich ohne weiteres mit Hilfe der von H. Schlichting im VDI-Forschungs- 
heft 447 (1955) angegebenen Formeln zur Profiltheorie zeigen und kann in analoger Weise auch fir die Profile 
des Voith-Schneider-Propellers nachgewiesen werden. 
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Tabelle 4 (Beispiel 6a) 


PORE ees Propeller 1 Propeller 2 
Fiigelprotds Tonal gt Sox K2y 
0° 4,25 | — 65,8°| —0,22 — 0,84 | —0,38 | 4,33 |— 69,7°| —0,55 — 2,25 | —0,83 
30° 5,04 | — 42,5°) —0,44 —1,49 | —1,62 | 5,03 |— 46,8°| —0,87 — 3,18 — 2,98 
60° 5,26 |"— 24,4°| —0,65 —1,40 | —3,09 | 5,27|— 27,8°| —1,00 —2,47 | —4,68 
90° 4,65 |— 4,2°| —0,25 —0,09 | —1,17 ]4,73|— 85°} —0,63 — 0,44 — 2,93 
120° 4,60} -+ 21,6°} + 0,56 — 0,94 + 2,38 | 4,52}+ 15,0°; —0,02 + 0,02 — 0,08 
150° 4,59 |}+ 48,9°) + 1,27 — 4,40 + 3,84 | 4,05) + 41,4°] -+.0,47 — 1,25 + 1,42 
180° 4,53 | -+ 76,3°) -+ 1,34 — 5,89 + 1,44 |3,51|-+ 73,2°) + 0,72 — 2,41 -+- 0,73 
210° 4,24 | + 104,4°| -+ 0,77 — 3,15 —0,81 | 3,23 | + 108,7°| + 0,80 — 2,43 — 0,82 
240° 3,76 | + 136,6°| + 0,02 — 0,04 —0,05 | 3,14 | + 145,4°| + 0,53 — 0,95 — 1,37 
270° 3,43 | + 173,6°} —0,50 | + 0,19 + 1,71 | 3,25 | —177,1°| + 0,11 + 0,02 — 0,36 
300° 3,06 | —143,0° | —0,58 — 1,06 + 1,41 | 3,17 |; —136,5°; —0,14 — 0,31 + 0,33 
330° 3,41 | — 99,3°| —0,30 | —1,02 + 0,17 | 3,60 |— 98,9°| —0,23 | —0,81 + 0,13 
r | 
ae ee 3 1 | | 
7 a 
06 i. | 
o# | i | 
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02 + + | 
0 be 
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Abb. 3. Beispiel 6. Zirkulationsverteilung 1’, (Propeller 1) und J’, (Propeller 2). 
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Abb. 4. Beispiel 6a. Zirkulationsverteilung I’, (Propeller 1) und J, (Propeller 2). 
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In Abb. 3 und 4 sind die Zirkulationsverteilungen der behandelten Beispiele aufgezeichnet. 
In Tab. 5 sind die resultierenden Gesamtkrifte der behandelten Voith-Schneider-Propeller zu- 
sammengestellt. (Als Mittelwerte der drei aus den Tab. 3 und 4 berechenbaren Propellerstellungen 


0, 2/6 und 7/3.) 


Tabelle 5 


Zusammenfassung 


Beispiel Propeller K(P) 
x 
2. K(P) 


Richtung BE 
| Richtung px 


1 et + 0,20 BLA 4.177,8° : 

6 5 4299 pg CY 5.14 aT a 10,18 al tat 
1 6:71 + 138 6.83 + 169.2° ’ 

be 2 Se ag | ean 6,08 Sede og leer ease 


Die Ergebnisse von Beispiel 6 zeigen, daB zwar die zusammengefaBte Gesamtkraft beider 
Propeller das Doppelte der Kraft des entsprechenden Einzelpropellers (Beispiel 5) betrigt. Jedoch 
bestehen in den Einzelheiten des Strémungsbildes der beiden nebeneinanderliegenden Propeller 
immerhin gewisse Unterschiede. 

Diese Unterschiede sind natiirlich im Fall des Doppelnachstrompropellers Beispiel 6a wesent- 
lich ausgepragter; denn hier sind ja die Anstrémverhaltnisse beider Propeller ganz verschieden. 
Die Nachstromziffer betragt nach (13) A=0,276, und der mittlere Nachstrom in y-Richtung 
gemaB (14) + 0,182. 


4, Seitwartsfahrt. Fiir zwei nebeneinanderliegende Voith-Schneider-Propeller ist auch die Seit- 
warts- und Schragfahrt von Interesse. Wir wollen auf dieses Problem in seiner vollen Allgemein- 
heit mit der Vielzahl von Parametern nicht eingehen. Lediglich die Behandlung des wichtigen 
Falles der Seitwartsfahrt sei kurz angedeutet. 

Dabei mégen die beiden Propeller wie in Abb. 1 liegen; sie werden in Richtung der positiven 
y-Achse mit der Geschwindigkeit v, angestrémt, bzw. sie bewegen sich in Richtung der negativen 
y-Achse. 

Es sei ferner A = —a/2, d.h. also —r'/R = h(t) =1j/Rsinwt fiir die Fliigel beider Pro- 

eller. 
j Neuartig ist jetzt, daB die von den Fliigeln des Propellers 1 abgehenden freien Wirbel durch 
den Propeller 2 hindurch flieBen. Das Geschwindigkeitsfeld der freien Wirbel des Propellers 1 


ist in einem beliebigen Punkt z= re‘? gegeben durch 


co 2a 
i 2 as dd dy 
ia he —~ Sone 3 [ felon epee: (15) 
0 0 
wahrend man fiir den Propeller 2 hat 
o 2a d 
a 8 dd dy 
; —— — aed = = 16 
ee de 200 3 PP ls he?’ sy =, RB (16) 
0 0 


(4fliigeliger Propeller). Die Randbedingung an den Propellerfliigeln ist prinzipiell genau wie friher, 
und der Rechnungsgang an sich auch. Jedoch sind die durch (15) und (16) gegebenen Geschwindig- 
keiten teilweise neuartig, so daB wir kurz die aus ihnen entstehenden Ausdriicke angeben wollen. 
Die Geschwindigkeitskomponenten der freien Wirbel beider Propeller geben auf der linken Seite 
der aus der Randbedingung an den Fliigeln von Propeller 1 entstehenden Integralgleichung mit 


y= Ry den Ausdruck 


o ln 
: ey | | r(2) sin eave cos D + h(t) [x sin ® + cos (P == 0) —1] dd dy 
moet ' F 7 Tos x" — cos (P — 8) — x sin D + x sin O 


+ 


47 a) 


oO 2% 
1 eo | | T ( 4 SE Dieses (9 rao) CO8'D HY) Uiete 0) si: Pres (P= 0) = Ug gy 
0 0 
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Der entsprechende Ausdruck (17') auf der linken Seite der aus der Randbedingung an den Fliigeln 
des Propeller 2 entstehenden Integralgleichung wird formal aus (17) erhalten durch Vertauschung 
von /', mit J’, und 6, mit (— 6,). Der Charakter der beiden Ausdriicke (17) und (17’) ist jedoch 
vollig verschieden. Dieses ist am besten aus den folgenden Integralformeln ersichtlich, die fiir die 
Auflésung des Integralgleichungssystems nach dem friiheren Iterationsverfahren benétigt werden. 
Es gilt mit 6)>2 


ae [ p Hn? — 9) = 8) gon Pe MAL a Oo) sin Ds 08 (Bt) aaa 
2 oo ‘ iL 4k = (x + 5)? — cos (Sb — #) —(y + 6o) sin D + (x + 69) sin B 


(18) 
A =D 
+ ip (Tiga) < (0=G=2n), 
Ni (Ea cen Oye 
1 [fe p sin (DB) —(%—A,) 08 ® + HC) [Cz —A4) sin ® + c08(P—H)—N gy gy 
ae Toe pee ye coe (Dae ee (y —6,) sin 8 
ea : (19) 
i le hee : : @ 
= [cos (n — 1) b—cos(n+ 1) @ +hsin (n —1) © + hsin (n +1) P] 
n (e~ ee Jy)” eee (3 
+fir 0S @=< 2; —fir 2x @O=22, 
al 
08 7) ji 
06 
7 
ea, 
OY 1 IL 
02 t r 
0 ie =F = = 
a 
if 
-O4 7 aah i Ts 
ot 
06 2a ba es 7 «On Int 5 ia 0 4 ed ra 
3 6 6 3 2 3 o 6 3 Z 
Abb. 5. Beispiel 7. Zirkulationsverteilung I, (Propeller 1) und IP, (Propeller 2). 
i ping in (@ — sehet cos D + h(t) [y sin © + cos (SG — 9) — 1] do dy 
CON 1 + + %— cos (@—#)—y sin ® + x sind 
1 (20) 


last h) ees De fun gs OS z., 


(1 + ih) ei ® _ = og (n +1) 0+ h= sin (n +1)0 fir 02027. 


Fiir den Beweis der Integralformeln (18), (19), (20) wird auf Ziff. 5 verwiesen. 

Nach der Theorie von Ziff. 4 wurden die Zirkulationsverteilungen )(t) und [°,(t) der Fliigel 
eines seitwartsfahrenden Voith-Schneider-Doppelpropellers mit den folgenden Daten berechnet 
(Beispiel 7): Je 4 Fliigel, Fliigeltiefe 2 = 20°, Steigung ro/R=0,25, Fortschrittsgrad v,fOoR=Ye 
Abstand 6)=3. Das Ergebnis der Rechnung ist in Abb. 5 dargestellt. Man erkennt deutlich, 
wie die Fliigelzirkulation und damit auch die Fligelkrafte bei Propeller 2 herabgemindert werden 
infolge des Durchflusses der von Propeller 1 stammenden freien Wirbel. 
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5. Einige wichtige Integralformeln. Den Beweis der in Ziff. 2 und 4 mitgeteilten Integral- 
formeln fiihren wir in einer verallgemeinerten Form, indem wir nunmehr auch freie Wirbel zu- 
lassen, deren AbfluBrichtung mit der x-Achse einen beliebigen Winkel @ einschlieBt. In Ziff. 2 
war © = 0, und in Ziff. 4 war O = 7/2. 


a) Die freien Wirbel des gerade betrachteten Propellers (EKigenpropellers) geben in der Rand- 
bedingungsintegralgleichung AnlaB zu dem Integralausdruck 


o 27 


ein? oe | | ein(9— 0) sin (DP — #) nt x sin (f= 0) + h(t) [x cos (© — O) + cos (SG — 9) — 
22 1 aici seb akiod a 
noe 1 + > x* — cos (G — 4) — x cos (S — O) + x cos (8 — @) 


zl dd) dy 


tee eee oh: , 
= (i a h) e(n—1)® gO fay ue =o= 6 < Aish 


— (i me! h) e—t(n—]) P e! Qn—1)0__ ft (i ae h) ei (n | 1) @ et 
nm 


—l)j.. : é 
se : ) (i ae h) e—i(tn+ 1) 2 ei (2n +1)0 fiir 3 ea @D a) = 4 : 


Die Richtigkeit dieser Integralformel entnimmt man unmittelbar aus dem Ergebnis von Integral- 
formel (I, 9). Damit ist als Sonderfall O = 7/2 auch Formel (20) bewiesen. 


b) Die freien Wirbel des Nachbarpropellers geben in der Randbedingungsintegralgleichung 
fiir n = 1 AnlaB zu dem Integralausdruck 


sin (® — 9) + x sin (® — 9) — dy cos ® + A(t) [% cos (P — O) + 6p sin B + cos (P — B) — 1] dd dy 
1-4 i | 52 cos (P — #) — x cos (®W — @) + x cos (A— O) + 05 x sin O — 6, sin ® + dy sin O 


Wir berechnen das Integral I, unter der Voraussetzung, daf die freien Wirbel des Nachbarpropellers 
nicht durch den betrachteten Propeller hindurchflieBen. Wie man sich durch Aufzeichnen der 
entsprechenden komplexen Vektoren in der (x, y)-Ebene leicht klar macht, bedeutet diese Voraus- 
setzung, da immer 

lid) Hy e?—e®|>1 
sein muB. In dem Integral [, vollziehen wir zunachst die Integration tiber # nach der Residuen- 
methode mit der Substitution z=e'’; Integrationsweg wird der Einheitskreis. Es folgt 


(ec *o x —e '? — id) (2 — 4) (2 —%) 


x a dy. 


Dabei ist 
1 @ 


a u 
e— iP —i9O 3 g Ga Ne 
—wve + 105 


Zz 
« 


1 


Da also immer |z,| > 1 d.h. |2,| <1 ist, erhalt man weiter nach einigen elementaren Umfor- 


mungen 


t@ 


(i + hen? al Ghent e 
iP, == YF —i®@ sic betel dy eae) a ana 
‘ (e —vre =— 10) VSN (e + 10) 


Damit ist im Sonderfall 9 = 0 auch die Integralformel (3) und als Sonderfall @ = /2 die Inte- 
gralformel (18) bewiesen. AuBerdem ist Formel (4) bewiesen, wenn man bei I, lediglich das Kr- 
gebnis der Integration iiber J betrachtet und y = 0 setzt. 
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Der Beweis der Integralformeln (5), (6), (7), (8), (9), (10), (11) wird ganz ahnlich gefiihrt, so 


da® wir darauf verzichten kénnen, ihn hier wiederzugeben. 


ce) Fiir den Fall, da® die freien Wirbel des Nachbarpropellers durch den betrachteten Propeller 
hindurchflieBen, ist das Integral I, im allgemeinen nur noch sehr umstandlich auswertbar. Hier 
wollen wir uns deshalb lediglich auf den in Ziff. 4 bendtigten Sonderfall O = 7/2, d.h. auf den 
Beweis der Integralformel (19) beschranken. 


Die Integration iiber 7 vollziehen wir wieder mit Hilfe der Residuenmethode. Substitution 
z—e'®, Integrationsweg Einheitskreis. Dann nimmt das Integral iiber # in (19) die Form an 


1 fp ee 22(x dp) cos D +h [22 *? + e® + 22 (y—6)) sin P—22] 1 
2 AER) (—ix—e'? + id) (@ — 4) @ — 2) 


(i + h) Ca 1 
Ey Ge tye ee oe tne re 
—(i—h) e'® (c°® —iy +76,)"7 fur), ized (2, = ls 


fis Jia ieeed (= uy = ei + id), 


wie man nach einer elementaren Zwischenrechnung leicht bestatigt. Die Grenze |z,| = |z,| = 1 
bedeutet 


1 + (y—6,)? > 2(y—6,) sin @ =1, also y=06, oder y—d,=—2sn@. 


Fiir die nun folgende Integration tiber y haben wir zwei Falle zu unterscheiden: 


1. Es ist sin OD > 0, 0< O <q; dann ist |\z,/> 1 fir 0 = y <6), und |z,\ << 1 fie oy 
<6) + 2 sin ®, und |z,| > 1 fir y > 6) +2 sin ®. Also wird das Integral (19) 
59 6) +2 sin (o-) 


FAB) BE Gi F , 
{ G+ hye *? dy (i —h) ef ® (e§® —i y + 155)" dy — 
0 


SG 5 Feed 
(ern ty te)? . Beran, 
do 6) + 2sin D 


Gab h)e dy 
(P+ txts yh 


Ihnen) ie 
= ( UE 6 , Cos (n 1)@ cos (n +1) 84 > asin (n—1)® +=hsin(n +1). 


n (Cm eee Oye 


: ees ist sin @ <0,” < @ < 2a; dann ist |z,| > 1 fir 0 < ¥ <6, + 2 sin @, und |z,| <1 
fir 6) + 2sin®<y< dy) und |z,|>1 fir y~ > dy. Also wird das Integral (19) 
6) + 2sin D On (ee) 
eur hye *? dy ; UD) (pt Dian i §.)2— 
("4 xia) (i—h) e ® (e' ®—iy +. 10)) tdy— | 


0 6) + 2 sin D 5p 


(i+ hye‘? dy 
(ec 2h it Oe 


Lely ican 


n Cae 50 t 


2 2 2 oh oe 
= e081) +--cos (n + 1) ® ——hsin (n— 1) cs sin(n+1)@. 
Fir ®=0 und =z stimmen beide Werte iiberein. Damit ist Formel (19) bewiesen. 


6. Verfeinerte Theorie des Voith-Schneider-Propellers. In der gesamten bisherigen Theorie 
haben wir die Zirkulationsverteilung der Propellerfliigel nicht genau auf der Kontur des Fliigel- 
profils angebracht, sondern auf dem Umfang des Propellerkreises (vgl. Arbeit I, S. 379). Daraus 
ergab sich eine erhebliche Vereinfachung fiir die Kerne der aufzulésenden Integralgleichungen 
Wir wollen jetzt untersuchen, inwieweit sich die bisherige Theorie andert, wenn die Zickeulane ae 
verteilung der Propellerfliigel genau auf der durch 


a q—1 
u (v at Sie + 0) 


ey = Rly Peos(oe +t eta tall 
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gegebenen Konturkurve der Fliigelprofile angebracht wird}. Um das festzustellen, werden wir die 
aus Ziff. 4 der Arbeit I bekannte Randbedingung am Propellerfliigel in der nunmehr verfeinerten 
Form ansetzen und die aus ihr entstehende an die Stelle von (I, 6) tretende Integralgleichung 
diskutieren. 

Zunachst ist es einleuchtend, daB die genannte Verfeinerung nur bei dem gerade betrachteten 
Fligelprofil (Aufpunkt) beriicksichtigt werden mu, wihrend fiir die ibrigen fiinf bzw. drei Pro- 
pellerfliigel in jedem Fall die bisherige Methode beibehalten werden kann. In Formel (I, 1) andert 
sich also nur der Summand q=1. Er lautet genommen in dem Fliigelaufpunkt z =r e'(@ + ©) 


i [ y(y, t) R dy 


em q re@Pt ot) i 


6 ely + wt) 


—X 
mit 
/ 


ges Bit 1p Res (ot +A), 9 = Ril —y feos (ot +2) : 


Den Wurzelausdruck des Bogenelementes haben wir in die Funktion y mit einbezogen, was ohne 
weiteres méglich ist. 

Fur das von den freien Wirbeln induzierte Geschwindigkeitsfeld behalten wir den Ausdruck 
(I, 3), (II, 1) unverandert bei; in der Randbedingung am Propellerfliigel setzen wir wie friiher 
r ~ R in den Formeln (I, 4). Denn bei der summarischen Naherungsformel fiir das Geschwindig- 
keitsfeld der freien Wirbel hat es keinen Zweck, den genauen r-Wert einzusetzen; dadurch wiirde 
die Integraldarstellung nur kompliziert, ohne daf§ damit eine echte Genauigkeitssteigerung ver- 
bunden ware. Die (I, 5) entsprechende Randbedingung am Propellerfliigel lautet somit 


worr’ + uyrcos (yp + wt) + ugr’ sin(p + @ t) 
= — [4 (R,¢g +t) +», (7,9 + w 1)] [rsin (p + wt) —r’ cos (p + wd)] (21) 
—[u,(R, p+ ot) +u,(7,9 + @2)][reos(p+o@t) +r’ sin(y + wt)] 
Dabei ist . 
r= R|l—o eos (ot +) 


und wie friiher 
r’ =—r,cos(wt +A). 


Nach Einsetzen der entsprechenden Werte fiir die Geschwindigkeitskomponenten erhalt man 


aus (21) die folgende zu (I, 6) analoge Integralgleichung fiir die Zirkulationsverteilung y(y, t): 


@ r, cos (wt + A) L yp - cos (w t + A) — Uy cos (p + w t) ! —Yy * cos (wt + | 


6 
, 1 ¥ - ly 
+ Uy 7B cos (wt + A) sin (p + wt) 4 repos (wt +A) : ard 
= 


co 22 
/ f 1 ft sin(p + wt) + sin(g + wt— #) 
=a Daeg cos (wt +e | [ a l _% ) = 
eee Leg tos (P 7 Ot — 0) — 7 cos (7 1-01) ey 0980 
x di dy 
1 Ubrigens kann die ganze Theorie ohne weiteres auf allgemeinere Profilformen, z. B. 
r=Rit-+ > gy? ——p z cos (mt -+ A) (« = konst.) 
iibertragen werden. Setzt man dann * =— Z cos(wt-+ A) -+ t¢ in (I, 6) ein, so ergibt sich aus den fruheren 


Formeln, da8 fiir 7 < 4 praktisch nur zu Ay ein Anteil der GroBe — x oa Rv binzutritt. Dadurch wird die 
Gesamtkraft des Propellers nicht geandert; der Schwerpunkt der Kraftwirkung verlagert sich jedoch fur 
t+ > 0 auf den hinteren Teil des Fliigelkreises, und fiir t << 0 auf den vorderen. Die bisherige Profilkontur mit 
+ = 0 ist insofern ausgezeichnet, als sie fast genau die Kriimmung der Kreiskontur hat; sie entspricht somit 
wegen des Uberwiegens der Umfangsgeschwindigkeit in der Anstrémung gegen das Profil dem Fall der Platte in 
homogener, gerader Anstromung. 
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ae d + (22) 
al |? leik = hee FONT ET ee Ot db dy 
x?— cos (p + ot —8)—4 cos (p + wt)—4 cos B 


+“ 0s ( (wt +A) 


ers 


1 ‘ eg ten (pi— 9) — gcos(p—Y) g 
+3q | 7, t) r+ 92 —2recos (vy —/) Ve 
Bevor wir uns der Auflésung der Integralgleichung (22) zuwenden, haben wir ihren (dem 
Summanden q=1 entsprechenden) Hauptkernanteil auf der rechten Seite zu untersuchen und 
umzuformen. Zunachst wird dessen Nenner 


r? + 9? — 2 ro cos (y—yp) = 2 R® [1 — cos (p —y)] [1 + yp) C08 (wt +A) 


y? + y?— 2 py cos (p—Y) (75 2 
as es or en bee cos (ot +2). 


Wegen |p—yi S2x50,42 kann 


1 1 
cos (p —y) = 1— (pv — py? + ag (Y — 
gesetzt werden, so dab 


gy ee 2_ 2m cos (P—Yy) _ a 


2—2cos(p —Y) + PY +7 ne yp)? 


wird. Andererseits hat man fiir den Zabler 


rr +rosin (~—y) —r’ 0 cos (p —yp) = — R* sin (py — y) i + . (p v al cos? (@ t + A) 


— Re" cos (w t + A) ! — cos (p —y) —93 —* cos (wt +A) + ya cos (wt + Aj cos (p —v)| 


} v 2 1 \9 
+ R? sin (p —y) i + yp) F608 (mt +A) 4 l t+ Py + 75 (yp vy) (| cos? (wt +4), 
so daf der betrachtete Kernanteil sich auf die Form bringen 1aBt 


rr’ +resin(y—y)—r’ e cos (p —y) 1 p—y 
=> ctg 


r? + 0? — 2recos(p— yp) 2 2 
1 3) 
, py —y cos (p —) (1 r 19 @ vy) sin (p — ) 
1—~ cos(wt + A) - 
th aK R ee 
ee cos (wt + A) — — ae = 


1—(p+y "8 cos (wt-+A) + (I ee +5 wy) (2 o cost (wtf A) 


Unter Verwendung der ersten Reihenglieder fiir cosinus und sinus erhalt man fiir den Bruch im 


Zahler des zweiten Kernausdruckes den Wert 1/6 p + 5/6 y, so daB sich ergibt 


rr’ +rgosin(~ —y)—r’ 0 cos (py —1 1 ) — 
(vy —y) GCOS DEY re eee 


a7; 
r+ o—2ro cos (p—y) Dela a mee, 2 


cos (wt + A) 


5 1 1 Nella 

ol Tae bree (p wie) eorta 

(Re) cote ba) e + a He & : ; 
L— fp +9) poo (at + Dx ( Font yo (8) Oe 
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Die ersten beiden Glieder sind die aus der bisherigen Theorie bekannten, wahrend das dritte jetzt 
bei der verfeinerten Theorie zusatzlich auftritt. Wie eine nahere Untersuchung mit Zahlenwerten 
bis ro/R < 1 zeigte, kann dieses dritte Glied noch weiter vereinfacht werden, so daB der unter- 
suchte Kernausdruck endgiiltig folgende Gestalt erhalt: 


rr’ +reosin(y —y)—r’ @ cos (p —p) 1 p—yp 1 
2 2 i AN Bis ps aera 
r? + 0? —2rocos(m—yp) Ke) 


—- > Re 008 (mt +A) 
1 
TAS 

eee > 
+3 (3 


Setat man (23) in die Integralgleichung (22) ein, so ergibt sich, wenn wir gleich zu der der Integral- 
gleichung (I, 14) entsprechenden Form iibergehen und den Druckmittelpunkt wy) nach (I, 40) 
einfiihren, 


(23) 


f 
ro \2 


r 
vi 
5 / 
= (3 cos* (wt + A) SP + : y z cos (wt + | 


[2 wr) cos(wt+/)—2 uy cos (py + wt)] jie Z cos (mt +A)| +2 uy A cos (wt +A) sin (p + w t) 


6 6 
(m) ; role 2 (m) (m) , Gare (m) 
++ Ds c, HOt ot) R cos (wt + A) De caer 4 eb Me t mt) 1 3 cos (wt +A) | Ap 


h=—O w=—O6 


ra) aaa (m) . y mY) /p\2 2 IN IPS 1 li 
+A, eoia a py eas o—sio4| oe) i) cos Me = Je ? | ; Wo - cos (w t +A) 
ia) 
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Die Integralgleichung (24) gilt fiir einen 6fliigeligen Propeller; die analoge (II, 21) entsprechende 
Integralgleichung des 4fliigeligen Propellers ist aus (24) sofort zu entnehmen, so daf wir sie nicht 
explizit aufschreiben. 


Vergleicht man (24) mit der entsprechenden friiheren Integralgleichung (I, 14), so sieht man in 
iibersichtlicher Form die Glieder, die durch die Verfeinerung der Theorie neu hinzugekommen 
sind. Sie stehen samtlich auf der linken Seite; durch die vorgenommenen Umformungen haben 
wir erreicht, daB der Kern der Integralgleichung (24) der gleiche ist wie derjenige von (I, 14). 
Ohne weiteres ist klar, daB fiir die Auflésung von (24) genau die gleiche Methode verwendet werden 
kann wie fir (I, 14). Auf eine explizite Weiterverfolgung des Lésungsweges kann daher hier ver- 
zichtet werden. 


Die durch die Verfeinerung der Theorie neu hinzugekommenen Glieder sind unbedeutend, 
solange die Fliigelsteigung nicht zu gro ist. Teilweise kompensieren sich diese Glieder auch noch 
gegenseitig.. So ist im allgemeinen 
2 


(uy —or,) baw. 3 1 © : (u; —@r,) bei 4 Fligeln, 


1 
% 


rr) 


wie sich an unseren bisherigen Beispielen gezeigt hat. 


Das durch den dritten Anteil des Kernausdruckes (23) bedingte Zusatzglied auf der linken 
Seite von (24) hat mit unserem bisherigen a-Wert und r,/R = 1/2 im Durchschnitt die GroBe 
0,006 I(t); es ist also fiir ro/R =1/2 noch praktisch bedeutungslos. 


Im ganzen kann man auf Grund einer iitberschlagigen Abschatzung sagen, daB sich die Zir- 
kulationsverteilung bei den Beispielen der bisherigen Arbeiten fiir r,/R = 1/2 durch die verfeinerte 
Theorie um durchschnittlich etwa 2% bis 3% andern wiirde. Dieses kann in Anbetracht der 
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Naherungsdarstellung des Geschwindigkeitsfeldes der freien Wirbel und der sonstigen in der Be- 
handlung als ebene Strémung liegenden Vernachlassigungen noch als unerheblich angesehen werden. 
Fiir gréBere Fliigelsteigungen rj/R Z?/,; muB fir die Bestimmung der Zirkulationsverteilung 
dann aber doch die Integralgleichung (24) der verfeinerten Theorie zugrunde gelegt werden. In 
diesem Fall ist auch das genaue Bewegungsgesetz der Propellerfliigel etwa (11,10) zu beriick- 
sichtigen (vgl. Arbeit IT, S. 160). 


Aus dem Institut fiir Angewandte Mathematik der Deutschen Akademie der Wissenschaften zu Berlin. 


(Eingegangen am 14, Oktober 1957.) 
Anschrift des Verfassers: Dr.-Ing. habil. Wolfgang-Hermann Isay, Berlin-Dahlem, Schweinfurthstr. 90. 
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Dieses Buch gibt einen Uberblick iiber die Berechnung von Schallvorgingen. 

Es wird zunachst die gebiindelte Strahlung in groBer Entfernung von dem Strahler oder der Strahlergruppe 
behandelt, nachdem die ungebiindelte Strahlung kurz erértert wurde. 

Richt- und Strahlungsfaktor werden definiert. 

Bevor auf das Schallfeld in der Nahe des Strahlers eingegangen wird, ist klargestellt, wann ein Aufpunkt 
zum Nahfeld oder Fernfeld gerechnet werden kann. An Hand verschiedener Beispiele wird das Schallfeld 
in der Umgebung zweier Strahler als Kurven konstanter Amplitude und konstanter Phase dargestellt. 
Das Schallfeld verschiedener Kolbenmembranen wird berechnet und graphisch wiedergegeben. 

Weiter wird eine Reihe von Kugelstrahlern berechnet und die Stérung des Schallfeldes durch eine starre 
Kugel und einen starren Zylinder bestimmt. Die Schallfeldverzerrung um Kugel und Zylinder wird eben- 
falls durch eine Ubersicht von Kurven konstanter Schalldruckamplitude ausgedriickt. 

Alle Formeln sind so angegeben, da es dem Leser jederzeit méglich ist, das Schallfeld fiir einen gewiinschten 


Fall ohne groBe Schwierigkeiten zu berechnen. 
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Dieses Buch bietet eine geschlossene, lehrbuchartige Darstellung der Verfabren der plastischen Berechnung 
biegesteifer Stabwerke aus elastisch-plastischem Material. Es wendet sich an Stahlbau-Statiker, Dozenten 
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Erlauterung des plastischen Versagens einfacher Stabwerke werden die Traglast-Satze angegeben, auf die 
sich die folgend erlauterten Verfahren der plastischen Berechnung fiir einfache proportionale Belastung 
griinden. Die Ermittlung elastisch-plastischer Ausbiegungen wird behandelt und anschlieBend die ver- 
schiedenen Faktoren, die von Einflu8 auf das volle plastische Moment sind. Die Erzielung einer vor- 
geschriebenen Tragfahigkeit bei minimalem Materialaufwand ist egos der Minimalgewichtsbemes- 
sung; die Prinzipe dieser Weiterentwicklung der einfachen Tragfahigkeitsbemessung werden an Hand 
eines geometrischen Analogons veranschaulicht und Lésungsverfahren angegeben. Theorie und Verfahren 
der Traglastberechnung fiir variable wiederholte Belastung werden dargelegt. — Viele Zahlenbeispiele im 
Text machen die Darstellung leicht verstandlich; Ubungsaufgaben mit Lésungen dienen der weiteren 
Vertiefung des Verstandnisses. Umfassende Literaturangaben erméglichen ein weiteres Eindringen in 
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